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MECHANICAL  CONSTITUTIVE  MODELS  FOR  ENGINEERING  MATERIALS 


PART  I:  INTRODUCTION 

Backgrovmd 

1.  Development  of  mechanical  constitutive  models  (defined  as 
load-deformation  or  stress-strain  relationships)  for  engineering  ma- 
terials has  received  considerable  attention  in  recent  years » particu- 
larly in  the  field  of  geotechnical  engineering.  The  primary  reason  for 
such  efforts  is  the  fact  that  with  the  advent  of  high-speed  electronic 
computers  and  the  development  of  new  methods  of  nxanericaQ.  analysis,  a 
variety  of  complex  engineering  problems  can  be  solved  provided  realistic 
constitutive  relationships  for  the  materials  of  interest  are  available. 
Stress-strain  relationships  for  a nximber  of  materials,  such  as  soil, 
rock,  euid  concrete,  are  often  nonlinear  even  when  the  magnitudes 

of  the  strains  involved  are  small.  This  type  of  nonlinear  behavior, 
referred  to  as  physical  nonlinearity,  has  been  the  subject  of  investiga- 
tion at  the  U.  S.  Array  Engineer  Wateivays  Experiment  Station  (WES)  since 
early  I96O;  special  emphasis  has  been  placed  on  modeling  the  mechanical 
behavior  of  eai'th  materials.  During  the  fall  of  1971,  an  elementary 
coui'se  on  raeclanical  constitutive  relationships  was  offered  at  the 
Vicksburg  Graduate  Center,  WES,  and  a series  of  lecture  notes  was  pre- 
pared for  use  by  the  students  taking  this  course.  The  purpose  of  the 
lecture  notes  was  to  acquaint  the  students  with  some  of  the  basic 
physical  concepts  and  mathematical  tools  availa-^le  for  developing  con- 
stitutive relationships.  The  lecture  notes  were  purjK'sely  kept  to  an 
elementary  level,  and  were  prepared  vith  the  fov-nu ' ation  of  constitutive 
relations  for  earth  materials  in  mind. 

Objective 

2.  The  objective  of  this  report  is  to  dociraent  the  lecture  notes 
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in  a format  that  can  he  used  for  engineering  training  throughout  the 
Corps  of  Engineers,  U.  S.  Army,  or  as  materieJ.s  for  self-study  and 
reference  pruposes. 

Scope 

3.  Some  of  the  basic  mathematical  tools  necessary  for  the  develop- 
ment of  constitutive  relationships  are  presented  in  Part  II.  Included 
in  Part  II  are:  a brief  discussion  of  indicia!  notation,  matrix 
algebra,  development  of  basic  equations  related  to  eigenvalue  problem, 
the  Cayley-Hamilton  theorem,  and  Cartesian  tensors  (with  emphasis  on 
second-order  tensors).  A number  of  numerical  examples  are  included  in 
this  part  of  the  report  in  order  to  help  the  reader  to  better  understand 
the  subject  matter.  Part  III  includes  a summary  of  appropriate  equa- 
tions from  continuum  mecheuiics  required  for  this  elementary  presenta- 
tion of  the  subject  of  constitutive  relationships.  Constitutive 
equations  of  elastic  materials  are  developed  in  Part  IV.  The  so- 
cailed  incremental  constitutive  equations  are  discussed  in  Part  V. 
Constitutive  equations  of  simple  viscoelastic  materials  are  discussed 
in  Part  VI.  Constitutive  equations  of  plasticity  are  contained  in 


PART  II:  MATHEMATICAL  PRELIMINARIES 


P 

h 


I 
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It.  Some  of  the  basic  mathematiceQ  tools  necessary  for  treatmen^. 
und  understanding  of  the  physical  concepts  to  be  presented  in  the  en- 
suing p€irts  of  this  report  are  developed  in  this  part.  The  development 
is:,  kept  to  an  elementary  level  and  is  confined  to  orthogonal  Ceu'tesian 
coordinate  system.  In  order  to  establish  a common  basis  of  terminology 
and  notation,  both  indicial  and  matrix  notations  are  briefly  discussed. 
However,  indicial  notation  is  used  for  most  of  the  presentations 
throughout  this  report  in  order  to  keep  the  nxanber  of  equations  to  a 
minimum. 


Indicial  Notation 


5.  The  development  of  indicial  notation  is  based  on  a n\imber  of 
agreements  motivated  by  miniaturization  of  a large  system  of  equations 
or  variables.  For  example,  if  three  variables  are  denoted  by  , 

Xg  , and  Xg  , we  can  simply  denote  them  by  X^  , where  the  subscript 

i is  called  an  index  and  we  agree  that  it  takes  on  values  1,  2,  and  3 

(three-dimensional  geometry).  Similarly,  the  system  of  equations 

, Aj  = Xg  + Yg  , and  A^  ■=  X^  ♦ Y^  can  be  expressed  as 

A^  = Xj  + Y^  . An  index  which  is  not  repeated  in  any  single  term  is 

called  a free  index.  Thus,  the  index  i in  X^^  and  ~ 
a free  index.  Furthermore,  a free  index  must  appear  in  every  term  of 
an  expression.  S^/stems  which  depend  on  one  free  index,  such  as  X^ 
and  A^  , are  called  systems  of  first  order.  The  terms  X^  , Xg  , 
and  X^  are  called  the  components  or  elements  of  the  system.  A first- 
order  system,  therefore,  has  three  components.  Systems  which  depend 
on  two  free  indices,  such  as  A.j  , are  called  systems  of  second  order. 
Since  the  indices  take  on  values  1,  2,  and  3,  a second-order  system 
has  nine  components.  Similarly,  we  can  define  systons  of  third  order 
which  depend  on  three  free  indices  and  have  twenty-seven  components, 
e.g.,  A. . In  this  report,  however,  we  will  be  dealing  mainly  with 
first-  and  second-order  systems. 
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6.  If  an  index  appears  twice  in  a term  it  is  called  a dxmimy  index. 
For  example,  the  index  i in  A.^  is  a dummy  index.  By  agreement,  a 
dummy  index  implies  that  the  tern  is  to  be  suanmed  with  respjct  to  this 


index  over  the  range  of  the  index.  Thus 
X,Y,  ==  X,Y,  + X^Y^  + yi^Y^  , and  S.. . = C. 


^ii  ^ ^1  ^22  ^33  * 

'mm^ij  ' ^“^ll  ^22  ^33  ^^ij 

i and  J in  the  last  expression  are 

free  indices.  The  particular  letter  used  for  the  dummy  index  in  an 


’ 2 2 3 3 

It  is  noted  that  the  indices 


ij 


operation  is  immaterial;  thus.  A..  = A = A__  , 

li  pp  mm 


X.Y.  = X Y = X Y 
i 1 p p mm 


and  S.  = C E.  = C,  ,E,,  . This  characteristic  of  dummy  indices  is 
ij  mm  ij  kk  ij 

very  useful  for  manipulating  several  expressions  that  have  common 
indices.  For  example,  consider  the  following  expressions 


A = B C 
m r mr 


(1) 


B = D E 
r mr  m 


(2) 


In  the  first  expression  the  index  m is  free  and  the  index  r is  a 
dummy.  In  the  second  expression  the  index  r is  free  and  the  index  m 
is  a dummy.  The  index  r in  Equations  1 and  2 is  called  a connecting 
index.  If  we  substitute  the  second  expression  into  the  first  expres- 
sion and  use  the  same  letters  for  indices,  we  obtain 


A = D E C 
m nr  m mi* 


(3) 


Equation  3 is  meaningless  since  it  is  not  consistent  with  the  rules 
(agreement;')  of  indicial  notation;  the  index  m appears  three  times 
on  the  right-hand  side  of  this  expression.  To  obtain  the  correct  ex- 
pression we  must  first  overhaul  the  dummy  index  m in  Equation  2. 
Using  the  useful  characteristic  that  the  particular  letter  used  for  a 
duBsay  index  is  icEaaterial,  we  can  write 


B ^ 
r 


D E 
nr  p 


(M 
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Substituting  Equation  U into  Equation  1 

A = D E 
m pr  p mr 

Equation  5 is  notationally  correctj  there  is  no  question  as  to  which 
index  is  the  free  index.  Expanding  the  dunrny  indices  p and  r over 
their  range.  Equation  5 takes  the  following  form 


we  obtain 

C (5/ 


A =D,EC,+DoEC„+D_ECo 
m pi  p ml  p2  p m2  p3  p m3 


= ^lA^ml  " °12^1^m2  " ^13'^l^m3 


^21^2*^1111  * ^22^2^X02  * ^23^2^m3 


31  3 ml  32  3 m2  33  3 m3 


(6) 


Equation  6 (or  Equation  5)  has  three  components.  The  first  component, 
for  example,  becomes 


* “sASl  ♦ '>82^12  * Wu 


(7) 


which  is  quite  long  in  comparison  with  the  compacted  indicial  fora. 

7.  Another  agi*e«nent  in  establishing  indicial  notation  is  the  use 
of  commas  in  the  subscripts  to  represent  partial  derivatives.  Thus, 
we  agree  that 


F 

3X, 


(8a. 


3U^ 

3X-^iJ 


(8b) 
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Similarly, 


3U  3U 

p SL._® 

nk  ■ 3X^ 


U U , 
m,n  m,k 


(8c) 


In  Equavion  3o^  m is  a dummy  index  and  n and  k are  free  indices. 
Expanding  tr.e  dummy  inder  m , '.quatioi  8c  takes  the  following  form 


E u = U,  U,  . + 
ak  l,n  l,k 


U,  Uo  V + 
2,n  2,k 


^S.nS.k 


(9) 


Equation  9 (or  Equation  8c)  has  nine  components.  For  example,  the  E^^ 
component  becomes 


* 


^2,1^2, 3 


^3,1^3, 3 


(10) 


8.  In  indicial  notation  the  condition  of  symmetry  of  a second- 
order  system  is  denoted  by 


The  condition  of  skew-symmetry  is  denoted  by 


Equation  11  results  in  conditions 


®12  “ ®21 


®23  “ ®32 


«31  = "l3 


conditions 
of  syraaetry 


(11) 


(12) 


(13) 
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whereas  Equation  12  indicates  that 


^11  * ^22  ~ ^33  ~ ’ 


*^12  ~ "^21 


^23  ■ "^32 


^31  "^13 


conditions  of 
skew-symmetry 


Using  the  above  conditions,  an 


brie  (i.e.,  neither  symmetric  nor 


skew-symmetric)  second-order  system  can  be  expressed  as  the  sm 

of  a symmetrical  system  l/2^T^j  + and  a skew-symmetrical  system 

■ ’^ji)  ’ 

9.  In  using  indicial  notations,  one  often  deals  with  quantities 
that  have  no  free  index.  Such  quantities  are  referred  to  as  scalars 
or  zero-order  systems.  For  example,  the  following  quantities  are 
scalars 


A B 
on  nm 


D D D 
an  np  jan 


It  is  noted  that  ail  indices  in  Equation  l6  are  dummy  indices.  The  ex- 
panded form  of  the  last  expression  in  Eq\;ation  l6,  for  example,  becomes 


9 


D D D 
mn  np  pm 


=D,  D D,+D„D  D«+D^D 
In  np  pi  2n  np  p2  3n  np  p3 

= D^.D,  D , + D , + D-_D„  D - 

11  Ip  pi  12  2p  pi  13  3p  pi 

+ D*-  D-  D y.^  + D/s«D«  D « + D/s*D*  D ^ 

21  Ip  p2  22  2p  p2  23  3p  p2 

^3Ap^p3  ^32^2p^p3  ^33^3p^p3 

^ ^ll^lAl  * ®11°12^21  * ^11^13°31 


+ D /^D^-  /nD^/kD^-v-  + ^D-«D/k- 

12  21  11  12  22  21  12  23  31 


^13^31^11  * ^13^32^21  ■*■  ^13^33^31 


^21^iA2  ^21^12^22  * ^21^13^32 


+ D22^2i’^12  ^22^22^22  ^22^23^32 


* ^23^3^2  ^23^32^22  '*'  ^23^33^32 


+ D^,D,  ,D,  _ + D«-D, D- _D-- 
31  11  13  31  12  23  31  13  33 


+ D-„  + D-~D_-D,-,~  + 

32  21  13  32  22  23  32  23  33 


^33°31^13  ^33^32^23  ^33^33^33 


(IT) 


The  compactness  of  indicial  notation  is  once  again  demonstrated  hy  the 
above  expansion. 


Matrix  Algebra 

10,  another  convenient  method  for  representing  a large  number  of 
equations  or  quantities  is  through  matrix  notation.  A matrix  is  an 
array  of  numbers  or  components  of  a system.  For  example,  the  components 
of  a first-order  system  can  be  arranged  as 


10 


f'- 


{X}  = 7 X, 


(l8a) 


[X]  = [x^  X2  X^] 


(18b) 


Equation  l8a  represents  a S-by-l  (3  rows  and  1 column)  column  matrix 
whereas  Equation  l8b  represents  a l-hy-S  (l  row  and  3 columns)  row 
matrix.  Similarly,  the  components  of  a second-order  system  can 

be  arranged  as 


hi  ^12  ^13 


■ hi  •^22  ^^23 


hi  ^32  ^33 


Equation  19  represents  a 3-by-3  square  matrix.  We  are  mainly  interested 

in  3-by-3  matrices  in  this  report.  Some  useful  types  of  matrices  are: 

a.  Diagonal  matrix  in  which  all  elements  other 
than  those  on  the  diagonal  are  zero. 


BiiO  0 


[B]  = ° ^22  ° 


(20a) 


0 0 B, 


b.  Unit  matrix  in  which  all  off-diagonal  elements 
are  zero  and  every  diagoncLl  term  is  unity. 


..  ■ : J 


Ji:' 

P 

% 


i 


I 


I 


dl 


m 


ru  = 


10  0 
0 10 
0 0 1 


(20b) 


£.  Symmetrical  matrix  in  which  off-diagonal  terms 
are  symmetrical. 


[B]  = 


®11  ^12  ®13 


Bi2  Bg^ 


\3  ®23  ®33 


(20c) 


d.  Skew-symmetrical  matrix  in  which  every  diagonal 
term  is  zero  and  off-diagonal  terms  are  skew- 
symmetric  . 


[Cl  = 


r 0 c,  „ 

12 

13 

-Ci2  0 

^23 

■^13  ■■^23 

0 

(20d) 


In  indie ial  notation  the  counterparts  of  Equations  20c  and  20d  are  given 
by  Equations  11  and  12,  respectively.  Similarly,  in  indicial  form 


Equation  20a  can  be  expressed  as  B.  = 0 for  i ^ J . 

11.  The  transpose  [A]  of  a square  matrix  [A]  is  obtained  by 


completely  interchanging  every  row  with  its  corresponding  column; 


r 1* 

[A]  = 


^11  ^21  ^31 


AAA 
12  "22  "32 


[ ^3  ^23  ^33  J 


(21) 


12 


In  indicial  notation  A..  = A..  . In  view  of  Equations  11  and  12,  in 

ij  ji  n 

the  case  of  a symmetric  matrix  B . , = B . . , and  in  the  case  of  a skew- 
symmetric  matrix  C^j^j  = -C^j  . 

12.  Matrices  obey  certain  prescribed  rules  of  matrf.x  axgebra. 
Addition  or  subtraction  of  matrices  having  the  same  number  of  rows  and 
the  same  number  of  coltimns  is  accomplished  by  adding  or  subtracting 
corresponding  elements.  For  example,  consider  the  following  3-by-3 
matrices : 


\l 

®12 

^3 

II 

^21 

&22 

®-23 

_^31 

®’32 

^33_ 

'>11 

'*12 

^13  ■ 

[b]  = 

'*21 

'*22 

'’S3 

'^31 

”33 

(22a) 


(22b) 


Two  3'by-3  matrices  can  be  obtained  by  adding  or  subtracting  matrices 
[a]  and  [b]  ; thus. 


[c]  = [a]  + (b3  = 


^1  ' ’’U  ^2  * h2  ^3  * ^3 


+ 


®21  ^21  ^22  ^22  ^23  '*£3 

^3-1  * Si  ®32  ^32  ^33  * ^33 


(23a) 


[d]  = ta)  - [b]  = 


- b. 


®12  " ^2  ^3  “13 

®21  ''  ^21  ^2  ” ^22  ^23  " '*23 


^ " '*31  *^32  " '*32  ®33  ~ '*33 


(23b) 
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In  indicial  'notation  the  second-oi’der  systems  [c]  and  [d]  can  he  ex- 
pressed hy  ® a^^j  + hj^j  and  d^j  = a^j  - h^j  . A matrix  can  he 
multiplied  hy  a number  k hy  simply  mxiltiplying  every  element  in  the 
matrix  hy  k . Two  matrices  can  he  multiplied  together  if  they  are 
comformahle,  i.e.,  if  the  number  of  columns  of  the  first  matrix  is  eqxtal 
to  the  number  of  rows  of  the  second.  A p-hy-q  matrix  and  a q-hy-s 
matrix  axe  conformable  and  can  he  multiplied  together.  The  result  of 
multiplication  is  a p-hy-s  matrix.  For  example,  consider  the  multipli- 
cation of  matrices  [a]  and  [b]  given  in  Eqtiation  22: 

[a][hl  = [e]  (2U) 

The  matrix  [e]  is  a 3-hy-3  matrix  whose  components  are  obtained  from 
the  following  rule,  expressed  in  indicial  notation,  governing  matrix 
multiplication: 


®ij  “ ®ik\j 


(25) 


(e.g.,  = agj^h^^  + ^22^23  ^23^3  * ®22  ” ^21^2  ®22^22 

+ ^22^22)  • Eqxiation  25  it  should  he  noted  that  [a][h]  ^ lh][a]  . 

For  further  examples  of  matrix  multiplication  consider  the  following: 


[ar  = t®3 


[al^  = (nj 


(alCbKal  = (pi 


(artbl  » (ql 


[a]2(bl2  = tr] 


(althUal  [h]  = Is) 


(26) 


lU 


Using  the  rule  governing  matrix  multiplication  (Equation  25)  it  follows 
that  the  components  of  matrices  [m]  , [n]  , [p]  , iq]  , [rl  » and 

[s]  take  the  following  forms 


“ij 

""a 

hi 

hi 

hi 

hi 


^ik®kj 

^ik\f^fj 


^ik\f^fj 

hk\t^ti 


®ik\f^fg^gj 

^k\fVgj 


^ik\f^gf^gj 


(27) 


» 

Note  that  in  the  last  expression  in  Equation  27  the  definition  = a^ 
is  invoked.  Furthermore,  it  should  he  noted  that  in  Equations  25  and 
27  the  indices  i and  J are  the  only  free  indices. 

13.  The  sum  of  the  diagonal  terms  of  a square  matrix  is  called 
the  trace  of  the  matrix  and  is  denoted  by  tr  (e.g.,  trace  of 
[a]  = tr[aj  « a^^^  + a^g  + a^^)  ♦ In  indicial  form, 

trta]  = a^^^  (28) 

Similarly,  in  view  of  Equations  25  through  27, 


tr{(al[bl)  <= 
tr(ai^  = 
tr(a]^  = 
tr([al[hUa])  = 


^k\i 


^ik\i 

^ik^kf^fi 


^ik^kf'^fi 


(29) 
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All  indices  in  Equation  29  axe  dnsmy  indices,  indicating  that  the  trace 
of  a loatrix  is  a scalar.  Also,  tr([a]['b])  = trCthHa])  even  though 
[alih]  ^ [b][a]  . This  can  be  verified  by  expanding  the  indicial  form 

®ik\i  “ ^ik\i  * 

lU.  The  determinant  of  a square  matrix  is  denoted  by  det[a]  , or 
simply  la]  , and  is  expressed  as  (for  a 3-by-3  square  matrix) 


®11  ^2  ®13 
®21  ^22  ®23 
^31  ^2  ^3 


^1^  ^22^33  " ®23^32^  “ *12^^21^33  “ ^23^31^ 

+ a^3(a2i^32  “ ^22^31^ 


It  is  noted  that  the  determinant  of  a matrix  is  also  a scedar.  In  con- 
jxmction  with  the  determinant  of  a matrix  we  define  the  minor  and 
cofactor.  The  minor  of  an  element  a.,  of  the  matrix  [a]  is  obtained 

XV  ^v 

by  deleting  the  i “ row  and  j " column  and  forming  the  determinant  of 
the  remaining  terms.  For  example,  the  minor  of  a^^^  element  is 
given  €is 


minor  of  a^g  = 


®21  ®23 

V^3 


“ ^21^33  ’ ®23^ 


(31) 


The  cofactor  of  an  element  a^j  is  the  minor  of  that  element  with  a 
sign  attached  to  it  according  to  the  following  criterion 

cofactor  of  a^^j  = (-1)^**^  minor  of  a^^j  (32) 

Thus,  the  cofactor  of  a^g  el<mient  is  given  as 

1^2 

cofactor  of  a^g  “ ^*^21®^3  " ^23*^31^ 

« * ^3®31^ 

In  view  of  the  definition  of  minor  and  cofactor  the  determinant  of 
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the  matrix  [a]  can  he  expressed  as 

|a|  = a^j^Ccofactor  of  a^^)  *■  a^^gCcofactor  of  a^g) 

+ of  ^13^  (3U) 

It  shoxild  he  pointed  out  that  Equation  3^  is  not  unique  in  calculating 
the  determinant  of  the  matrix  [a]  . The  same  final  products  will  re- 
s\ilt  from  expansion  on  columns  or  other  rows,  e.g.. 


a I = a^gC cofactor  of  a^g)  + agg( cofactor  of  agg) 


+ a^gC cofactor  of  a^g) 


(35) 


15.  Finally,  we  define  the  inverse  [a]“^  of  a square  matrix  (a) 
such  that 


tal'^la]  . taltar^  = tU 


(36) 


The  inverse  matrix  is  given  hy 


(3T) 


where  the  matrix  [A]  , called  the  adLloint  of  [a]  , is  determined  hy 
replacing  the  elements  of  [a]  hy  their  corresponding  cofactors i thus. 


[A] 


agga^S  - a^gOgg 

^31^23  “ ®21^3 

.^1^2  " ^*^2 


^2®13  ■ *12^3 
*11*^3  " ^1*13 
^^2  * *11^32 


®12*^3  " *22®13 
®21®13  “ ^1*23 
*U®22  “ ^2^ 


(38) 


From  Equation  37  it  follows  that  the  inverse  exists  provided  |a|  #0  . 


Solutions  of  Linear  Algebraic  Equations 


16.  Consider  a set  of  linear  algebraic  equations 
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(39) 


^ ^2*2  * ^13^3  “ 

^1*1  *^22^2  ■*■  ^23*3  * ^2 
^31*1  ^32*2  ^33*3  " ^3 

In  indlclal  notation  Equation  39  c£m  te  expressed  as 


Vi  ■ "i 


(ItO) 


In  matrix  notation  Equation  39  takes  the  foUoving  form 


[a]{x)  = {k} 


(Ul) 


vhere  [a]  is  a square  matrix  of  coefficients,  {x}  is  a column  matrix 
of  unknowns,  and  {k}  is  a column  matrix  with  known  elements.  The 
objective  is  to  solve  for  the  elements  of  the  column  matrix  {x}  . Pre- 
multiplying both  sides  of  Eqxiation  Ul  by  [a]”^  results  in 

[a]'^[aKx}  = [ar^{k>  (U2) 


or,  in  view  of  Equation  36, 

M4{x)  = (x)  = la]"^{K}  (U3) 


From  Equation  U3  it  follows  that  once  the  inverse  of  the  coefficient 
matrix  is  determined  the  solution  for  (x)  is  obtai*>ed  by  perfoxming 
the  indicated  matrix  multiplication. 


Eigenvalue  Problem 


17.  In  a number  of  engineering  problems  the  following  systea 
of  algebraic  equations  is  often  encountered 


([el  -XriO){x}  • (0) 


(«*!*) 
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where  the  elements  of  the  column  matrix  {x}  are  the  unknowns  to  he 
determined,  X is  a scalar  parameter,  and  {0}  is  a null  column  matrix 
(all  elements  being  zero).  According  to  Equations  U3  and  37,  a non- 
trivial solution  of  Equation  UU  exists  only  if  the  determinant  of  the 
coefficient  matrix  vanishes,  i.e.. 


ail  - X 

^2 

^3 

^21 

®22  " ^ 

^23 

» 0 

(45) 

^31 

®33  ■ ^ 

The  expansion  of  the  above  determinant  yields  the  following  cid>ic  equa- 
tion in  X : 


X^  - I„X^  ■¥  ll\  - = 0 

a a a 


(U6) 


where 


la  * •'Sm 


(U7a) 


II  K sum  of  the  minors  of  the  diagonal  elements  of  [a) 
& 


III 


a 


®22  ®23 

^2  ^3 

IM 


^1  ®13 
*^1  *^3 


^L1  ®12 
^21  ^2 


(U7b) 

(l»7c) 


Equation  U6  is  called  the  characteristic  equation  of  the  matrix  (a)  . 

18.  The  three  roots  Xj^  , X^  , and  X^  of  the  characteristic 
equation  are  called  the  characteristic  values  or  eigenvalues  of  (a)  . 
For  evei*y  eigenvalue  X^  (assiiming  that  all  three  roots  are  distinct) 
Equation  46  is  satisfied  and  hence  Equation  44  has  nontrivial 
solutions: 
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{x(X^)} 


(W) 


Every  such  solution  of  {x}  is  called  a characteristic  vector  or 
eigenvector  of  [a]  . The  eigenvectors  {x(X^)}  corresponding  to 
eigenvalues  X.  can  be  grouped  together  to  foxn  a s(iuare  matrix,  re— 

^ A 

f erred  to  as  a modal  column  matrix « i.e.. 


[x]  = 


*11  *12  *13 
*21  *22  *23 
*31  *32  *33 


(U9) 


For  each  eigenvalue  and  the  corresponding  eigenvector.  Equation  UU  can 
be  written  as 

[a]{x(X.)}  = X^{x(X^)}  (50) 

Since  [a]  is  a 3-by-3  matrix.  Equation  50  can  be  expressed  in  the 
following  form  for  all  of  the  eigenvalues  X^  , X^  , and  X^  . 


'^1  ^2 

*11  *12  *13 

*11  *12  *13 

'Xj^  0*  0 ’ 

*21  *22  *23 

8 

*21  *22  *23 

0 Xg  0 

•^31  *^2  ^3 

*31  *3?  *33 

• m 

*31  *32  *33 

1 J 

0 0 x^ 

m ^ 

Equation  50  is  also  satisfied  if  each  eigenvector  is  multiplied  by  an 
arbitrary  constant  Cj^  , i.e.. 


lalcj{x(X^)}  » c^X.{x(X^))  (52) 

Therefore,  an  eigenvector  is  indeterminate  to  the  extent  that  it  can  be 
multiplied  by  an  arbitrary  constant.  Selecting  an  eigenvector 


(53) 


appropriate  to  eigenvalue  , the  corresponding  modal  column  matrix 
hecomes 


^11  ^12  ^13 

®1*11  *^2*12  °3*13 

[ul  = 

'‘Sl  ^22  ^23 

= 

CiXgi  CgXgg  C^Xg^ 

.^31  ^32  ‘‘33. 

C1X31  CgXgg  ^3X23^ 

■*11  *12 

^3‘ 

•Cj  0 0 1 

*21  *22 

*23 

0 

CVJ 

a 

0 

.*31  *32 

*33, 

L“  ° =3j 

(5U) 


It  is  obsei*ved  from  £q\iation  5^  that  a modal  column  matrix  is  indeter- 
minate to  ihe  extent  that  it  can  he  postmultiplied  hy  a diagonal  matrix 
of  arhitrazy  constants  c^  . Nov  utilizing  the  modal  column  matrix 
t Equation  51  can  he  exprrssed  as 


[al[yl  a IuKaJ  (55) 

where  fAj  is  a diagonal  matrix  with  elements  » Xg  , and  X^  . 
Premultiplying  Equation  55  hy  we  obtain 

luT^taKu]  = tAj  (56) 

Poatmultiplying  Equation  55  by  gives 

[a]  = (wirAjlur^  (57) 

From  Equation  56  it  is  observed  that  the  modal  column  matrix  [v] 
which  is  found  by  grouping  the  eigenvectors  of  [a]  diagonalizes  the 
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matrix  [a]  . Furthermore,  the  elements  of  the  diagonalized  matrix  are 
the  eigenvedues  of  [a]  . This  diagonalization  process  is  8ui  important 
part  of  the  eigenvalue  problem  and  its  significance  will  he  realized 
when  dealing  with  second-order  systems. 

19.  As  an  example  of  an  eigenvalue  problem,  consider  the  follow- 
ing system  of  equations 


(2  - X)x^  - Xg  + x^  = 0 

-8x^  + (3  - X)x2  + Tx^  = 0 

-8x^  ~ ^2  * - X)x2  = 0 


In  matrix  form  the  above  system  of  equations  is  expressed  as 

r 


2-11 
-837 
-8  -1  11. 


10  0 
0 10 
0 0 1 


or  (see  Equation  UU) 


([a]  - xrUXx}  = {0} 


The  chenracteristic  equation  of  [a]  is  given  as  (See  Equations 
and  U6) 


2 - X -1  1 

-8  3 - X 7 

-8  -1  11  - X 


-X^  ♦ 16X^  - 68X  + 80  = 0 


where  it  is  noted  that  I = I6  , II  = 68  , and  III  = 00  . Solution 

A A ft 

of  the  characteristic  equation  yields  the  following  eigenvalues  for  the 
matrix  [a]  : 


Xi  = 2 , Xg  » h , X3  a 10 
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For  each  value  of  X there  exists  three  homogeneous  equations.  For 
X = Xj^  = 2 we  have 


-Qxj^  + Xg  + Tx^  = 0 
-8x^  “ *2  ^*3  “ 

where  it  is  noticed  that  x^  = l,  XgSl,  ^3“^  nontriArial 

solution.  The  eigenvector  corresponding  to  X^  then  hecomes  (see 
Equation  53) 

Similarly*  for  X = Xg  = U , 

ivCXg)}  = Cg 

and  for  X = X^  = 10  , 

{p{X^)}  = c^ 


The  zaodal  coliuan  matrix  hecomes  (see  Eqvtation  5^) 


Cl  CgO 

h]  = Cl  -Cg  C3 


Using  Equation  56  it  can  "be  verified  that  the  modal  column  matrix  trans- 
forms the  matrix  [a]  into  a diagonal  matrix  with  elements  Xi  , , 


and  X^  } i.e., 


[ci  lire, 


2 0 0 


Cl  -Cg  c^  -8  3 T Cl  -Cg  Cg  = 0 0 


°1  ^^2  ^^3 


-8  -1  11  Cl  Cg  c^  00  10 


Cayley-Hamilton  Theorem 

20,  The  Cayley-Hamilton  theorem  plays  an  important  role  in  ex- 
pressing higher  powers  of  square  matrices.  It  simply  states  that  a 
square  matrix  satisfies  its  own  characteristic  equation.  The  result 
of  the  theorem  is  given  here  without  proof.  Let  [a]  he  a 3-'by-3  matrix 
and  its  characteristic  equation  be  given  as  (see  Equation  1*6 ) 


X^  - I X^  + II  X - III  =0 


If  [a]  satisfies  its  characteristic  equation  it  follows  that 


[al^  = lllJU  - II-[a]  + I Ja]^ 

& a <1 


Note  that  the  constant  III^  is  multiplied  by  a unit  matrix  flJ  • 


From  Equation  59  it  follows  that 


[a]^  = [al^La]  = ^ 6a  " 


i 

n 


Similarly, 


ta]5  = [al^a]  * (yil^  - 

* * (^a  - 

It  is  clear  from  the  examples  given  in  Equations  60  and  6l  that  using 
the  Cayley-Hamilton  theorem. (i.e..  Equation  59)  we  can  express  any 
power  of  [a]  greater  than  3 in  terms  of  [a]  and  [a]  . Accordingly, 

a polynomial  representation  of  [a]  , i.e.. 


[g]  = f([a])  = + k^[a]  + kgCa]^  + +...+  k^[a]°  (62) 


v 


if- 

ii: 


where  kQ  , k^  , . . . , k^  are  constants , can  he  expressed  as 

[g]  = nQ^I4  + rtj^la]  + 


where  the  coefficients  , and  are  now  polynomial  functions 

of  , and  III^  . 

21.  For  an  illustrative  example  of  the  Cayley-Hamilton  theoran, 
consider  the  following  matrix 


[c] 


12  0 
3 -1  -2 


LI  0 -3 


I 


The  characteristic  equation  of  [cl  is  given  as 

♦ 3X^  - 71  - 17  = 0 

where  it  is  ’^ted  that 

1^,-3 

n,  = -7 

nic  “ 
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i 

1 


We  substitute  [c]  for  X in  the  characteristic  equation  and  multiply 
the  constant  term  hy  a unit  matrix  tlj  , i.e.. 


1 

2 

0 

3 

1 

2 

0 

2 

1 

2 

- 

0 

3 

-1 

-2 

+3 

3 

-1 

-2 

-7 

« 

-1 

-2 

_ 1‘ 

0 

-3 

1 

0 

-3. 

ul 

0 

-3. 

1 

0 

0 

’3 

Ik 

12 

21 

0 

-12' 

17 

0 

1 

0 

= 

27 

-11 

-38 

+ 

-6 

21 

2lt 

0 

0 

1 

13 

-6 

-31 

-6 

6 

27 

" 7 lit 

o' 

17 

0 

m 

0 

0 

0 

. o' 

- 

21 

-7 

-lit 

- 

0 

17 

0 

= 

0 

0 

0 

7 

0 

-21 

0 

0 

17 

0 

0 

0 

resulting  in  a noil  matrix  [O]  . 

Cartesian  Tensors 

Cartesian  coordinate 

22.  Let  us  consider  the  orthogonal  Cartesian  coordinate  system 
Xj^  (Figure  l)  with  unit  vectors  i^^  , ig  » and  i^  along  the  , 
Xg  , and  x^  axes,  respectively. 


h 


Figure  1.  Orthogonal  Cartesian  coordinate 
system  x. 


From  elementary  vector  analysis  the  dot  products  of  these  unit  vectors 
are  given  as 


i *1  =1  •!  = •*  •!  =1 

1 1 2 2 “d  ^3 


^1  * ^2  ■ ^1  * ^3  " ^2  * ^3  “ ° 


or,  in  indicia!  notation. 


( 0 p jt 


This  product  is  denoted  by  6^^  and  is  known  as  the  Kronecker  delta; 


thus. 


i • i = 6 
p r pr 


iov^ 

U P == 


The  counterpart  of  6 in  matrix  notation  is  the  unit  matrix  flj 

PJ* 

(see  Equation  20b).  From  Equations  65  and  66  it  follows  that 


6 =3 

pp 


6 6 =3 
pr  rp 


Transformation  matrix 


23.  The  vector  V with  components  (x^^  , x^  , x^)  in  the  Xj^ 

coordinate  system  (Figure  2)  can  be  expressed  in  vector  forai  ar 

*2*2  * *3*3  ° *P*P  ■ 

If  we  fix  the  origin  and  rotate  the  axes  forming  a new  coordinate  system 
x^  , with  corresponding  unit  vectors  i^  (Figure  2),  then  the  vector 


/ f-rvv 


‘I'-T-'V'! 


V with  components  (xj^  , , xp  in  the  primed  (rotated)  system  can 

he  denoted  hy 

V = x^i;  (69) 

*2 


Figure  2.  Orthogonal  Cartesian  coordinate 
systems  x^  and 

In  view  of  Equations  66  and  69 


*p'p 


x'i* 

S 8 


(70) 


The  dot  product  of  Equation  70  with  i.  results  in 


1 


, Equation  71  reduces  to 


Since  i • i,  = 6 . 

p K pk 


x’i' 


s s 


i 


k 


By  the  definition  of  dot  product. 


i;  . ij^  = cos(x;  . Xj^) 


(72) 


(73) 


where  co8(x*  , x.  ) is  the  cosine  of  the  angle  between  the  x'  and 

S Iv  8 

Xj^  axes.  We  denote  cos(x^  , x^)  hy  a^j^ 

In  view  of  Equations  73  and  7^,  Equation  72  takes  the  form 


*k  = *^sk*^; 


(75) 


Similarly,  the  dot  product  of  Equation  70  with  i^  results  in  the 
following  relation 


v> 


s s 


(76) 


Equation  75  relates  the  components  of  the  primed  system  (rotated)  to 
the  components  of  the  unprimed  system.  Equation  76  relates  the  com- 
ponents of  the  unprimed  system  to  the  coe^nents  of  the  primed  system. 
The  matrix  a^j^  is  called  the  transformation  matrix  and  consists  of 
the  following  table  of  direction  cosines: 


2U.  Equations  75  and  76  can  now  be  utilized  to  establish  certain 
properties  of  the  transfonnation  matrix.  Differentiating  Equation  75 
with  respect  to  yields 

Since  \ i = » i*e. f x^  ^ = 1 , x^^g  “ ^ * ©"tc.*  Equation  77 

reduces  o 


*ki  ° Sk^B.l 


From  Equation  76,  x'  = a x , and  thus 

S Scl  91 


*8,i  “ ^sa*ta,i  “ ®sm^mi 


(78 


(79 


Substituting  Equation  79  into  Equation  78  results  in 


\i  ° *sk*^sm®ni 


(80 


In  view  of  the  definition  of  » 


Equation  80  reduces  to 


it  i'  iHiiiiiinri  iVnn  itiilM 


Similarly,  "by  differentiating  Equation  76  and  following  the  same  pro- 


cedure we  get 


®kp®ip  ^ki 


Equations  8l  and  82  describe  the  basic  properties  of  the  transformation 
matrix.  Expanding  Equation  8l  yields 


^21  ^31  " ^11  “ ^ 

2 2 2 

^2  * ^22  *32  ^ ^ 

2 2 2 

*13  * ^23  ^ ^33  = ^ 


*11  *12"^  *21  *22  * *31*32  " ^12  ~ ° 
*11  *13  * ^21  *23  * *31  *33  “ 

*12  *13  *22  *23  *32  *33  ' ° 


(83a) 


(83b) 


Similarly,  expanding  Equation  82  yields 


9 9 2 

*11  *12  *13  “ ^ 

4l  *22  *23  “ ^ 

2 2 2 

a^  ♦ a^  ♦ a33  = 1 


(8iia) 


*11  *21  *12  *22  * *13  *23 
*11  *31  * *12  *32  * *13  *33 
*21  *31  * *22  *32  * *23  *33 


(84b) 


Equations  83a  and  84a  indicate  that  the  sum  of  the  squares  of  the 
elments  of  any  column  or* cow  of  the  transfonaation  matrix  is  unity  and 
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called  normalization  conditions.  Equatiotxa  83b  and  8^b  Indicate  that 
the  sum  of  the  products  of  corresponding  elements  in  any  two  distinct 
columns  or  rows  is  zero  and  called  orthogonality  conditions . Through 
algebraic  manipulations  of  Equations  83  and  8U  it  can  also  be  shown  that 

la(^  = 1 (85) 


25*  For  a nvmerical  example  of  a transformation  matrix  coz^ider 
the  following  rotation  (Figure  3)  of  the  Xj^  coordinate  system: 


h 


/ 


Figure  3.  Rigid-body  rotation  of  x. 
coordinate  system 


The  transformation  matrix  a^  associated  with  this  rotation  can  be 
constructed  easily: 


0 1 


32 


I 

1 

i 

I 


I 


It  can  readily  be  verified  that  the  above  matrix  satisfies  Equations  83  | 

through  85.  j 

First-order  tensor  | 

26.  If  with  a coordinate  transforniation  = ^3*s  Equa-  1 

tion  76),  the  three  quantities  in  the  xinprimed  coordinate  system  | 

transform  to  three  quantities  A^  in  the  primed  reference  frame  by  i 

K.  = 


then  A is  a tensor  of  the  first  order.  We  already  know  that  A is 
a vector.  Therefore,  a vector  is  a tensor  of  first  order.  Within  the 
context  of  indicial  notation,  a first-order  tensor  is  a first-order 
system,  i.e,,  it  has  one  free  index.  Any  quantity  whose  value  does  not 

V 

chemge  with  coordinate  transformation  is  called  a tensor  of  order  zero 
or  a scalar  (see  Eq\iation  16),  A scalar  is,  therefore,  invariant  to 
rigid-body  rotation  of  the  coordinate  system.  Considering  the  scalar 
product  of  A^  with  itself  we  obtain 


^Tshs.  ” ®ks^s\p^p 


Since  = ^gp  (see  Equation  61),  Equation  87  becomes 


A^A^  = A_A_6 


s p sp 


In  view  of  the  definition  of  6 , Equation  88  reduces  to 

sp 


KK  - A-A 


s s 


(87) 


(88) 


(89) 


Equation  89  indicates  that  the  sum  of  the  square  of  the  elements  (com- 
ponents) of  a first-order  tensor  (vector)  is  invariant  to  rigid-body 
rotation  of  the  coordinate  axes.  This  quantity  is  the  only  invariant 
associated  with  a first-order  tensor.  The  magnitude  or  length  of  the 
vector  A is  given  as  a/  A A and  is,  therefore,  invariant  to  rigid- 

8 » 8 S 

body  rotation  of  the  coordinate  system. 


I 

'j 

'1 

i 

•1 
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27.  For  an  example  of  transformation  of  a first-order  tensor, 
consider  the  vector  Aj^  with  components 


in  the  coordinate  system.  The  magnitude  of  the  vector  is 

= ^(5)^  + (10)^  + (2)^  = VliT 

If  the  coordinate  system  undergoes  a rigid-hody  rotation  as  shown  in 
Figure  3,  the  components  of  the  vector  in  the  rotated  system  can  he 
calculated  from  Equation  86,  i.e., 

*i  ° * “l2*2  * “l3*3  • ^ <5)  + 4 ^ 

^2  “ *21^1  * ®22^2  * “23*3  “ ”2 

^ ■ ^1^1  * ^2^  ^ 

It  is  noted  that  the  magnitude  of  the  vector  is  not  affected  by  the  co- 
ordinate transformation,  i.e.. 

Second-order  tensor 

28.  Consider  two  first-ordsr  tensors  u^  and  v^  associated 
with  coordinate  system  . Since  u^  and  v^^  are  first-order 
tensors  we  may  write  (see  Bqxiation  7^) 


^ ' \i% 


(90a) 


■ '-'ri 


(90b) 


Combining  tiae  vectors  and  we  can  construct  the  second-order 
system  u^Vj  , which  we  me^  call  the  array  , i.e., 

Vi  V2  Vs' 


Vj  = '^ij  = 


Vi  V2  Vs 
Vi  V2  Vs 


In  view  of  Eqmtion  90  the  product  u^v^  can  he  written  as 


(91) 


■^ij  = Vj  = = \i^mj"X 


(92) 


Equation  92  provides  the  array  of  nine-numher  . Denoting  the  array 

u'v'  hy  T’  . Equation  92  becomes 
am''  nm  * ^ 


T.  . = T’ 

ij  ni  mj  nm 


(93) 


where  T'  is  referred  to  the  primed  coordinate  system.  Similarly, 
nm 

starting  from  ul  = a.  u and  v!  = a.  v^  , we  can  derive 
i in  n i im  m 


TI,  = a.  a,  T 
iJ  in  Jm  nm 


(91*) 


Any  quantity  that  transforms  according  to  Equation  9l*  is  called 

a second-order  tensor.  Within  the  context  of  indicia!  notation,  a 
second-order  tensor  is  a second-order  systen,  i.e.,  it  has  two  free 
indices.  Accordingly,  the  addition,  subtraction,  and  multiplication 
of  second-order  tensors  are  governed  by  the  rules  expressed  in  Equa- 
tions 23  throiigh  27.  In  matrix  notation  the  transformation  laws 
(Equations  93  and  9l»)  are  expressed  as 


[T]  = [al  (T’lta) 


[T'l  =>  la][THa] 


(95) 

(96) 
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vhere  [a]  is  the  transpose  of  [a]  . 

29.  The  second-order  tensor  is  an  extremely  imi)ortant  tensor  in 
mechanics  and  will  he  used  extensively  in  this  report.  In  partictilar, 
we  are  interested  in  second-order  symmetric  tensors  such  as  stress  and 
strain  tensors.  It  was  shown  in  Equation  89  that  there  is  one  invariant 
associated  with  a first-order  tensor  (vector).  In  the  case  of  a second- 
order  symmetric  tensor,  however,  there  are  three  independent  quantities 
that  remain  constant  with  respect  to  coordinate  transformation.  These 
independent  invariants  are 


- tr[Tl  . Tjj 

iij  = trdl"  = 
iiilj  = 


From  Equation  9^  it  follows  that 


tr[T' ) = T^^ 


a.  a.  T^ 
in  im  nm 


(97) 

(98) 

(99) 


(100) 


According  to  the  property  of  the  transformation  matrix  (Equation  81), 

a.  a.  - 6 , and  Equation  100  becomes 

in  im  nm 

tr[T')  = T (101) 

Tmi  Tim 

In  view  of  the  definition  of  6 , Equation  101  reduces  to 

nm 

tr[T'l  = « tr(Tl  (102) 

indicating  that  tr[T]  is  an  invariant.  Similarly,  from  Equation  9^* 

*ik  * ‘in’W'™  • ’'id  ' “kp^isV  • “* 
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= a.  .a,^  T 


in^lon^ nm^kp^is^ps 
^in®'is'^nin\ri\p'^ps 


(103) 


Again  using  the  property  of  the  transformation  matrix  (Equation  8l), 
^n^is  = ^ns  » V\p  = V ’ Equation  103  becomes 


tr[T']^  = 6 T 6 T 
ns  nm  mp  ps 


(lOU) 


In  view  of  the  definitions  of  <S  and  6 , Equation  lOlt  reduces  to 

ns  mp 


tr[T']^=T  T = tr[T]‘ 
sm  ms 


(105) 


indicating  that  tr(T]-  is  an  invariant.  Using  the  same  procedure  it 
can  be  shown  that 


tr[T']^  = tr(T]^ 


(106) 


indicating  that  tr(T]'^  is  also  an  invariant. 

30.  The  three  invariants  of  the  second-order  tensor  (l^  , 11^  , 
III^)  can  be  related  to  the  coefficients  in  the  characteristic  equation 
of  the  tensor  (Equation  U'r).  By  algebraic  manipulation  it  can  be  shown 
that 


II^  = - 2II.J,  (107a) 

IIj  = i (4  - iiij)  (107b) 


Nov,  using  the  Cayley  ^Hamilton  theorem  (£q\iation  ^9)  in  indicieJ.  form, 

X eC*  % 


- ^ViJ  ^ Wnj 


(108) 
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and  taking  the  trace  of  the  tensor  (putting  i = j)»  we  obtain 


w 

i 


T.,T,  T , = 3III„  - IlmT,.  + ImT.  T , 
ik  km  mi  T T is.  T in  ni 


(109) 


In  view  of  Equations  97 » 98 » 99*  and  107,  Equation  109  results  in 


= 3III„  - 


iiIt  = 


1 


1 


= :r  + •?•  I, 


T-^T  ^ T 


(llOa) 

(not) 


Equations  107  and  110  indicate  that  the  coefficients  in  the  character- 
istic  equation  of  the  tensor  stre  also  invariant. 

31.  For  an  i3.1ustrative  example  of  transformation  of  second-order 
tensors,  consider  the  following  tensor  associated  with  an  x,  co- 

iC 

ordinate  system: 


1;  1 2 

Ty  = 16  0 

.2  0 8_ 

From  Equations  97,  98,  and  99  we  have 

= tr[T]  * 18 
il^  = tr[T]^  = 126 
ni^  » tr[T]^  = 966 

Also,  from  Equation  U7, 


IIt  = 99 
III^  = 160 

where  it  is  noted  that  Equations  107  and  110  are  satisfied.  If  the 
coordinate  system  x^^  undergoes  a rigid-body  rotation,  such  as  the  one 

38 


shown  in  Fig\zre  3,  the  components  of  the  tensor  in  the  (rotated) 
syst^  can  he  calculated  from  the  transformation  law  of  second-order 
tensors  (Equation  9^*  The  transformation  matrix  associated  with  the 
coordinate  rotation  in  Figure  3 is  given  as 


•1*1 

2 2 

2 2 


0 

0 


0 0 1 


From  Equation  9^  it  follows  that 

“ ^in^jl\l  ^in^J2\2  ®in^J3'^n3 

^ ^il^jAl  ®i2®jl'^21  ^i3®J1^31 

^il^j2'^12  ^i2^j2'^22  ^i3^J2'^32 

^ ^i2^j3'^23  ^ ^3*^J3'^33 


Substituting  for  the  components  of  • and  a^j^  , we  obtain 


6 1 


t:  = 1 U -y'S 

-i/2  8 


Now,  utilizing  Equations  97»  98,  99*  and  U7  we  obtain 


1^,  = 18  = 
IT^,  = 126  = 

mij,  = 966  = 


39 


r 


I 

it 

I 

\ 

%■ 

V 

I 


1; 


m 


= 99  = 


III^,  = 160  s III^ 


indicating  the  inveiriant  nat\ire  of  these  quantities. 

32.  We  now  proceed  to  establish  some  useful  relationships  for 
second-order  tensors.  Consider  a second-order  symmetric  tensor  whose 
elements  in  the  coordinate  system  are  given  as  . Using 

the  transformation  law  of  second-order  tensors,  we  seek  a transformation 
matrix  that  will  transform  into  a diagonal  form  (i.e., 

= 0 for  i j)  associated  with  an  x^  coordinate  system.  The  axes 
x^  are  called  the  principal  axes  (or  principal  directions)  of  the 
tensor  and  the  elements  of  T’ ^ are  called  the  principal  values  of  the 
teusor.  A diagonalization  process  was  previously  demonstrated  for 
3-by-3  matrices  in  conjunction  with  the  eigenvalue  problem.  It  was 
shown  that  the  modal  coliimn  matrix,  which  is  found  by  grouping  the 
eigenvectors  of  a square  matrix,  diagonalizes  the  matrix  as  indicated 

by  Equation  56.  Furthermore,  it  was  shown  that  an  eigenvector  is  in- 

« 

determinate  to  the  extent  that  it  can  be  multiplied  by  em  arbitrary 
constant.  If  the  arbitral^  constant  is  chosen  to  be  the  inverse  of  the 
length  or  magnitude  of  the  eigenvector,  then  the  eigenvector  is  said  to 
be  normalized.  The  modal  column  matrix  of  normalized  eigenvectors  is 
called  a normalized  modal  column  matrix.  Denoting  the  normalized  modal 
column  matrix  by  [y]  , the  diagonalization  relation  (Equation  56)  fv>r 
the  matrix  [T]  can  be  written  as 


m = [yl"^[T][y] 


(111) 


In  view  of  the  properties  of  a treinsformation  matrix,  the  transformation 
law  of  a second-order  tensor  (Equation  96)  can  be  written  as 

-1 


[T«] 


k] 


[T][a] 


(112) 


Comparison  of  Equation  112  with  Equation  111  indicates  that  the  trans- 
pose of  the  normalized  modal  column  matrix  is  the  transfonaation  matrix 


ko 


ma 


which  transforms  T^j  into  a diagonal  form.  Furthermore,  the  elements 
of  the  diagoneiLized  matrix  are  the  eigenvalues  of  . In  the  case  of 

second-order  symmetric  tensors,  the  eigenvaltt>^8  (principal  values)  are 
always  real.  It  should  he  noted  that  the  normalization  of  eigenvectors 
is  necessary  in  order  to  conform  with  the  nurmalization  conditions  of 
the  transformation  matrix  (Equations  83a  and  81^a) . 

33.  For  a numericel  example  of  diagonal! zation  of  a second-order 
symmetric  tensor  consider  the  tensor  whose  elements  in  the 

coordinate  syston  are  given  as 


T.,  = 


-2  2 10 
2 -11  8 
10  8 -5 


The  characteristic  equation  of  is  given  as 


-2  - X 

2 

10 


2 

-11  - X 

8 


10 

8 

-5  - X 


= (X  - 9)(X  + 9)(X  + 18)  = 0 


The  eigenvalues  of  T^^j  are,  therefore. 


Xi  = 9 i Xg  = -9  i X^  = -18 


Next,  we  determine  the  normalized  eigenvectors  for  T^^  . For  X = X^ 
we  can  write  down  (see  Equation  UU) 


-llx^  + 2X2  + lOXj  = 0 
2Xj^  - 20X2  + Sxj  = 0 


lOx^  + 8xg  - lUx^  = 0 


Solving  the  above  system  of  equations  and  considering  the  normalization 


1*1 


condition  of  the  eigenvector  (i.e.,  + + the  normalized 

eigenvector  corresponding  to  heccmes 


The  normalized  modal  column  matrix  then  hecomes 


Ifa 


As  was  stated  previously,  the  transformation  matrix  which  transforms 
into  a diagonal  form  is  the  transpose  of  the  normalized  modal  column 
matrix.  This  can  he  verified  hy  using  [p]  as  the  transformation  matrix 
[a]  in-Etjuation  96,  i.e.. 


[T*]  = [a][T][a]*  = [uf[T][ul 


-2  2 10  X - T T 


2 -11  8 i I-  - 4 


10  8 -5  f i - 


Performing  the  above  matrix  operation  we  obtain 


900  0 0 

[T*]  = 0 -9  0 = 0 Xg  0 

0 0 -18  0 0 X, 


3I*.  Consider  three  second-order  symmetric  tensors  A.,  , , 

ij  mn 

and  C . Using  the  Cayley-Hamilton  theorem  it  was  shown  previously 
rs 

that  a polynomial  representation  relating  the  components  of  two  tensors 
takes  the  form  given  in  Equation  63*  In  indicial  notation  Equation  63 
is  expressed  in  the  following  form 


hs  “ ^ Vik\j 


(113) 


The  counterpart  of  Equation  113  expressing  the  components  of  one  tensor 
in  terms  of  the  components  of  two  other  tensors  was  derived  by  Rivlin 
and  Ericksen.  The  Rivlin-Ericksen  eq\iation  given  here  without  proof 
has  the  following  form 


hi  " '^ij'^mn  * ^rs^  ' ^O^iJ  * Wi  * ’’2®ik®kJ 
^3^iJ  * ^l*°ik^kj  * ^5^®ik^kJ  ^ik\j^ 
**■  ’’6^®ik\pSj  *^ikV®pJ^ 

^8^®ik®kpSt^tJ  * ^ik^®pt®tj^ 


ink) 


where  the  coefficients  Hq  « ...»  hg  are  polynomial  functions  of 

the  invariants  of  and  C and  the  following  Joint  invariants 

BUI  rs 


= B vCv 
1 ab  ba 


”2  “ \bS>c^ca 


n_  = B V B,,  C 
3 ab  DC  ca 


(115) 


"U  “ ®ab\c^ 


cd^da. 


It  is  noted  that  when  dependence  on  C disappears*  Equation  llU  re- 

IPS 

duces  to  Equation  113*  Equations  113  and  llU  are  the  bases  for  most  of 
the  presentations  in  the  ensuing  parts  of  this  report. 
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PART  III:  SUMMARY  OP  BASIC  COHCEPTS  FROM  CONTINUUM  MECHANICS 


Stress  Tensor 


3^.  In  Cartesian  coordinate  system  , ve  define  the  stress 


tensor  at  a iwint  as 


a,,  = limit  -r^ 
A^  -*■  0 ^i 


(116) 


where  F is  force  in  the  coordinate  direction  J and  A.  is  the  area 

J X ^ 

normal  to  i axis  on  which  the  force  Fj  acts.  Figure  U depicts  the 


r |/' 

> 

“Si  / 


Figure  U.  Stress  components 

positive  directions  of  the  components  of  the  stress  tensor.  In  the 
absence  of  distributed  body  or  surface  couples  the  stress  tensor  is 
sysDetrical*  i.e.*  o^j  * °J1  * Accordingly*  the  state  of  stress  at  a 
point  can  be  described  by  six  independent  stress  components. 

Invariants  of  stress  tensor 

36.  Stress  tensor  is  a second-order  syMetric  tensor  and  it  obeys 


1»5 


the  transfonnation  lav  given  in  Equation  9^,  l.e., 


®’in^jm^  nm 


(117) 


where  is  referred  to  (rotated)  coordinate  system.  As  was 

shown  in  Part  II,  a second-order  tensor  has  three  independent  invariants 
(Equations  97,  98*  and  99).  In  the  case  of  stress  iansor  we  define 
these  invariants  as 


= a 


nn 


o “ 2 °ik®ki 
3 ^ik^^km'^mi 


(118) 

(119) 

(120) 


Stress  deviation  tensor 

37.  Stress  tensor  can  he  expressed  as  the  sum  of  two  second- 
order  symmetric  tensors  in  the  following  manner 


(121) 


where  the  tensor 


hi  ' ®ij  ’ 3 V^iJ 


(122) 


is  referred  to  as  the  stress  deviation  tensor  and  o 6.,/3  is 

nn  ij' 

called  the  spherical  stress  tensor.  An  important  property  of  the  stress 
deviation  tensor  is  that  its  trace  is  equal  to  zero,  i.e.. 


^ii  “ 


'ii 


- 0 = 
nn 


(123) 


The  stress  deviation  tensor,  therefore,  has  only  two  independent  invari- 
ants. We  denote  these  invariants  as 


1*6 


I* 


k 

Si 


“ 2 ° 2 

(I2I*) 

5 “S  ' 1 

(125) 

The  invarieuits  of  stress  and  stress  deviation  tensors  can  be  related  by 
using  Equation  122.  In  view  of  Equations  122  and  ll8, 

^2  = I 

. “ 2 (^ik  " 3 *^l^ik)  (‘^ki  " 3 *^1^1) 

' ■|(°ik\i  " 3 ‘^I'^ik^i  ■*■  "9  ^ik^kij 


(126) 


Since  = 3 , , and  = 21^  , Equation  126 

becomes 


Similarly*  it  can  be  shown  that 


73  = ^3-!V2*^ 


(127) 


— j3 
27  *^1 


(128) 


Principal  stresses 

38.  The  three  principal  values  of  stress  tensor  are  referred  to 
as  principal  stresses  and  are  denoted  by  (using  the  principal  directions 
as  reference  axes) 


[0]  = 


0^  0 0 

0 Og  0 
0 0 0- 


(129) 


It  should  be  pointed  out  that  the  ordering  of  the  principal  stresses  in 
Equation  129  does  not  imp"':'  that  the  numerical  value  of  is  greater 
than  Og  :ir  . As  discussed  in  Part  II,  the  three  principal  values 
are  the  roots  of  the  characteristic  equation  of  the  tensor 

1»7 


- I X^  + II  X - III„  = 0 
0 0 0 


(130) 


where 


^0  ~ *^1  ” *^iin 


‘^22°'23 

''ll°13 

°11®12 

s 

+ 

+ 

‘^32^33 

°31°33 

®21°22 

= I”! 


(131a) 


(I31h) 


(131c) 


The  two  coefficients  II^  and  III^  are  usually  denoted  by  and 
t respectively,  and  can  be  related  to  , Jg  , and  by  using 
Equations  lOT,  110,  ll8,  119,  and  120. 


= = f (4  - ^2) 


= Illg  « J3  - ^2^1  i *^1 


(132) 

(133) 


The  invariants  of  stress  deviation  tensor  can  also  be  expressed  in  terms 
of  , 7g  , and  . In  view  of  Equations  127,  128,  132,  and  133, 


we  obtain 


j' = J3  - i * If 


(13M 

(135) 


Piincipal  stress  space 
and  octahedral  stresses 

39*  Since  the  three  principal  stresses  are  orthogonal,  they  form 
a three'>dimensional  space  called  the  principal  stress  space  (Figure  5)> 
Of  particular  interest  in  the  principal  stress  sx^e  are  the  octahedral 
planes.  The  direction  cosines  of  a normal  to  an  octahedral  plane  are 
(Figure  5) 


2i8 


- •-  - -1  -• 


/ N=NORMAI  . TO  OCTAHEDRAL 
! PLANE 


‘^oct"'3 


fl-5l*44' 


Figure  5.  Principal  stress  space 


cos  (N  , a^)  = cos  (N  , Cg)  = cos  (N  , a^)  = cos  (5^°  UL') 


(136) 


The  normal  and  shear  stresses  on  octahedreQ.  planes  are  denoted  as 

and  , respectively.  The  magnitude  of  and  can 

be  determined  from  the  transformation  law  of  stress  tensor^ 


= f---l 


0 0 


0 Og  0 


0 0 oJf~ 


3 ^°1  '"2 


(137) 


■ ' 

kv;  '.'M 


and 


= (o^  - 0^)^  + - <^3)^  + (^2  ■ ^^38) 

Using  Equations  121  and  II8  it  can  be  shovn  that  for  a general  state  of 
stress 

-i 

^oct  3 

vt 


(139) 

(ll*0) 


Equations  139  and  lUo  indicate  that  the  octahedral  stresses  are  also 
invariant.  The  octahedral  space  versus  0^^^)  is  commonly  used 

for  plotting  stress  paths  for  various  laboratory  tests.  In  this  report 


we  will  use 


Vi/S 


T . versus  o . 
oct  oct 


space  (i.e.,  "ersus  J^/3) 


for  defining  stress  paths. 


Examples  of  sim- 
ple  states  of  stress 

1*0.  The  following  states  of  stress  are  often  utilized  in  the 
laboratory  in  order  to  determine  the  stress-strain  properties  of  a 
material : 


Spherical  or  hydrostatic  state  of  stress. 


0 


3 
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b.  Uniaxial  state  of  stress. 


0 0 1 


0 0 0 


0 0 0 


c_.  Cylindrical  state  of  stress. 


0 0 

"ij  = "3  ° 

0 0 0. 


°2  = °3 


d.  Triaxial  state  of  stress. 


‘^l  I 


0 0 
0 0 
0 0 0., 


°3  “TA- 

/Ji 


• ‘^2  °1 


e.  Pure  shear. 


0 0,2  0 

^ij=  '^Pl®  0 

0 0 0 


Note  that  in  examples  a through  ^ all  stresses  are  principal  stresses . 
Stress  paths  associated  vith  the  above  states  of  stress  can  be  readily 
defined  in  the  versus  Jj/3  space.  The  stress  path  associated 
vith  spherical  or  hydrostatic  state  of  stress  is  shovn  in  Figure  6a,  It 
is  noted  that  for  spherical  state  of  stress  is  zero.  The  stress 
path  associated  with  uniaxial  state  of  stress  is  shovn  in  Figxire  6b. 

For  uniaxial  state  of  stress  VTj  = o,/>^  and  J,/3  = o,/3  result- 
ing in  the  expression  V^=  yj  for  the  stress  path. 
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a.  SPHERICAL  OR  HYDROSTATIC  b.  UNIAXIAL  STATE  OF  STRESS 

STATE  OF  STRESS 


3 


c.  CYLINDRICAL  STATE  OF  STRESS  d.  CYLINDRICAL  STATE  OF  STRESS 

(CONSTANT  Jj/3)  (CONSTANT o^) 


3 

«.  PURE  SHEAR 


Figure  6.  Stress  paths  associated  with  simple  states  of  stress 
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Figures  6c  and  6d  depict  special  stress  paths  associated  with  cylindrical 
state  of  stress.  In  Figure  6c  the  material  is  first  loaded  hydro- 
statically and  then  sheared  while  J^/3  is  kept  constant.  In  Figure  6d 
the  material  is  first  loaded  hydrostatically  and  then  sheared  hy  in- 
creasing while  keeping  constant.  Since  for  cylindrical  state 
of  stress  ^3^/3  = (c^  + 20^)73  it  follows 

that  the  expression  for  the  stress  path  of  Figure  6d  becomes 

=s  (J3^/3  “ ^2^  • stress  path  associated  with  pxire  shear 

test  is  shown  in  Figxire  6e.  In  the  case  of  p\ire  shear  J.^/3  = 0 . In 
the  actual  laboratory  coordinate  systen,  the  stress  components  and 
Og  = Og  associated  with  cylindrical  state  of  stress  are  usually  de- 
noted by  0 (axial  stress)  and  o =0,.  (radial  stress),  respectively, 
a r 0 

In  the  case  of  triaxial  state  of  stress,  the  stress  components  , 

Og  , and  are  denoted  by  , and  , respectively.  For 

p\xre  sheeir  the  only  nonzero  stress  component  is  generally  denoted 

by  T . 


Strain  Tensor 


Ul.  Let  us  consider  a cylindriced  specimen  of  length  and 
extend  it  to  length  I . The  ratio  is  defined  as  the  stretch 

X . 


X = 1/Iq  (lUl) 

The  question  is,  what  is  the  axial  strain  in  the  specimen?  There  are 
several  measures  of  strain  that  cem  be  used  to  determine  the  axial 

3 

strain  e in  the  specimen.  These  measures,  named  after  Cauchy,  Green, 
Hencky,  Almansi,  and  Swainger,  respectively,  are: 

= X - 1 (lU2a) 

e°  = i(x2-l)  (lUa) 
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e®  = ilnX 


(ll+2c) 


e'^  = i M - i-. 


(llt2d) 


e®  = 1 - i 


(l42e) 


In  order  to  demonstrate  the  difference  hetween  the  various  measures  of 
strain  given  in  Equation  l42,  let  S.  = (i.e.,  let  the  length  of 

specimen  he  doubled).  The  stretch  X = 2 in  such  case  and  from  Equa- 
tion 1^2  it  follows  that 


e = 100^ 


= 695^ 

= 37. 5^ 
e®  = 50?  , 


(1U3) 


As  observed  from  Equation  11*3 » for  a stretch  of  X = 2 , the  difference 
hetween  the  various  measures  of  strain  is  quite  appreciable.  Now  let 
£ = 1.25£q  » which  gives  a stretch  of  = 1.25  . -i-n  view  of  Equa- 
tion ll*2,  for  X = 1.25  the  various  measures  of  strain  become 


= 25? 


e®  = 22? 
= 18? 
G®  e 20? 


ilhh) 


It  is  observed  in  this  case  that  the  difference  between  various  measures 
of  strain  is  not  as  appreciable  as  was  the  case  for  X = 2 . If  the 


stretch  X is  further  reduced,  say  X = 1.1 
in 


Eqxxation  lk2  will  resvilt 


(11+5) 


Therefore,  for  small  deformations  (infinitesimal  strain  theory)  the 
various  measures  of  strain  will  yield  approximately  the  same  results. 

Our  interest  here  is  also  within  the  framework  of  infinitesimal  strain 
theory  and  we  adopt  the  Cauchy  measure  of  strain  for  further  analysis. 

U2.  In  order  to  determine  strain-displacement  relations  and  de- 
fine the  infinitesimal  strain  tensor,  we  consider  a particle  P with 
position  vector  x^^^  in  the  x^  coordinate  system  as  shown  in  Figure  T. 


Figure  7.  Particle  displacement  in 
coordinate  system 
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We  assume  that  the  particle  undergoes  displacement  and  assumes  a 
new  position  vector  as  depicted  in  Figure  7.  From  Figure  7 we  ceui 
write 


+ Ui  = x^ 


ilk6) 


or 


= Xi  - x^ 


(11+7) 


Since  x^  is  a function  of  x^  , i.e.,  x^  = differen- 

tiate Equation  ll+6  with  respect  to  x^  ; thus. 


^Ij  ’^i.J  " *i.j 


(148) 


The  terms  x.  . and  u.  , are  called  the  coordinate  gradient  and  dis- 

1 ^ J 1 ^ J 

placement  gradient  matrices,  respectively.  The  displacement  gradient 
laatrix  can  he  expressed  as  the  sum  of  a symmetrical  system  and  a skew- 
symmetrical  system  (see  Equation  15) 


(Ik?) 


The  first  term  in  Equation  1U9  is  symmetrical  and  is  called  the  infini- 
tesimal strain  tensor  ; thus. 


(150) 


The  second  term  in  Equation  l49  is  skew-symmetric  and  is  called  the 
rotation  tensor  ; thus. 


“ij  = 2 S.J  - "J.i* 


(1.1) 


Equation  150  relates  the  components  of  infinitesimal  strain  tensor  with 
components  of  displacement  vector. 
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43.  To  demonstrate  the  application  of  Equation  150,  consider  a 
rod  of  length  extended  to  length  I as  shovn  in  Figure  8. 


Figure  8.  Rod  in  uniaxial  extension 


The  boundary  conditions  associated  with  displacement  u^  in  the  x^ 
direction  are 


» 0 at  X,  * 0 


= £ - iQ  at  x^  = 


For  a homogeneous  state  of  strain  to  exist  in  the  rod,  the  displacement 
must  be  a linear  function  of  x^  . Thus, 

u^  = Cx^  (153 

where  C is  a constant.  In  view  of  Equation  152,  Equation  153  becomes 

■ -T-^  Xi  <151. 


Substituting  Equation  154  into  Equation  150  we  obtain 


£ - £, 


"u  • ? <“1,1  ♦ "i.i‘  ° ^.l  “ TT 


which  is  the  Cauchy  measure  of  strain  (see  Equation  lU2a)., 

Invariants  of  strain  tensor 

Strain  tensor  is  a second-order  tensor  and  obeys  the  trans- 
formation law  given  in  Equation  9^,  i.e.. 


e! . = a,  a,  e 
ij  in  ran 


(156) 


where  is  referred  to  the  coordinate  system.  There  are» 

therefore,  three  independent  invariants  associattid  with  the  strain 
tensor.  As  for  the  invariants  of  stress  tensor,  we  define  the  invari- 
ants of  strain  tensor  as 


^nn 

ie  e 

2 ik  hi 

3 ^ik'^km^^mi 


(157) 

(158) 
(159) 


Strain  deviation  tensor 

1*5.  Strain  tensor  can  he  expressed  as  the  sum  of  two  symmetric 
tensors  in  the  following  manner 


+ 


1 

3 


e 6. , 
nn  ij 


(160) 


where  the  tensor 


nn^ij 


(161) 


is  referred  to  as  the  strain  deviation  tensor  and  e ^6. ,/3  is  called 

nn  ij 

the  spherical  strain  tensor.  As  for  invariants  of  stress  deviation 
tensor,  we  deflre  the  invariants  of  strain  deviation  tensor  as 
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(163) 


liiHiliMiw 


V.  =i 


The  invariants  of  strain  deviation  tensor  can  also  he  expressed  in  terms 
of  the  invariants  of  strain  tensor  as  follows: 


Tl  - I _ i. 

^2  ^2  6 ^1 


^3 


— 2 — 2 
I - — II  +~I^ 
^3  3 -^1  2 27  1 


(161:) 

(165) 


Principal  strains 

U6.  The  three  principal  values  of  strain  tensor  fire  referred  to 
as  principal  strains  and  are  denoted  by  (using  the  principal  directions 
as  reference  axes) 


[ej  = 


0 


(166) 


The  princip-al  strains  are  the  roots  of  the  characteristic  equation  of 
strain  tensor 

+ II  A - III  = 0 (167) 

& O C 


where 


"1  = 


nn 


(l68a) 


^22  ^23 

^11  ®13 

*^12 

♦ 

+ 

®32  ®33 

^31  ^"33 

*^21  *^22 

III  = |e| 


(168b) 


(l68c) 


it'. 


r- 


! 


I 

I 


I 


The  two  coefficients  II  and  III  are  denoted  hy  I-  and  I-  , 

E E i ^ ■ 

respectively,  and  are  related  to  the  invariants  I^  , . Ig  , and 
as  follows: 


(169) 


I 


3 


^2^1 


(170) 


The  invariants  of  strain  deviation  tensor  can  also  he  e:q>ressed  in  terms 
of  I^  , Ig  , and  I^  : 


I*  = — I^  - I 
2 3 1 2 


(171) 


— II 
3 12 


(172) 


Examples  of  simple 
states  of  deformation 

U7.  The  following  states  of  deformation  are  often  utilized  in 
the  laboratory  in  order  to  determine  the  stress-strain  properties  of 
the  material; 

a.  Uniform  dilatation. 


0 


^ij  ' 3 ^ij 


b.  Uniaxial  state  of  strain. 
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0 O' 
0 0 0 


0 0 0 


£.  Cylindrical  state  of  strain  (sg  = e^). 


0 0 


0 Eg  0 

0 0 e. 


d.  Trieixisil  state  of  strain. 


0 0 
0 Eg  0 

0 0 e- 


£.  Simple  shearing  deformation  (no  volume  change). 


0 0 
ty  = 0 0 

0 0 0 


Note  that  in  examples  a through  ^ all  strains  are  principal  strains, 

Strain-Rate  Tensor 

1*8.  The  time  derivative  of  infinitesimal  strain  tensor  is  re- 
ferred to  as  rate  of  infinitesimal  strain  tensor,  or  simply  strain-rate 
tensor,  l thus. 


'ij  dt  '^ij^ 


(173) 


where  d/dt  indicates  differentiation  with  respect  to  time.  In  view  of 
Equation  150,  the  strain-rate  tensor  takes  the  form 


ink) 


where  v^  = components  of  velocity  vector. 

Invariants  of 
strain-rate  tensor 

U9.  Strain-rate  tensor  is  a second-order  symmetric  tensor  and, 
like  the  stress  and  strain  tensors,  it  obeys  the  transformation  law 
given  in  Equation  9^.  Similarly,  we  define  the  invariants  of  strain- 
rate  tensor  as 


i,  * I.  = e 

1 e nn 

(175) 

= — iT  = — i i 

2 t 2 ik  ki 

(176) 

^ m.  = ^ t e . 

3 e 3 ik  km  mi 

(177) 

Strain-rate  deviation  tensor 

50.  Strain-rate  tensor  can  he  expressed  as  the  sum  of  two  sym- 
metric tensors  in  the  following  manner 


^nn^ij 


(178) 


where  the  tenser 


3 ^nn^ij 


(179) 


is  called  the  strain-rate  deviation  tensor  and  i 6j,/3  is  called 

nn  ij 

the  spherical  strain-rate  tensor.  We  define  the  invariants  of  strain- 
rate  deviation  tensor  as 


(180) 


2 2 " E 2 


^3  ' I ' 3 ^ 


(181) 


The  invariants  of  strain-rate  (ieviation  tensor  can  also  be  expressed  in 
terms  of  the  invariants  of  strain-rate  tensor: 


— — 1 *2 

I*  = I - ^ I 

2 2 Z 1 


(182) 


— — 2 • — 2 *3 

T*=T  - — I I + — I 

^3  3 3 ^1  2 27  1 


(183) 


Principal  rates  of  strain 

51.  The  three  principal  values  of  strain-rate  tensor  are  denoted 
by 


H]  = 


0 0 


0 

0 0 


^2  ° 


(181*) 


and  sore  called  the  principal  rates  of  strain.  The  principal  rates  of 
strain  are  the  roots  of  the  characteristic  equation  of  strain-rate 
tensor 


+ IlgX  - = 0 


(185) 


where 


e 1 nn 


(186a) 
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(1861)) 


^22  ®23 

‘11  '13 

+ 

+ 

t & 

^32  '"as 

£ i 

31  33 

t i 

21  22 

III^  = |£|  (l86c) 


The  two  coefficients  and  III^  are  denoted  by  ig  and  » re- 

spectively, and  are  related  to  the  invariants  , Ig  , and  as 
follows : 


(187) 


^2^1 


(188) 


The  invariants  of  strain-rate  deviation  tensor  can  also  be  expressed  in 

t * • 

terms  of  , Ig  , and  : 


(189) 


— « 1 • • ? 

I'iT  - — II  + — I 

3 3 3 1 2 27  1 


(190) 


Equations  of  Continuity  and  Motion 


52. 


The  motion  of  any  continuum  is  governed  by  the  following  laws: 

a.  Conservation  of  mass. 

b.  Conservation  of  energy. 

£.  Balance  of  linear  momentum. 

Balance  of  angular  momentum. 
e.  Principle  of  inadmissibility  of  decreasing  entropy. 


6l» 


^-  raiMilwiiw«iBti.iiw^ 


These  lavs  constitute  the  hasic  axioms  of  continuum  mechanics.  In  the 
absence  of  distributed  couples,  the  balance  of  angular  momentum  leads  to 
the  symmetry  of  stress  tensor,  mechanical  energy  is  the 

only  form  of  energy  to  be  considered  in  a problem  (as  is  the  case  in 
this  report),  the  above  principles  lead  to  the  continuity  equation 


It  (pVi)^i  = 0 


(191) 


and  the  equations  of  motion 


“iJ.J  * ■ f*! 


(192) 


where  p = mass  density  , v^  = components  of  velocity  vector  , 

~ components  of  body  force  , and  a^  = components  of  acceleration 
vector  . Equations  191  and  192  are  applicable  to  all  materials. 

Constitutive  Equations 

53.  Equations  191  and  192  constitute  four  equations  that  involve 
ten  unknown  functions  of  time  and  space:  the  mass  density  p , the 
three  velocity  components  v^  , and  the  six  indept  t stress  com- 
ponents . The  body  force  components  f^  art  n quantities  and 
the  acceleration  components  a^  are  expressible  in  terms  of  the  veloc- 
ity components  v^^  . Obviously,  Equations  191  and  192  are  inadequate 
to  determine  the  motion  or  defonaation  of  a medium  subjected  to  ex- 
ternal disturbances,  such  as  surface  forces.  Therefore,  six  additional 
equations  relating  the  ten  unknoi/n  variables  p , v^  , and  are 

required.  Such  relationships  are  referred  to  as  constitutive  equations, 
which  relate  the  stress  tensor  to  deformation  or  motion  of  the 

mediimi.  As  was  pointed  out  previously.  Equations  191  and  192  are 
applicable  to  all  materials,  whereas  constitutive  equations  represent 
the  intrinsic  response  of  a particular  material.  Furthermore,  a con- 
stitutive equation  provides  a matheiaatical  description  or  definition  of 
an  ideal  material  rather  than  a statement  of  a universal  lav.  The 
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general  form  of  a constitutive  equation  may  be  expressed  by  the  func- 
tional form  (considering  only  mechanical  effects) 


where  = time  derivative  of  stress  tensor  . Equations  191,  192,  and 
19!*  (or  Equation  193)  , therefore,  constitute  ten  equations  in  ten  un- 
knowns and  will  lead,  in  conjxmction  with  kinematic  relations  given 
by  Equations  150  and  1?^,  to  a complete  description  of  the  boundary- 
value  problem.  In  addition  to  the  above-mentioned  equations,  boundary 
conditions  in  terms  of  boundary  displacement  and/or  surface  forces  must 
also  be  specified  to  completely  define  a particialar  problem  of  interest. 

5I+.  In  order  for  constitutive  equations  to  describe  physical  ma- 
terials adequately,  the  functional  forms  f^^  or  g^j  must  remain 
invariant  with  respect  to  rigid  motion  of  spatial  coordinate  . This 

requirement  stems  from  the  fact  that  the  response  of  a material  is 
independent  of  the  motion  of  the  observer.  Furthermore,  the  functionals 
f..  or  g . . must  be  consistent  with  the  general  princdples  of  con- 

1«J  1 J 

servation  or  balance  of  mass,  momentum,  and  energy. 

55.  We  adopt  Equation  19*+,  relating  four  second-order  symmetric 

tensors  s , £ , 0 , , and  d , , as  a basis  for  development  of 

rs  ’ mn  ’ ab  ’ cd  ’ ^ 

various  constitutive  equations  in  the  following  parts  of  this  report. 


PART  IV;  CONSTITUTIVE  EQUATIONS  OP  ELASTIC  MATERIALS 


56.  For  an  elastic  material,  the  state  of  stress  is  a function  of 
the  current  state  of  strain  only.  PurtheiTnore , an  elastic  material  re- 
turns to  its  initial  state  after  a load-unload  cycle  of  deformation  (no 
permanent  strain).  The  stress  tensor  can,  therefore,  he  expressed  in 
terms  of  strain  tensor 


c. . = F ,(£  ) 

ij  Ij  mn 


(195) 


where  = elastic  response  function.  Two  different  procedures  have 
been  utilized  in  order  to  determine  the  response  function  F^^^  for  iso- 
tropic materials.  The  first  procedure,  referred  to  as  Cauchy's  method, 
is  based  on  the  Cayley-Hamilton  theorem  (Equation  59 )•  The  second  pro- 
cedure, referred  to  as  Green's  method,  is  based  on  conservation  of 
energy.  Both  of  these  methods  are  dealt  with  in  this  part  of  the  report. 


Cauchy's  Method 


The  response  function  F. , in  Equation  195  can  be  expanded 
as  a polynomial  in  the  strain  tensor  , i.e., 

a. . = a«  + a, e. . + a_e.  e , + a«e.  e e . +...  (196) 

ij  0 ^ ij  2 im  mj  ^ im  mn  nj  ' ^ ' 

where  ag  , a^  , ...  a^  are  real  coefficients.  Utilizing  the  Cayley- 
Hamilton  theorem  we  can  express  Equation  196  in  the  following  form  (see 
Equations  62  and  63} 


“y  “ *o®ij 


* t'ij 


P^im^mj 


(197) 


where  , and  ♦g  elastic  response  coefficients  which  are 

polynomial  fimctions  of  strain  invariants.  Equation  197  is  refeiTed  to 
as  the  Cauchy  elastic  c^  astitutive  equation.  Alternately,  for  an 
elastic  material  we  can  express 
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where  , euid  Vg  are  elastic  response  coefficients  which  are 

polynomial  functions  of  stress  invariants.  From  Equation  197  it  follows 
that  for  isotropic  elastic  materials  the  initial  state  of  stress  is 
hydrostatic,  i.e.-. 


°1J ' Vij  ° ° > 

Also,  using  the  transformation  law  of  a second-order  tensor  (Equa- 
tion 9^)j  it  can  be  shown  that  Equation  197  is  form  invariant  with  re- 
spect to  rigid  motion  of  a spatial  coordinate  system,  i.e., 

o'  = a .a  ,c. . 
mn  mi  nj  ij 

= ^.a  .a  ,6.,  + A^a  .a  .e. , + A^a  .a  ,e,,  e,  , 

^0  mi  nj  ij  ^1  mi  nj  ij  ^2  mi  nj  ik  kj 

* -M'  + (200) 

0 mn  1 mn  2 mk  kn 

where  e'  is  referred  to  the  primed  (rotated)  coordinate  system.  We 
mn 

can  now  utilize  Equation  197  to  develop  various  types  of  isotropic 
elastic  constitutive  equations. 

Linear  elastic  material 

58.  For  linear  elastic  materials  the  response  coefficient 
vanishes.  The  response  coefficient  is  a constant  and  is  a 
linear  function  of  the  first  strain  invariant.  Assuming  that  the  initial 
state  of  stress  is  zero,  the  constitutive  eqiiation  of  linear  elastic  ma- 
terial can  be  written  as 


o, , = AI-5. . + Be. . 
ij  1 ij  10 


(201) 


where  A 6ind  B are  material  constants.  In  order  to  determine  the 
physical  meaning  of  the  materiel  constants  A and  B let  us  consider  a 
simple  shearing  deformation  defined  by 
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° ^12  °1 


=1J  “ V “ 


0 0 0 


For  this  state  of  deformation.  Equation  201  reduces  to 


0 ° 


Oi^  = 0 0 


0 0 oj 


(202) 


(203) 


Since  is  half  the  shearing  strain  (see  Equation  150),  it  follows 

that  B is  two  times  the  shear  modulus  which  we  define  as  G ; thus. 


B = 2G 


(204) 


Next,  we  consider  imiform  dilatation  defined  hy 


^00 


0 3-  0 


(205) 


0 0 3-1 


For  this  state  of  deformation.  Equation  201  becomes  (invoking  Equa- 
tion 204) 


’tj  “ (“i  * “ r)‘ 


(206) 


Taking  the  trace  of  (let  i » J ) we  obtain 


69 


_ '^1  _ . 2GV 


(207) 


Equation  207  relates  pressure  to  volumetric  strain  (l^).  The 

slope  of  the  press\ire-volumetric  strain  relation  is  defined  as  hulk  mod- 
ulus K ; thus. 


A + 


or 


(208) 


(209) 


The  constant  A is  usually  denoted  as  X and  is  referred  to  as  the 
Lame  constant.  In  terms  of  the  shear  and  hulk  moduli,  the  constitutive 
equation  of  linear  elastic  material  (Equation  201 ) then  becomes 


(210) 


The  expression  recognized  as  the  strain  deviation 

tensor  E^^  (Equation  l6l);  thus, 

"ij  = '“ihj  * 

Prom  Equation?  122,  207,  and  208,  it  follows  that  the  constitutive  equa- 
tion of  linear  elastic  materials  can  also  he  written  as 


(212a) 


Kl^  (212h) 

Equation  212  indicates  that  for  linear  elastic  materials  volumetric 
strain  is  caused  hy  hydrostatic  stress  only,  and  that  the  shearing  re- 
sponse of  the  material  is  independent  of  pressure. 
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59*  Using  Equation  212  we  can  readily  express  the  strain  tensor 
in  terms  of  stress  tensor 


E 

ij  2G 


X 

1 3K 


(213a) 


(213b) 


or,  using  Equation  l6l, 


•^1  ®ij 

P = -i.  5 + -iu- 

^ij  9K  °ij  2G 
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60.  We  will  now  proceed  to  examine  the  behavior  of  linear  elastic 
materials  xinder  various  states  of  stress  and  deformation.  Let  us  first 
consider  uniaxiel  state  of  stress,  a common  laboratory  test,  defined  by 


""a  ' 


0 0 


0 0 0 
0 0 0 


(215) 


For  this  state  of  stress.  Equation  2lU  results  in 


/3KjJ_G\ 

\ 9KG  )^1 


(^1 


r2G  - 3KV 

^ idKG  ' Jri 


(216) 


Equation  2l6  indicates  that  under  uniaxial  state  of  stress 


”l  3K  ♦ G *^1 


(217) 
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2 


(218) 


The  ratio  under  uniaxial  state  of  stress  is  referred  to  as 

Yo\ing*s  modulus  E , and  the  ratio  of  radial  strain  to  axial  strain  is 
called  Poisson's  ratio  v ; thus  (for  an  incompressible  elastic  material 

V = 1/2) 


3K  + G 


- 3K  - 2G 
~ 6K  + 2G 


(219) 

(220) 


61.  Another  common  laboratory  test  is  the  uniaxieQ.  strain  test 
defined  by 


0 0 
0 0 0 

0 0 0 

m 


(221) 


For  this  state  of  deformation.  Equation  211  results  in 


(K  + UG/3)e3^ 
0 
0 


0 

(K  - 2G/3)e^ 
0 


0 

0 

(K  - 2G/3)e^ 


(222) 


From  Equations  222,  219,  and  220  it  follows  that  under  uniaxial  strain 
condition 


»1  ' = (I  f '1 

”2  = “3  = (#^)'’i ' (r^)"i 
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It  is  noted  that  Og  is  the  radial  stress  required  to  prevent  radial 
strain.  The  ratio  under  uniaxial  state  of  strain  is  referred  to 

as  the  constrained  modulus  M ; thus, 

“ = “ (1  av)  '"=5) 

Using  Equations  223  and  224  ve  can  determine  an  expression  for  the 
stress  path  associated  wich  the  uniaxial  state  of  strain  in  the  yfTi 
versus  J^/3  space: 

V^=^J,/3  (226) 

In  terms  of  Poisson's  ratio  v , the  equation  of  the  stress  path  hecomes 

VST.  j^/3  (JJ7) 

62.  Next,  let  us  consider  the  behavior  of  linear  elastic  mate- 
rials under  condition  of  pleuae  strain  defined  by 


‘•ij  “ 


^11  ®12 


^21  ^22 


'12 


'22 


For  condition  of  plane  strain.  Equation  211  res’ilts  in 
(iC  + |G)ej^^  2Ge_  0 

- f ) 

♦ (k  . I G)e2j 

0 (K-f)(< 


2Ge: 


21 


‘11 


"U  ^ ^22 


(228) 


(229) 
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Prom  Equation  229  it  follows  that 


^^33  “ (^K  + 2g)^°11  ®22^  ~ '*^'^11  ■*■  °22^ 


(230) 


Equation  230  gives  the  magnitude  of  stress  necessary  to  maintain 

plane  strain  condition,  i.e. , = 0 . 

63-  The  counterpart  of  plane  strain  is  the  condition  of  plane 
stress  defined  by 


‘^ll  °12  ° 


°ij  ■ °21  °22  ° 


0 0 0 


For  plane  stress  condition.  Equation  2lU  results  in 

X !i2 

V9K  30/11  20 

* ^)‘^22 

^21  1 \ 

'^ij  " 20  \9K  ■ 11 

(9K  3g)®22 


(9K  ■ '^)^®11  ®22^ 


Prom  Equation  232  it  follows  that 


“ (9K  ■ ^)^°11  ^2?J  “ “ 1 - V ^^11  * ^22^ 


(231) 


(232) 


(233) 


Equation  233  gives  the  magnitude  of  strain  produced  by  condition 

of  plane  sti*ess,  i.e.,  = 0 . 

61t.  The  constitutive  equations  of  linear  elastic  material  ex- 
pressed in  terms  of  various  combinations  of  elastic  constants  are  given 
in  Table  1 for  ready  use. 


Nonlinear  eleistlc  material 


65.  Constitutive  eq,uations  for  various  classes  of  nonlinear 
elastic  material  can  be  developed  from  the  genersJ.  form  of  the  Cauchy 
elastic  constitutive  equation  (Equations  197  and  198).  Before  we  de- 
velop constitutive  equations  for  various  classes  of  nonlinear  eleistic 
materiail  it  would  be  beneficial  to  examine  the  significance  of  the 
second-order  terms  ^im'^mj  Equations  197  and  198.  Con- 

sider a simple  shearing  deformation  of  amount  2y  defined  by  the  follow- 
ing strain  tensor 


‘0  Y 0 

Y 0 0 
0 0 0 


(231*) 


For  this  state  of  deformation,  Equation  197  results  in  the  following  ex- 
pression for  the  stress  tensor 


10  0 

0 Y 0 

' 2 

Y 0 0 

''ij  “ '*’0 

0 10 

0 0 1 

+ 

Y 0 0 

0 0 0 

^*2 

2 

0 Y 0 

0 0 0 

• a 

. 

(235) 


Froa  Equation  235,  the  shearing  stress  and  shearing  strain  are  related 
by 


0 


12 


♦lY 


(236) 


and  the  normal  stresses  are  given  as 


''ll  = °22 


♦0  V 


(237a) 


(237b) 


Equation  237  indicates  that  to  maintain  a simple  shearing  defonoation 


76 
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(Equation  23^),  normal  stresses  must  te  applied  to  the  houndaries  of  the 
specimen.  Since  two  of  the  normal  stresses  are  unequal.  Equation  237 
predicts  the  occurrence  of  normal  deviatoric  stresses 


^11  “ ^22  " 3 *^2^ 


(238) 


on  the  shearing  planes.  This  is  a direct  consequence  of  the  second- 
order  term  in  Equation  197  and  is  a departure  from  the  linear 

theory  where  (fi^  = 0 . We  now  consider  the  counterpart  of  simple  shear- 
ing deformati.  1,  i.e. , simple  shearing  stress,  and  show  that  Equa- 
tion 193  will  predict  volume  change  for  this  state  of  stress.  Consider 
a simple  shearing  stress  of  amount  r defined  by  the  following  stress 


tensor 


0 T 0 


aij  = T 0 0 


(239) 


0 0 0 


For  this  state  of  stress.  Equation  198  results  in  the  following  expres- 


sion for  the  strain  tensor 


f‘ " i 


0 T 0 


0 o] 


Lo  0 1 


(2U0) 


to  0 0_ 


Lo  0 


From  Equation  2U0  it  follows  that 


^ii  = ^1  = 3^0  ^ 2^2^ 


(?»^1) 


which  indicates  that  simple  shearing  stress  is  accompanied  by  volume 
change.  Also,  from  Equation  2^0  it  follows  that  there  are  normal  devia- 


toric strains 


= =*22  = I ’2"' 


(2U2) 


eissociated  with  the  volume  change.  The  occurrence  of  deviatoric  strains 
is  a direct  consequence  of  the  second-order  term  iii  Equa- 

tion 198.  From  these  two  examples  we  can  conclude  that  in  the  case  of 
nonlinear  elastic  materials  volumetric  strains  are  caused  by  both  the 
hydrostatic  and  shearing  stresses.  Also,  the  shearing  response  of  the 
material  is  dependent  on  the  hydrostatic  state  of  stress.  In  order  to 
further  demonstrate  these  coupiing  effects  we  consider  a combined  state 
of  hydrostatic  and  simple  shearing  stress  given  by  the  following  stress 
tensor 


a 


ij 


P T 0 
T P 0 
OOP 


(2l;3) 


where  P = superimposed  hydrostatic  stress.  For  this  state  of  stress. 
Equation  198  results  in  the  following  expression  for  the  strain  tensor 


10  0 

P T 0 

2 2 

P + T 2Pt  0 

'ij  = ^0 

0 10 

0 0 1 

T P 0 

OOP 

w*  m 

2 2 

2Pt  P*^  + t 0 

2 

OOP 

The  volumetric  strsdu  then  becomes 


(2UU) 


e^i  = 1^  = + 34-^?  + 'i'gCSP^  + 2t^) 


(2U5) 


and  the  shearing  strain  the  form 

^12  ■ '^’l^  ^'^2^ 

Equations  2l*5  and  2it6  once  again  illustrate  the  coupling  which  exists 
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between  the  lydrostatic  stress  (or  volumetric  response)  and  the  shearing 
response  of  nonlinear  elastic  materials. 

66.  We  will  now  proceed  to  develop  constitutive  equations  for 
various  classes  of  nonlinear  elastic  material  within  the  framewc.k  of 
the  Cauchy  elastic  constitutive  equation  (Equations  197  and  I98).  The 
simplest  form  of  nonlinear  elastic  materieil  is  the  second-order  stress- 
strain  relation  where  terms  in  strain  up  to  the  'second  power  sn:*e  re- 
tained in  the  stress-strain  relationship.  To  derive  the  constitutive 
relationship  for  second-order  elastic  material  we  will  start  from  Equa- 
tion 197  and  express  the  response  coefficients  (Ji^  , , and  in 

appropriate  polynomials  of  strain  invariants.  Accordingly,  for  second- 
order  material  we  can  write 

<j>0  = + C2I2  (2h7a.) 

(2U7b) 

4>2  = Cg  (2hTc) 

where  C^  through  are  material  constants  which  must  be  determined 
experimentally.  Since  the  constitutive  equations  of  second-order  elas- 
tic material  must  degenerate  to  the  first-order  equation  (Equation  210) 
if  second  power  terms  in  strain  are  neglected,  the  material  constants 
and  Cj^  should  be  replaced  by  (K  - 2G/3)  and  2G  , respectively. 

The  constitutive  equation  of  second-order  elastic  material  then  becomes 

=13  = [<“  - * C31JSJ, 

. (ao  . Cjijjcy  . (aw) 

67.  Equation  2U8  contains  six  material  constants.  The  physical 
meaning  of  these  constants  and  the  manner  in  which  they  can  be  deter- 
mined from  laboratory  test  resvilts  can  be  doaonstrated  by  examining  the 
behavior  of  second-order  elastic  materials  under  various  states  of 
stress  and  deformation.  Let  us  first  consider  a simple  shearing 


T9 


deformation  defined  by  Equation  202.  For  this  state  of  deformation 
= 0 and  Ig  = Equation  21*8  results  in  the  following  rela- 

tions for  the  components  0:!“  stress  tensor 


’12  “ 

2Geig 

(249a) 

'll  “ 

*^22  “ * ^6^ ^12 

(249b) 

'33  ■ 

^3^12 

(249c) 

Equation  2l*9a  indicates  that  a second-order  elastic  stress-strain  rela-' 
tionship  predicts  a linear  relation  between  shearing  stress  and  shesiring 
strain.  Equations  2l*9b  and  249c  give  the  magnit  de  of  normal  stresses, 
as  a function  of  shearing  strain,  required  to  maintain  shearing  deforma- 
tion. It  is  noted  that  the  normal  stresses  are  not  uniform,  thus  re- 
sulting in  normal  deviatoric  stresses 


S 


1 

3 


2 

12 


(250) 


on  the  shearing  planes.  The  significance  of  the  material  .onstants 
and  the  combination  (C^  + C^)  is  realized  f;-om  Equations  249c  and  249b, 
respectively. 

68.  We  next  consider  uniform  dilatation  definfd  by  Equation  205. 
For  this  state  of  deformation  = 1^/6  and  Equation  248  results  in 
the  following  relationship  between  pressure  and  volumetric  strain 

•^  = KI^  + + I j ^ (251) 

Equation  251  describes  a parabolic  stress-strain  relationship.  The 
matei'ial  constants  K and  the  combination  ^Cg  + + C^/3  + Cg/9^ 

can  be  determined  from  experimental  data  by  curve  ficting  techniques. 

It  is  noted  that  K is  the  initial  slope  of  the  pressure-volumetric 
strain  curve  ail  is  a positive  constant.  If  the  combination  of  the  ma- 
terial constants  in  the  parentheses  is  also  positive,  the  stress-strain 
curve  will  be  concave  to  the  stress  axis.  If,  on  the  other  hand,  this 
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combination  is  negative,  the  stress-strain  curve  vill  be  concave  to  the 
strain  axis. 

69.  We  will  next  consider  uniaxial  state  of  strain  defined  by 
Equation  221.  For  this  state  of  deformation  and  - e^/2 

and  equation  2il8  results  in 

0^  = (k.  I o)c^  . (c^  . i Cj  . C5  . 

.,  = .3  = (K-|aX.(c,4)c^ 

Equations  252  and  253  also  describe  parabolic  stress-strain  relation- 
ships. The  combinations  of  material  constants  + 0^/2  + 

and  (C^  + ^2/2)  can  be  determined  from  experimental  data  by  curve  fit- 
ting techniques.  As  was  pointed  out  previously,  the  shape  of  the  stress- 
strain  curves  predic  by  Equations  "^52  and  253  depends  on  the  sign  of 
the  combination  of  the  material  constants  in  the  parentheses.  The 
stress-strain  ciirves  will  be  concave  to  the  stress  aids  if  these  combina- 
tions are  positive,  and  concave  to  the  strain  axis  if  they  axe  negative. 
Since  is  the  only  nonvanishing  strain  component  in  uniaxial  strain 
configuration,  we  can  use  Equations  252  and  253  to  relate  stress  dif- 
ference - Og  'to  strain  difference  , i . e . , 

- Og  = 2Ge^  + (C^  + Cg)e^  (25^*) 

Again,  the  shape  of  the  stress-strain  curve  is  determined  from  the 
sign  of  the  combination  (C^  + C^)  of  the  material  constants. 

TO.  More  complicated  states  of  stress  and  deformation,  such  as 
plane  strain  and  triaxial  stress  conditions,  can  also  be  studied  within 
the  framework  of  second-order  stress-strain  law.  Such  states  of  stress 
generally  lead  to  lengthy  mathematical  expressions  between  stress  and 
strain  components.  For  example,  consider  cylindrical  state  of  strain 
defined  by 


(252) 

(253) 
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0 0 


0 eg  0 


0 0 e. 


^255) 


For  this  state  of  deformation  = e^  + 2e2  and  Ig  = 1/2 + SCg^ 
and  Equation  2kd  results  in  the  following  expressions  for  the  components 
of  stress  tensor 

* j °)=i  * (=2  * r * S * - 1 <^h 

* * ^3)^2  ■*■  (^^2 


02  = 03 


/r  2 , 

,\  / ^3\  2 

2G\ 

V"l  (^2  2 j^l 

(2K  . -} 

(l.c. 


+ I 4Cg  + + 2C^  + Cg 


;)4  * ("S  * =5>i^2 


(257) 


Various  stress  paths  may  he  employed  in  a lahoratoiy  test  maintaining  a 
cylindrical  state  of  strain.  The  most  common  stress  path  used  with  this 
state  of  strain  is  to  keep  the  lateral  stress  constant  while  in- 

creasing (Figure  6d).  For  this  stress  path  it  is  possible,  in  prin- 
ciple, to  solve  for  Eg  in  term  of  , using  Equation  257  (since  o^ 
is  a constant),  and  then  substitute  the  resulting  expression  into  Equa- 
tion 256  to  develop  a relationship  between  and  . Other  stress 
paths  such  as  constant  J^/3  path  (Figure  6c)  and  uniaxial  stress  test 
(Figure  6b)  can  also  be  considered. 

71.  i'jllowing  the  same  procedure  we  can  develop  and  analyze  more 
complicated  nonlinear  elastic  stress-sti'ain  laws.  Let  us  consider  a 
third-order  law  where  terms  in  strain  up  '■  the  third  power  are  retained 
in  the  stress-streiin  relationships.  Accordingly,  the  response  coeffi- 
cients iJiq  , , and  ({ig  (Equation  1.97)  for  a third-order  stress- 

strain  law  take  the  following  forms 


(258a) 


♦o  = (■'  - 1 " S'a  " v!  * VA  * V: 


(j»  = 2G  + C.I 


S'"!  '*’  ^10^1 


(258b) 


^2  = "6  "l2'^l 


(2‘^8c) 


where  Cj  through  additional  material  constants  which 

must  he  determined  experimentally.  A third-order  elastic  stress-strain 
law  fo emulated  within  the  framework  of  Equation  197  (Cauchy’s  method), 
therefore,  contains  twelve  material  constants.  In  view  of  Equation  197 
the  constitutive  equation  of  third-order  elastic  material  becomes 


"ij  “[(''- 1 “) 

+ ^20 


h " =2^1  * hh  * v!  * VA  * <=9^1  "y 


3]®!, 


. C,I,  . C,„ir  * . (C,  . C„I,  k,_.. 
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"10  1 


12n'  im  mj 


(259) 


It  is  noted  that  if  third-order  terms  in  strain  are  neglected,  Equa- 
tion 259  reduces  to  Equation  248  (constitutive  equation  of  second-order 
elastic  materials ) . It  was  pointed  out  previously  that  a second-order 
elastic  stress-strain  relationship  predicts  a linear  relation  between 
shearing  stress  and  shearing  strain.  Nonlinear  relation  between  shear- 
ing stress  and  shearing  strain  is  due  to  third-  or  higher-order  terms  in 
strain  tensor.  This  phenomenon,  which  is  a departure  from  second-order 
effect  (see  Equation  2h9a),  can  be  demonstrated  by  examining  the  behav- 
ior of  third-order  elastic  materials  under  simple  shearing  deformation 

defined  by  Equation  202,  For  this  state  of  deformation  1=0, 

— 2 — 

Ig  = > and  ~ 0 Equation  259  results  in  the  following  re- 


lation for  the  shearing  stress  a 


12 


^^12  ^ ^^^12  ^ll'^12 


(260) 


Equation  26O  is  a third-order  equation  in  shearing  strain  e 


12 


The 


behavior  of  third-order  stress-strain  law  under  various  states  of  stress 
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and  deformation  can  al*'^'*  be  studied  similar  to  the  second-order  law.  In 
the  next  section  we  will  consider  less  complicated  and  perhaps  more 
useful  forms  of  elastic  stress-strain  laws  referred  to  as  quasi-linear 
elastic  material. 

^uasi-linear  elastic  material 

72.  Nonlinear  elastic  stress-strain  laws  are  too  complicated  for 
application  in  all  engineering  problems.  In  many  engineering  problems 
only  an  approximate  or  gross  behavior  of  the  material  under  considera- 
tion needs  to  be  modeled.  For  this  reason,  we  will  develop  a number  of 
simple  sur ess-strain  laws  for  simulating  the  gross  behavior  of  a nianber 
of  materials.  We  will  start  with  Equation  197  (or  Equation  198)  by 
making  the  assumption  that  the  response  coefficient  is  zero.  The 
basic  constitutive  equation  then  becomes 


a 


♦bij 


(261) 


where,  as  before,  <{|q  and  (j>^  are  polynomial  functions  of  strain  in- 
variants. Equation  26l  is  usually  called  a quasi-line<sir  relation.  From 
Equation  26l  it  follows  that 


1 _ . ^ 1 . T 

3 " ^0  ^ 3 ^1^1 


(262a) 


^I^iJ 


(262b) 


In  view  of  Equations  12l*,  l62,  and  262  we  can  write 


♦o 
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1 
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(263a) 


(263b) 
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Substituting  Equation  263  in  Equation  26l  results  in  the  following  gen- 
eral constitutive  equation  for  quasi-linear  elastic  material 


ij  3 ij 


(264) 


It  is  only  necesseiry  to  postulate  (based  on  experimental  evidence)  math- 
ematical expressions  for  and  Jg  in  terms  of  strain  invariants  in 
order  to  utilize  Equation  264  for  any  material  of  interest.  The  inverse 
of  Equation  264,  resulting  in  strain-stress  law,  can  be  obtained  from 
Equation  198  by  assuming  that  S',  is  zero  and  following  the  above  pro- 

' C. 

eedure.  The  resiUting  relationship  becomes 


e 


iJ 


+ 


(265) 


To  use  Equation  265  we  need  to  express  and  in  terms  of  stress 
invariants.  For  example,  in  the  case  of  the  linear  elastic  materialr- 


(266a) 

(266b) 


and  Equations  264  and  265  reduce  to  constitutive  equations  of  linear 
elastic  materials  (Equations  210  and  2l4,  respectively). 

T3.  We  will  now  proceed  to  develop  constitutive  equations  for 
various  classes  of  quasi-linear  elastic  materials  which  are  of  interest 
for  engineering  application.  For  the  simplest  class  of  quasi-linear 
elastic  materied.,  in  which  there  are  no  couplings  between  the  dsviatoric 
and  volumetric  responses  of  the  material,  we  can  write 


(267a) 

(267b) 
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where  the  functions  f^  and  fg  must  he  determined  based  on  experi- 
mental evidence.  A good  approximation  for  a number  of  materials  (such 
as  clay  soils)  is  to  assume  the  following  relations  for  f^  and  fg 

J al 

r = (fo  * '‘b>®  - "’b 


(268b) 


where  P„  , a , , and  are  material  constants  which  must  be  de- 

termined  experimentally.  Equations  268a  and  268b  are  depicted  graphi- 
cally in  Figvire  9*  It  is  observed  from  Figure  9 that  defines  an 

initial  hydrostatic  state  of  si  ess  (for  materials  that  can  sustain  ten- 
sion, Pq  can  be  taken  to  be  zero)  and  P^  defines  the  maximum  hydro- 
static tensile  stress  that  the  material  can  sustain  before  it  fails 
(breaks)  at  such  a tension.  For  materials  that  cannot  sustain  tension, 

the  material  constant  F„  is  zero.  In  this  case  P_  defines  the  state 

B 0 

of  "ease”  or  the  initial  stress  state  of  the  material.  The  material 
constant  k^  is  proportional  to  the  inverse  of  the  initial  shear  modu- 
lus and  kg  is  the  inverse  of  the  ultimate  shear  strength  of  the  mate- 
rial. Substituting  Equations  268a  and  268b  in  Equation  26k  results  in 
the  following  quasi-linear  elastic  constitutive  equation 


al. 


- P. 


ijj 


(269) 


The  inverse  of  Equation  269  resulting  in  strain-stress  law  can  be  ob- 
tained by  inverting  Equations  268a  and  268b  said  substituting  the  result- 
ing expressions  for  I^  and  y[^  into  Equation  265  as  follows: 


1 , * P3' 


*1J  * 
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(270) 
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8.  PRESSURE-VOLUIIETRICSIllMli  RESPONSE  (cQ.268i) 


7^.  Equation  269  (or  Equation  270)  is  a simple  but  useful  consti- 
tutive equation  which  can  be  used  to  study  the  stress-strain  behavior  of 
a number  of  physically  nonlinear  materials.  It  contains  only  five  mate- 
rial constants  which  have  physical  meaning  and  can  easily  be  determined 
experimentally.  Having  determined  the  nvnnerical  values  of  these  mate- 
rial constants  we  can  use  Equations  269  and  270  to  predict  the  behavior 
of  the  material  xander  any  state  of  stress  and  deformation.  For  example, 
consider  condition  of  uniaxial  strain  defined  by  Equation  221.  For  this 
state  of  deformation  and  Ig  = e^/3  » and  Equation  269  results 

in 


oe. 


a,  = 


2 

3 "1 


^2 

k.  + — e- 


(271) 


ae. 


02  = O3  = (Pq  + Pg)e  - Pfi  " 
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3 ^ 


k-  + — e, 
1 ^ 1 


(272) 


Equations  271  and  272  predict  stress-strain  curves  that  may  initially  be 
concave  to  the  strain  axis  and  then  become  concave  to  the  stress  axis  as 
the  vertical  strain  increases.  Using  Equations  271  and  272  we  can 
determine  the  stress  path  associated  with  the  state  of  uniaxial  strain 
in  the  versus  J^/3  space,  i.e.. 


(273) 


Next  consider  the  behavior  of  the  material  \inder  condition  of  uniaxial 
stress  defined  by  Equation  215  (assuming  that  Pq  = O).  For  this  state 
of  stress  = o^  and  = o^/3  and  Equation  270  results  in 
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75*  A more  complicated  quasi-linear  elastic  material  model  can 
now  be  constructed  by  assuming  that  the  sheaoring  response  of  the  mate- 
rial is  a function  of  both  the  hydrostatic  and  the  deviatoric  stresses, 
while  the  volumetric  response  is  only  a function  of  pressure.  Accord- 
ingly, Equation  267a  is  still  valid  while  instead  of  Equation  267b  we 


can  write 


(276) 


Various  forms  of  Equation  276  can  be  utilized  to  construct  a material 
model.  A useful  form  for  materials  such  as  sand  can  be  developed  by 
using  Equation  268b  emd  assuming  that  the  ultimate  shear  strength  of  the 
material  is  a function  of  . For  a first-order  approximation  we 

can  assume  that  the  ultimate  shear  strength  is  a linear  function  of 
hydrostatic  stress;  thus. 


kg  ^2  ^3  3 


(277) 


where  kg  and  k^  are  material,  constants  that  must  be  determined  ex- 
perimentally. Utilizing  Equations  277,  268a,  and  268b  in  Equation  26U 
results  in  the  following  constitutive  relationship 
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It  is  noted  that  when  dependence  of  shear  strength  on  hydrostatic  stres 
disappears  (i.e. , when  = O)  Equation  278  reduces  to  Equation  269. 

76.  To  examine  the  significance  of  the  dependency  of  shear 
strength  on  hydrostatic  stress  let  us  consider  the  behavior  of  Equa- 
tion 278  xmder  cylindrical  state  of  strain  ^Equation  255).  For  this 
state  of  deformation  VII  = (e^  - and 

can,  after  arranging  terms,  obtain  the  following  relationship 
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(279) 


If  we  consider  a stress  path  where  is  kept  constant  during  the 

test  (Figure  6d),  Equation  279  -an  be  used  to  relate  stress  difference 
(o^  - Og)  to  strain  difference  (e^  - e^)  for  a constant  value  of  Og  . 

The  qualitative  behavior  of  Equation  279  is  depicted  in  Figure  10.  It 
is  obsei^ed  from  Figure  10  that  the  shear  strength  of  the  material  in- 
creases with  increasing  confining  stress  Og  . If  k^  is  set  to  zero 
in  Equation  279  it  is  noted  that  the  dependency  of  shear  strength  on  con- 
fining stress  disappears. 

77*  The  next  step  in  developing  more  complicated  quasi-linear 
elastic  stress-strain  relationships  is  to  eissume  that  volumetric  strain 
I^  is  caused  by  both  the  hydrostatic  and  deviatoric  stresses.  For  such 
material  we  can  write 
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Figure  10.  Response  predicted  by  Equation  279  for  cylindrical 

state  of  strain 


where  the  function  g must  be  postulated,  based  on  experimental  re- 
sults, for  any  material  of-  interest.  Equation  280  csin  be  simplified 
further  by  expressing  the  volumetric  strain  as  the  sum  of  two  com- 
ponents, i.e. , 


(28l) 


where  g,  = contribution  due  to  hydrostatic  stress  and  g = contribu- 
n s 

tion  due  to  deviatoric  stresses.  For  the  contribution  due  to  hydro- 
static stress  we  can  use  the  inverse  of  Equation  268a 


i.P 

3 ^ 


(282) 


As  a first-order  approximation,  for  the  contribution  due  to  deviatoric 
stresses  we  can  e3q>ress  g as 


gg  = aV^ 


(283) 


where  o is  a materieil  constant  that  can  he  positive  or  negative  de- 
pending on  whether  the  material  contracts  or  expeinds,  respectively,  dur- 
ing the  application  of  deviatoric  stresses.  In  view  of  Equations  28l 
through  283  the  relationship  for  the  volumetric  strain  (Equation  280) 
becomes 

+ aV^  (281;) 

During  a hydrostatic  test  (Figure  6a)  Jg  = 0 and  Equation  281;  reduces 
to  Equation  282.  To  formxilate  the  constitutive  equation  for  this  class 
of  quasi-linear  material  we  combine  Equations  268b  and  277  to  develop  an 
ej^ression  for  then  substitute  this  expression  and  Equa- 

tion 284  in  Equation  265*  The  resulting  constitutive  equation  becomes 

m 

+ aV^ 


78.  Equation  285  allows  for  the  dependency  of  sheer  strength  on 
hydrostatic  stress  and  the  coupling  of  volumetric  strain  and  deviatoric 
stresses.  If  the  material  constants  and  a are  set  to  zero  these 
cross-effects  will  disappear  and  Equation  285  will  reduce  to  Equa- 
tion 270.  It  sho\ild  be  noted  that  these  cross-effects  (in  particular 
the  coupling  of  volxmetric  strain  and  deviatoric  stresses)  are  due  to 
scalar  nonlinearity  (invr.riants)  and  are  different  from  the  second-order 
effects  discussed  in  the  development  of  constitutive  equations  for  non- 
linear elastic  materials.  To  illustrate  this  point  further,  let  us 
examine  the  behavior  of  Equation  285  xmder  a simple  shearing  stress 
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defined  by  Equation  239.  For  this  state  of  stress  (taking  P„  to  be 

2 ^ 

zero)  J^/3  = 0 and  ^2  “ ^ Equation  285  it  follows  that 


e . . =1,  = at 
11  1 


(286) 


indicating  that  simple  shearing  stress  is  accompanied  by  volume  change. 
However,  Equation  285  indicates  that  there  are  no  normal  deviatoric 


strain  components  associated  with  the  volume  change,  i.e.,  E^^  = E^ 


= = 0 . In  the  case  of  nonlinear  elastic  material,  on  the  other 


hand,  it  was  shown  that  there  are  normal  deviatoric  strains  associated 
with  volume  change  due  to  simple  sheeiring  stress  (see  Equation  2U2). 
Finally,  it  should  be  pointed  out  that  for  certain  stress  paths  Equa- 
tions 278  eind  285  will  produce  inelastic  stress-strain  response  (these 
stress  paths  will  be  discussed  in  Part  V in  conjunction  with  incremental 
stress-strain  laws ) . This  is  a consequence  of  the  dependency  of  the 
shearing  response  of  the  material  on  the  hydrostatic  state  of  stress 
(Equation  276).  To  avoid  such  possibility  the  use  of  Equations  278 
and  285  should  be  restricted  to  stress  paths  where  J^/3  and  do 
not  decrease,  i.e.,  they  remain  9QU3tant  or  increase. 


Green's  Method 


79.  In  order  to  develop  the  constitutive  equations  of  elastic  ma- 


terials based  on  Green's  method  we  first  state  two  of  the  fundamental 
laws  of  mechanics:^ 


The  first  law  of  thermodynamics,  which  is  a statement  of 
the  law  of  conservation  of  energy:  "The  work  that  is 
performed  on  a mechanics.:  system  by  external  forces  plus 
the  heat  that  flows  into  the  system  from  the  outside 
equals  the  increase  of  kinetic  energjr  plus  the  increase 
of  internal  energy." 


b.  The  law  of  kinetic  energy:  "The  work  of  all  the  forces 
(internal  and  external)  that  act  on  a system  equals  the 
increase  of  kinetic  energy  ol  the  system." 


Mathanatically,  we  can  express  the  first  law  of  thermodynamics  as' 


+ Q = AT  + AU 


(287) 
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where 


Wp  = work  performed  on  the  system  hy  external  forces 
Q = heat  that  flows  into  the  system 
AT  = increase  of  kinetic  energy 
AU  = increase  of  internal  energy 
The  law  of  kinetic  energy  can  be  stated  as 

W + VJ.  = AT  (288) 

e 1 

where  = the  work  performed  by  internal  forces  in  the  system.  In 
view  of  Equations  28?  end  288  we  can  write 

= Q - AU  (289) 

Since  we  are  only  dealing  with  mechanical  energy  weassxame  that  Q = 0 
and  Equation  289  reduces  to 


W.  = -AU 
1 


(290) 


80.  We  will  new  proceed  to  derive  the  constitutive  equation  of 
Green’s  elastic  materials.  If  V is  a material  volume  (region)  within 
a deformaole  body  and  S is  the  surface  enclosing  this  region,  and  if 
this  region  undergoes  an  infinitesimal  displacement  (the  symbol  6 

defines  a small  variation) , the  work  of  external  forces  can  be  expressed 
as 


6W 


" /£  * III,  ^ 


^6u^  dV 


(291) 


where  n^  = direction  cosines  of  the  outward  normal  to  surface  S . The 

first  integral  in  Equation  291  is  due  to  tractive  forces  on  S and  can 

5 

be  transformed  to  a volume  integral  by  using  the  Divergence  Theorem, 

i . e . , 

IL  ’’jl-'j®"!  = //I  * ///„ 


The  second  integral  in  Equation  291  is  due  to  body  forces  acting  on 
material  in  V . Combining  Equations  292  and  291,  the  work  of  external 
forces  becomes 


6W 


Now,  let  us  assume  that  during  the  displacement  6u^  the  material  vol- 
ume V is  in  equilibrium  and  the  change  in  kinetic  energy'  is  zero.  The 
equations  of  motion  (Equation  192)  then  take  the  form 


“UJ  * '■l  “ ° 


(29^+) 


which  are  referred  to  as  the  equations  of  equilibrium.  Since  the  stress 
tensor  is  symmetrical,  application  of  Equations  29it  in  Equation  293  re- 
sults in  the  following  expression  for 


6W 


= //i  ^ 


(295) 


Similar  to  Equations  11^9  throxigh  151,  the  infinitesimal  displacement 

gradient  (6u. ) . can  be  expressed  as 
1 »J 
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(296) 


In  view  of  Equation  296,  Equation  295  becomes 
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Since  ® (symmetry  of  stress  tensor)  and  is  skew- 

symmetric,  the  expression  is  zero  and  Equation  297  reduces  to 
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= ///, 
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(298) 


Invoking  the  assumption  that  AT  = 0 during  the  infinitesimal  displace- 
ment, and  since  Q = 0 , Equation  28?  can  be  written  as 


6W  * 6U  (299) 

e 

The  internal  energy  associated  vith  the  material  volume  V can  be  ex- 
pressed as 


U 


dV 


(300) 


where  Uq  is  the  internal  energy  per  unit  volume,  referred  to  as  the 
internal  energy  density.  In  view  of  Equations  298,  299 » and  300,  we  can 
write 


III 


dV 


V ir 


-III 


6U^  dV 


(301) 


which  leads  to 


(302) 


Since  the  internal  energy  density  function  depends  on  the  strain 
components  , the  variation  6Uq  due  to  6e^j  can  be  expressed  as 

% = ^ ‘=13  '3«) 

In  view  of  Equations  302  snd  303,  the  stress  tensor  takes  the  fol- 

lowing form 
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Equation  30it  is  referred  to  as  the  Green  elastic  constitutive  equation. 

8l.  For  isotropic  materials  the  strain  energy  function  Uq  mxxst 
be  invariant  and,  thus,  a function  of  strain  invariants.  Therefore,  for 
isotropic  materials  we  can  write 
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If  the  material  under  consideration  is  incompressible,  a 0 


and 


.)• 


For  such  material,  an  arbitrary  hydrostatic  state 


C ’ ^3) 

of  stress  may  be  superimposed  on  the  existing  state  of  stress  given  by 
Equation  30i*.  Using  chain  rule  of  differentiation.  Equation  305  can  be 
expressed  as 


31^  3e.^  3e.j  3e^j 


(306) 


Since 
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Equation  306  can  be  written  as 
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Comparison  of  Equation  308  with  Equation  197  indicates  that  the  Green 
and  the  Cauchy  elastic  constitutive  equations  have  the  same  form.  The 
difference  between  the  two  formulations  is  that  the  response  coeffi- 
cients 6q  , (}>^  , and  4*2  Equation  197  are  independent  whereas  the 
corresponding  response  coefficients  , bUQ/SIg  » and  BUq/SI^ 

in  Equation  308  are  not.  By  differentiating  the  response  coefficients 
in  Equation  308  with  respect  to  strain  invariants,  it  follows  that  the 
following  relationships  exist  between  these  coefficients 


97 
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(309a) 


(309b) 


(309c) 


The  consequence  of  the  above  restrictions  imposed  on  the  response  coeffi- 
cients will  be  realized  when  we  develop  second-  and  higher-order  elastic 
stress-strain  laws  using  Equation  308.  The  Green  elastic  material  can, 
therefore,  be  considered  as  a special  type  of  Cauchy  elastic  material 
where  the  response  coefficients  are  restricted  by  Equation  309. 

82.  The  inverse  of  Equation  308  (the  counterpart  of  Equation  198) 
can  be  determined  by  assuming  that  there  exists  a function  so  that 


"o  * ’’o  ° 


(310) 


Equation  310  holds  for  elastic  materials  where  application  of  a positive 
stress  increment  results  in  a positive  strain  increment  and  vice  versa. 
The  function  Tq  is  referred  to  as  the  complementary  energy  density 
function.  From  Equation  310  it  follows  that 
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Differentiating  Equation  311  with  respect  to  a yields 
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Since  Uq  is  a function  of  strain  it  follows  that 
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Comhining  Equations  313  and  312  results  in 
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In  view  of  Equation  30^  the  second  expression  in  Equation  3lH  is  zero. 
Equation  31^  then  becomes 
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Since 
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mn 


Equation  315  becomes 


i m or  J n 


a . j 


(316) 


mn  vo 


nttt 


(317) 


or 
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Equation  318  is  the  inverse  form  of  Equation  30U. 

83.  For  isotropic  materials  is  a function  of  stress  invari- 
ants given  by  Equations  118  through  120,  i.e. , 
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(319) 


Using  chain  rule  of  differentiation.  Equation  319  can  be  expressed  as 


Since 
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Equation  320  takes  the  following  form 
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(321a) 


(321b) 


(321c) 


(322) 


Equation  322  has  the  same  f^rm  as  Equation  198 • The  response  coeffi- 
cients in  Equation  322,  however,  are  not  independent.  It  can  readily  be 
shown  that  relations  similar  to  Equation  309  exist  between  these 
coefficients. 

8U.  The  complementary  energy  density  function  can  edso  be 
expressed  in  terms  of  strain  invariants  by  utilizing  Equations  308  and 
310.  Prom  Equations  308  and  310  it  follows  that 
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(323) 


Since  for  a given  material  Uq  is  known,  we  can  use  Equation  323  to  ex- 
press the  complementary  energy  density  function  in  terms  of  strain  in- 
variants. However,  to  obtain  an  inverse  constitutive  relationship 
(Equation  322)  we  need  to  express  in  terms  of  stress  invariants. 

This  can  be  accomplished,  at  least  in  principle,  by  first  expressing  the 
stress  invariants  in  terms  of  strain  invariants  using  Equation  308.^ 

The  resulting  expressions  can  be  inverted  to  obtain  strain  invariants  in 
terms  of  stress  invariants  and  then  substituted  in  Equation  323  in  order 
to  express  Tq  in  terms  of  * *^2  * *^3  * will  now  develop 

constitutive  equations  for  various  classes  of  elastic  materials  utiliz- 
ing Green’s  method  (i.e..  Equations  308  and  322), 

Linear  elastic  material 

85*  For  linear  elastic  materials  only  terms  in  strain  up  to  the 
first  power  are  retained  in  the  stress-strain  relationship.  It  then 
follows  from  Equation  308  that  the  strain  energy  density  function 
must  be  quadratic  in  strain  in  order  for  the  resulting  stress-strain  re- 
lation to  be  of  first  order.  Assuming  that  the  initial  state  of  stress 
is  zero,  Uq  takes  the  following  form 


^0  = ^1^2  * ^2^ 


(324) 


where  and  A^  are  material  constants.  Substituting  Equation  324 
in  Equation  308  results  in 


®ij  ^2^1^iJ  * ^l^ij 


(325) 
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Eqviation  325  is  identical  with  Equation  201,  indicating  that  hoth 
Green's  method  and  Cauchy's  method  result  in  the  same  stress-strain  rela- 
tion for  linear  elastic  materials.  In  view  of  Equations  20^*  and  209  the 
material  constants  and  become 

= I - I (326a) 


= 2G 


(326b) 


The  strain  energy  density  function  for  linear  elastic  materials  (Equa- 
tion 32U)  can  then  be  expressed  in  terms  of  shear  modulus  G and  bulk 
modulus  K as  follows: 


"o  = “ ('2  - 5 ^1)  * I 4 '327) 

It  is  noted  that  the  expression  - 1^/6^  is  the  second  invariant 
of  strain  deviation  tensor  1^  (Equation  l61»).  We  can  therefore  write 

Uq  = 2GI’  + I I^  (323) 

Since  for  linear  elastic  materials  J^/3  = KI^  and  = 2G^I^  (Equa- 
tion 266),  the  strain  energy  density  function  can  eilso  be  written  as 

Uo=>^V^  + |^ii  (329) 

where  and  J^I^/6  can  be  considered  as  energy  due  to  distor- 

tion and  volume  change,  respectively,  during  a deformation  process.  The 
strain  energy  density  f<inction  can  also  be  expressed  in  terms  of  stress 
invariants  or  various  ether  combinations  of  elastic  moduli  end  invari- 
ants. For  example,  the  counterpart  of  Equation  320  becomes 


(330) 
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In  view  of  Eqiiations  323  and  327  the  complementary  and  strain  energy 
density  functions  for  linear  elastic  materials  are  identical.  For 
linear  elastic  materials  we  can,  therefore,  write 


“o  = '■o  = 5 

86.  We  can  now  examine  the  nature  of  the  restrictions  that  must 
he  placed  on  the  linear  elastic  moduli  due  to  the  existence  of  strain 
energy  density  function  Uq  . Expanding  Equation  32^*,  we  can  write 

“o  = (f  * * 4 * 4)  * ^(4  * 'i3  * 4) 


^2(^11*22  ^22^33  * ^11^33) 

Equation  332  is  queuiratic  in  strain  and  can  he  ej^ressed  in  the  qua- 
dratic form 

6 6 
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^ i=l  j=l  ^ ^ 

where  and  €j  denote  the  six  independent  components  of  the  strain 
tensor.  The  matrix  c^j  ~ ^Ji  expressed  in  terms  of  and  Ag 
and  has  the  form 


2 + A2  Ag 


^2 


0 0 


0 0 


c.. « r2 


^ + A 0 0 


The  strain  energy  density  funcnion  Uq  is  a positive  qusntity.  Accord- 
ing to  the  theory  of  quadratics,  in  order  for  to  be  a positive  quan- 
tit.'y  al).  the  mi;  of  the  ^yeigonal  elements  of  c^j  must  be  iwsitive. 
Imposing  \hiv'  restriction  0‘.  tx*.,  -iiagonal  minors  of  the  matrix  c^j 
leads  to  the  following  in^  >iut.  If  ties: 

A,  > 0 (335a) 


+ > 0 


(335b) 


^ + 3A2  > 0 


(335c) 


In  view  of  Equations  326,  219,  and  220,  the  above  inequalities  impose 
the  following  restrictions  on  the  linear  elastic  moduli  G , K , E , 
and  V : 


G > 0 


(336a) 


K > 0 


(336b) 


E > 0 


(336c) 


-1  < V < I 


(336d) 


It  is  emphasized  that  negative  values  of  Poisson's  ratio  v have  not 
been  found  experimentally  for  isotropic  elastic  materials. 

Nonlinear  elastic  material 

07.  We  will  now  proceed  to  develop  constitutive  equations  for 
various  classes  of  nonlinear  elastic  materieds  within  the  framework  of 
Green's  method  (Equations  300  and  322).  As  in  the  ^ ^velopment  of  non- 
linear Cauchy  elastic  constitutive  equations,  we  will  t>cart  with  second- 
order  stress-strain  relation  where  terms  up  to  the  second  power  eu^  re- 
tained in  the  stress-strain  relationship.  The  strain  energy  density 
function  for  a second-order  stress-strain  law,  therefore,  must  be  cubic 


in  strain.  Assigning  that  the  initial  state  of  stress  is  zero,  the  most 
general  cubic  relation  for  takes  the  following  form 

Uq  = 2GI2  + (f  - 1)1^  + ^ ^U^l  ^ Vs 

where  A^  through  A^.  are  additional  material  constants  associated 
with  second  power  terms  in  the  stress-strain  relationship.  It  is  noted 
that  if  third  power  terms  in  strain  are  neglected  Equation  337  degen- 
erates to  the  corresponding  expression  for  linear  elastic  material  (see 
Equations  32U  and  326).  Substituting  Equation  337  in  Equation  308  we 
obtain  the  following  second-order  stress-strain  law 

. (20  . AjIJcy  . (338) 

Equation  338  contains  five  material  constants  whereas  its  counterpart 
based  on  Cauchy's  method  (Equation  2U8)  contains  six.  This  reduction  in 
materied  constants  is  a consequence  of  thermodynamic  restrictions  im- 
I)osed  on  the  response  coefficients  of  Green  elastic  material  (see  Equa- 
tion 309).  It  was  noted  that  in  the  case  of  linear  elastic  materials 
both  methods  resulted  in  the  same  equation  and  their  difference  was  not 
apparent.  The  difference  between  these  two  methods  becomes  more  pro- 
nounced when  considering  higher-order  nonlinear  elastic  and  quasi-linear 
elastic  materials. 

88.  Let  us  now  consider  the  stress-s* rain  relationship  for  thii^i- 
order  elastic  materials.  The  strain  energy  density  function  for  a 
third-order  elastic  material  must  contain  strain  terms  up  to  the  fourth 
power.  Again  assuming  that  the  initial  state  of  stress  is  zero,  the 
strain  energy  density  function  for  a third-order  elastic  material  takes 
the  following  representation 

= 2GT^  ♦ ( I - IV?  ♦ A,I,T„  ♦ A,  I?  A^T,  •*  A,I? 

0 2V2  3/1  312  “1  53  61 

+ A^l|  + A0I1T3  + A^I^g  (339) 
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where  Ag  through  are  additional  material  constants.  It  is  noted 
that  when  fourth  power  terms  in  strain  are  neglected  Equation  339  re- 
duces to  Equation  337  (strain  energy  function  for  second-order  material) 
Substituting  Equation  339  in  Equation  308,  we  obtain  the  following  third 
order  stress-strain  law 


= [(k  - I o)l^  . . l.Agl3  . Ajij  . 2AjI^l3]  Sy 

. (20  . AjI^  . OA^Ij  . A,l2)cy  * (Aj  . (3k0) 


Equation  3^0  contains  nine  material  constants.  The  counterpart  of  Equa- 
tion 3^0  based  on  Cauchy's  method  (Equation  259) » on  the  other  hand, 
contains  12  material  constants.  Therefore,  in  the  case  of  third-order 
stress-strain  law  the  effect  of  thermodynamic  restrictions  is  to  reduce 
the  n’lmber  of  material  constants  by  three.  In  the  next  section  we  will 
consider  quasi-linear  elastic  materials  within  the  framework  of  Green's 
method. 

Quasi-linear  elastic  materials 

89.  If  the  strain  energy  density  function  Uq  is  independent  of 
the  third  strain  invariant  I^  > Equation  308  reduces  to 


■ ai^  ^ij  ^ij 


(31+1) 


Ec.uation  3^1  is  the  counterpart  of  Equation  26l  (constitutive  equation 
of  quasi-linear  elastic  materials  based  on  Cauchy's  method).  In  the 
case  of  Equation  3^1 > however,  the  response  coefficients  are  restricted 
by  Equation  309a.  From  Equation  3^1  it  follows  that 


3U. 


31, 


(3J*2a) 


(3i*2b) 
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Using  Equations  12U,  l62,  and  3^2,  we  can  obtain  expressions  similar  to 
Equation  263  for  the  response  coefficients  in  Equation  31^1,  i.e.. 


(3^3a) 


(3l+3b) 


The  inverse  of  Equation  leading  to  strain-stress  relationship, 
takes  the  form 


!!o  , !!£ 

'IJ  - “ij  * “IJ 


(3UM 


90.  In  order  to  examine  the  nature  of  the  restriction  placed  on 
the  response  coeiTicients  by  Equation  309a,  we  substitute  Equation  3^+3 
in  Equation  309a  as  follows: 


(345) 


Therefore,  the  functional  forms  of  and  j must  satisfy  the 
above  differential  equation.  In  the  case  of  quasi-1 inear  elastic  mate- 
rial based  on  Cauchy's  method,  it  was  noted  that  the  functional  forms  of 
and  Ig  were  not  restricted.  If  we  consider  a material  for 

which  depends  on  I^  only  (i.e.,  volumetric  strain  is  caused  only 
by  hydrostatic  stress),  then  Eqiiation  345  reduces  to 


(346) 


Eqtxation  3^6  can  be  satisfied  only  if 


is  of  the  form 


(31*7) 


Equation  3^7  indicates  that  for  a quasi-linear  elastic  material  for 
which  depends  only  on  , the  shearing  response  of  the  material 
is  independent  of  hydrostatic  stress.  In  view  of  Equations'  3^1,  3^3 » 
and  3^7 j the  constitutive  equation  for  this  class  of  material  becomes 


°ij  (^ij  “ 3 Vij) 


(3W) 


The  functions  f^  and  can  be  postulated,  based  on  experimental 

evidence,  for  a given  material.  For  example.  Equation  269  is  a special 
case  of  Equation  3^S  where  the  functions  f^  and  fg  are  obtained  from 
Equations  268a  and  268b,  respectively. 

91.  Equation  3U8  is  the  simplest  form  of  equation  for  quasi- 
linear  elastic  material  in  that  there  is  no  coupling  between  the  devia- 
toric  and  the  volumetric  responses  of  the  material.  Equations  for  more 
complicated  forsis  of  quasi-linear  elastic  materials  can  be  developed  by 
expressing  the  strain  energy  density  function  as  a polynomial  func- 

tion of  and  Ip  , or  by  postulating  mathematical  expressions  for 

I^  and  J^/3  that  will  satisfy  Equation  3^5  and  will  include 
various  degrees  of  coupling  as  desired. 


PART  V;  INCRi2ffilWAL  CONSTITUTIVE  EQUATIONS 


92.  Incremental  constitutive  equations  are  often  used  to  describe 
the  stress-strain  behavior  of  materials  in  vhich  the  state  of  stress  is 
a function  of  the  current  state  of  strain  as  well  as  of  the  stress  path 
followed  to  reach  that  state.  The  general  form  of  the  constitutive 
equation  for  this  class  of  material  behavior  is  generally  expressed  as 


6 . = F , . ( e ,0  ) 

ij  ij  mn  * rs 


(3U9) 


where  F.  . is  a response  function.  Equation  3^9 > which  is  a special 
case  of  Equat’on  19^,  expresses  the  components  of  one  tensor  in  terms  of 
thfj  components  of  two  other  tensors.  Therefore,  in  view  of  Equation  III+, 
the  functional  form  of  F^j  takes  the  following  representation 


6 ik  kp  pj  ik  kp  py 


i I 

I ; 

i 

i 

s.-  f 


+ , 0 , + 0,,  0,  £ , ) 

7 ik  kp  pj  ik  kp  py 


+ 0 4.0a.,  a.,0.  i 

8 ik  kp  pt  tj  xk  xp  pt  tj 


(350) 


where  the  response  coefficients  , . . . , t)0  are  polynomial  func- 
tions of  the  invariants  of  and  and  the  following  joint 

mn  rs  -a  v 

invariants 


n.  = 6 .0. 

1 ab  ba 


(351a) 


II-  = e .0.  0 
2 ab  oc  ca 


(351b) 
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(351c) 


n.  = £ . e.  a 
3 ab  be  ca 


4 6vb  be  cd  da 


(351d) 


93.  We  can  simplify  Equation  350  by  assuming  that  the  materials 
of  interest  are  rate  independent.  To  eliminate  time  effects,  Equa- 
tion 350  must  become  homogeneous  in  time.  This  can  be  accomplished  by 
eliminating  all  terms  containing  second  and  higher  powers  of  in 

Equation  350.  Accordingly,  the  response  coefficients  ri0  , ng  , and 
rig  must  vanish,  , and  must  be  independent  of  and 

functions  of  s+-res8  invariant  alone,  and  j and  Hq  must  be 

of  degree  one  in  e . Imposing  the  above  restrictions  on  the  response 
mn 

coefficients  in  Equation  350,  we  obtain 


*ij  ““  ^O^ij  * '^3'^iJ  '^5^^ik'^kJ  '^ik^kj^ 


+ ri.,(e.iC,  a . + a.,o,  t.)  (352) 

7 ik  kp  pj  ik  kp  pj 


The  response  coefficients  r\^  , 


“ ®0^1 


’^3  ” ^3^1 


% = Hh 


, and  Hq  can  now  be  written  as 


(353a) 

* ®l.“l  * ®5"2 

(353b) 

+ B^n^  + BjHj 

(353c) 

where,  similar  to  , and  t\^  , the  response  coefficients 

a .....  Bq  are  independent  of  £ ^ and  functions  of  stress  invari- 
ants  alone.  Substituting  Equation  353  in  Equation  352.  we  obtain  the 
following  incremental  constitutive  equation  for  rate-independent 
materials 


no 


'i 


■ 

? 

1 

! 


= (Boll  ♦ SiHi  * e^niJSy  + ni^ij 


+ (Bjil  + S^Hi  t 


+ (Bgii  ♦ B^Hi  + 68'’2''’ik“: 


kj 


* "s'^iAj  * “ik^kj’ 


+ n,(e.,a,  a . + a.,  o,  e .) 
7'  ik  kp  pj  ik  kp  pj" 


(35l^) 


Equation  35^  is  the  most  general  incremental  constitutive  relationship 

for  rate-independent  material.  It  contains  twelve  response  coefficients 

which  axe  polynomial  functions  of  stress  invariants.  Since  each  term  in 

Equation  35^  contains  a time  derivative  d/dt  (i.e..  Equation  35^  is 
« 

homogeneous  in  time),  both  sides  of  the  equation  can  he  multiplied  hy 
dt  , resulting  in  the  following  differential  form 


<®0  * h 

* 

* * ®7 


ah 

de 

de 


%a  ^ 


V 0^.  + & 

ah  hfi 


,0,  B 
ah  ha 


de 

5 

8 


ah 

de 

de 


0.  0 )6. , 

he  ca  ij 

ah  °hc%B^°iJ 

ah  %c'’ca^^ik°kj 


. "l  lEjj  ♦ Ofcj  ♦ »ik  d'y) 

♦ "7  VpJ  * ”iAp 


(355) 


^ere  de.,  and  do.,  are  referred  to  as  the  strain  increment  and 
IJ  ^0 

stress  increment  tensors,  respectively.  From  Equation  355  it  is  ap- 
parent that  incremental  constitutive  equations  are  first-ordei*  differen- 


tial equations.  To  obtain  unique  solutions  to  these  equations  ve  need 
to  prescribe  some  initial  conditions.  The  integration  of  the  differen- 
tial equations  for  a given  stress  path  and  initial  condition  will  lead 
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to  stress-strain  relationships.  Varioxis  classes  of  incremental  con- 
stitutive equations  can  be  obtained  from  Equation  355  by  specifying  the 
highest  degree  of  stress  that  appears  on  the  right  side  of  the  equation. 
In  the  following  sections  we  will  develop  and  examine  various  classes 
of  incremental  constitutive  equations. 

Incremented  Constitutive  Equation  of  Grade-Zero 

9^.  If  the  right  side  of  Equation  355  is  independent  of  stress, 
the  incremental  constitutive  relationship  is  referred  to  as  grade-zero. 
In  this  case  = Bg  = 6^  ~ “ ^8  “ ° ’ 

and  and  6q  become  constants.  Thus,  for  grade-zero.  Equation  355 
reduces  to 


«0  * "l 


(356) 


It  is  noted  that  Equation  356  has  the  same  form  as  the  constitutive 
equation  of  linear  elastic  materied  (Equation  201).  In  order  to  in- 
clude the  lineetr  elastic  stress-strain  law  as  a special  case  of  Equa- 
tion 355,  the  material  constants  Bf,  (0,  0,  O)  and  n,  (O,  0,  O)  will 
be  replaced  by  K - ^ 2G  , respectively.  Accordingly,  Equa- 

tion 356  may  be  expressed  as 

toy  . (k  - I o)ae^  * 20  dCjj  (357) 


Equation  357  is  the  constitutive  equation  of  linear  elastic  material 

(see  Equation  210)  expressed  in  incremental  form.  In  order  to  obtain 

a relation  between  stress  and  strain.  Equation  357  must  be  integrated. 

If  we  adopt  the  condition  that  o^j  = 0 when  = 0 , integration  of 

Equation  357  results  in  the  sane  stress-strain  relations  as  predicted 

by  linear  (Elastic  constitutive  equation.  For  example,  consider  the 

condition  of  uniaxial  strain  (Equation  221)  where  dc  « dc,  . For 

nn  1 

this  state  of  deformation  Eqxiation  357  results  in 
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=(K-fo; 

^ de^  2G  de^ 

(358a) 

dOg  = (k  - 

I a)  dc. 

(358b) 

0^  = Og  = 0 when 

(359a) 
(359t) 

which  are  identical  with  the  relationships  predicted  hy  linear  elastic 
stress-strain  law  (Eqtiation  222),  Equation  357*  which  is  the  simplest 
form  of  incremental  constitutive  relationship,  therefore  does  not  mani- 
fest any  nonlinear  behavior. 

Incremental  Constitutive  Equation  of  Grade-One 

95.  If  terms  up  to  the  first  power  of  stress  are  retained  in  the 
right  side  of  Equation  355,  we  obtain  the  incremental  constitutive  equa- 
tion of  grade-one.  This  can  be  achieved  by  allowing  the  response  co- 
efficients , 00  , 6.^  , 00  , , 0jj  , and  0g  to  vanish,  and 

the  ronaining  response  coefficients  to  take  the  following  forms 


6o  = K - f 0 . 

(360a) 

8i“*2 

(360b) 

(360c) 

= 20  ♦ 

(360d) 

"5  • S 

(360e) 

Integrating  Equation  358  and  using  t£.c  condition  that 
= 0 , we  obtain 


U3 


where  through  , K , and  G are  material  constants  that  must 
he  determined  experimentally.  In  view  of  Eqxiation  36O,  the  incremental 
constitutive  equation  of  grade-one  becomes 


^°ij  = - 2 G + X^J^^  de^j^  + Xg  de^^  %a®ij 

^3  ^^nn^ij  Vl)  ^^ij  '^kj  °ik  ^^kj^ 


Equation  36l  contains  seven  material  constants.  It  is  noted  that  Equa- 
tion 361  reduces  to  incremental  constitutive  equation  of  grade-zero 
(Equation  357)  when  the  material  constants  X^  , Xg  , X^  , X^  , and 
X^  vanish.  For  a given  initial  condition  and  stress  path.  Equation  361 
can,  in  principle,  he  integrated  in  order  to  obtain  stress-strain  rela- 
tionships. The  differential  equations  generated  from  Equation  361  for 
var'x''U8  states  of  stress  and  deformation  are  coupled  first-order  equa- 
tions. It  is  not  possible,  in  general,  to  obtain  closed-form  solutions 
for  these  equations.  Therefore,  numerical  integration  schemes  must  be 
utilized. 

96.  Let  us  examine  the  behavior  of  Equation  361  under  hydrostatic 
state  of  stress  (Figure  6a).  For  hydrostatic  state  of  stress.  Equa- 
tion 361  results  in 

= dJ^  = 3(K-|G  + 4 XgJ^  dlj^ 

♦ X^^i  ♦ XjjJj^)  dl^  + I X^J^  dlj^  (362) 

Equation  362  can  be  written  as 

dJ, 

dll  =>  F r-z : : 2 - T '! 

Integrating  Equation  363,  we  obtain  the  following  relationship  between 
volunetric  strain  and  pressure 


llU 


It  is  of  interest  to  note  that  Equation  365  has  the  same  form  as  Equa- 
tion 268a,  which  was  postulated  to  describe  the  pressure-volumetric 
strain  behavior  of  quasi-linear  elastic  materials.  In  fact,  it  is  noted 
that  the  ratio  K/K^  in  Equation  365  conresponds  to  Pg  (the  maximum 
tensile  stress  that  the  material  can  surtain  before  it  breaks)  in 
Equation  268a.  It  is  important,  however,  to  realize  that  unlike  Equa- 
tion 268a,  Equation  365  is  the  outcome  of  a theory  (i.e.,  the  incre- 
mental constitutive  equation  of  grade-one). 

97.  In  order  to  examine  the  coupling  of  the  deviatoric  and 
volumetric  responses  of  the  incremental  constitutive  equation  of  grade- 
one  we  will  examine  the  behavior  of  Equation  36l  for  a state  of  simple 
shearing  defoxmation.  The  strain  increment  tensor  associated  with 
sie^e  shearing  deformation  is  given  as 


where  it  follows  that  de  = 0 . For  this  state  of  deformation,  Equa- 

nn 

tion  361  results  in  the  following  expressions  for  the  nonvanishing 
components  of  stress  increment  tensor 

da^^  = dOgg  = 2(X2  + X^)  de^2  0^2  (367a) 

da^^  = 2X2  *^^12  °12  (36Tb) 

^°12  ■ " °22  * ®33^  S^°ll  *^22^]  *^^12  (367c) 

Equation  367  is  a set  of  first-order  differentieil  equations  which  must 
be  integrated  to  relate  stresses  to  shearing  strain  e^2  • Without 
going  through  the  process  of  integration,  however,  we  csin  draw  certain 
conclusions  about  the  response  of  the  material  in  simple  shearing  de- 
formation, First,  to  maintain  a simple  shearing  deformation  (no  volxane 
change),  normal  stresses  must  be  applied  to  the  boundaries  of  the  speci- 
men. Second,  since  two  of  the  normal  stresses  are  unequal.  Equation  367 
predicts  the  occxurrence  of  normal  deviatoric  stresses  on  the  shearing 
planes.  It  is  recalled  that  the  same  behavior  was  predicted  for  non- 
linear elastic  material  (Equation  235).  As  was  pointed  out  previously, 
the  c-.currence  of  normal  deviatoric  stresses  on  the  shearing  planes  is  a 
secoad-oraer  effect  due  to  tensorial  nonlinearity  (in  this  case  the  last 
term  in  Equation  361)  and  Is  a deparf.ui'e  from  linear  theories  where 
such  nonlinearity  does  not  exist. 

Incremental  Constitutive  Equation  of  Grade-Two 

98.  If  we  retain  terms  up  to  the  second  power  of  stress  in  the 
right  side  of  Equation  355,  the  resulting  incremental  constitutive  equa- 
tion is  referred  to  as  grade- two.  Tims,  for  incremental  constitutive 
equation  of  grade-two  the  response  coefficients  6g  , 8.^  , and  8^ 

vanish,  and  the  roaaining  response  coefficients  in  Equation  355  take  the 
following  forms 

8q  = K - I G ♦ + XgJ^  ♦ X^Jg  (368a) 
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' WiA*  i,M,i  *»  w 
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&!  = Xg  + XgJ^ 

(368b) 

S 

(368c) 

(368d) 

h = '^11 

(368e) 

= ^2 

(368f) 

= 2G  + + X^^Ji  + 

(368g) 

n5  = X5  + x^^j^ 

(368h) 

\ “ ^16 

{368i) 

where  through  X^g  , K , and  G are  material  constants.  In  view 
of  Equation  369,  the  incremental  constitutive  equation  of  grade-two 
takes  the  following  representation 


° ("  ■ I ° * * '6''l  * Va  ) '^'nn  *ij 

♦ (»2  * Vl>  “katj  * ‘9  V’ca'lJ 

. (‘3  * HflJi)  “ij  * ‘u  ^'ab  VlJ 

♦ (20  ♦ . X.jJ^  ♦ XjljJj)  dCjj 

♦ (‘5  > ‘is'flX'i'u  ’kj  ’ “Ut 

♦ *'pj> 


(3(9) 
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Equation  369  contains  l8  material  constants.  It  is  noted  that  Equa- 
tion 369  reduces  to  Equation  361  (constitutive  equation  of  grade-one) 
when  Xg  throxigh  X^g  are  set  to  zero.  Similarly,  we  can  develop  in- 
cremental constitutive  equations  of  grade- three  (or  higher  grades)  by 
retaining  the  third  power  (or  higher  powers)  of  stress  in  the  right  side 
of  Equation  355.  However,  due  to  mathematical  complexities  of  incre- 
mental constitutive  equations  of  grade-one  or  higher,  they  are  seldom 
utilized  to  solve  actual,  engineering  problems.  For  this  reason,  a class 

of  incremental  constitutive  equations,  usually  referred  to  as  variable- 

T 8 

modx’li  constitutive  models,  * is  often  used  for  solution  of  many  en- 
gineering problems.  The  veuriable-moduli  models  are  relatively  simple 
in  that  they  do  not  contain  second-order  effects  due  to  tensorial 
nonlinearity. 


Variable-Moduli  Constitutive  Models 

99.  The  basic  constitutive  relation  of  the  variable-moduli  models 
is  given  by 

“ij  * *lj)  <3^0) 

It  is  noted  that  Equation  370  has  the  same  form  as  the  incremental  con- 
stitutive equation  of  grade- zero  (Equation  357).  In  the  case  of 
variable-moduli  models,  however,  the  equivalent  bulk  and  shear  moduli,' 

K and  G , respectively,  ere  assumed  to  be  functions  of  stress  invari- 

% 

ants.  Depending  on  the  functioned  foxms  of  these  moduli,  various 
classes  of  variable-moduli  models  can  be  constructed.  Since  Equa- 
tion  370  does  not  include  any  second-order  term,  or  terms  involving 
Joint  invariants  of  stress  and  strain  increment  tensors,  it  can  readily 
be  integrated  (for  a given  stress  path  and  initial  condition)  to  yield 
stress-strain  relationships.  Various  classes  of  variable-moduli  models 
are  examined  in  the  following  sections. 

Constant-shear-modulus  model 

100.  As  implied,  for  a constant-she-xr-aodulxis  model,  G is 
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constant  and  K is  a function  of  stress  invariants.  If  we  further 
assume  that  volumetric  strains  are  caused  only  hy  changes  in  pressure 
(i.e.,  there  is  no  coupling  between  volumetric  strain  and  deviatoric 
stresses)*  the  bulk  modulus  K becomes  a function  of  only.  As 

a first-order  approximation*  let  us  assume  that  K is  linearly  related 
to  pressure*  i.e.. 


K = Kq  + K,  ^ 


(3T1) 


where  Kq  (an  initial  bulk  modulus)  and  ^ are  material  constants. 
Substituting  Equation  371  in  Equation  370,  we  obtain  the  following  in- 


cremental constitutive  equation  for  a constant-G  model 


de  6, , + 2G 
nn  ij 


- 1 


de  6,-1 
nn  ij  / 


(372) 


Equation  372  is  also  a special  case  of  incremental  constitutive  equation 
of  grade-one  (Equation  36l).  That  is,  if  we  set  the  material  constants 
^5  * zero,  and  replace  K and  with  Kq 

and  K^/'3  , respectively.  Equation  36l  reduces  to  Equation  372. 

101.  We  will  now  proceed  to  examine  the  behavior  of  Equation  372 
under  various  states  of  stress  and  deformation.  Under  hydrostatic  state 
of  stress.  Equation  372  reduces  to 


(373) 


3Kq  > 


integration  of  Equation  373  (with  the  initial  conditions  that  = 0 
when  = 0)  results  in  the  following  relationship  between  pressure  and 
volumetric  strain 


(37l*) 
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Equation  37*+  is  of  the  same  form  as  Equation  365  (pressure-volumetric 
strain  relation  for  incremental  constitutive  equation  of  grade-one). 

102.  For  simple  shearing  deformation  (Eqxiation  366),  Equation  372 
predicts  the  same  behavior  as  for  linear  elastic  material  (Equation  203). 
Under  conditions  of  uniaxial  strain.  Equation  372  res\ilts  in  the  follow- 
ing expressions 


dOj  = dOg 


i02=(Ko^K,i)dCi-fGde^ 


(375) 


(376) 


Since  there  is  no  coupling  between  volumetric  strain  and  deviatoric 
stresses,  we  can  eliminate  J^/3  from  Equations  375  and  376  by  using 
Equation  37*+.  Substituting  Equation  37*+  in  Equation  375  and  noting  that 


for  •uniaxial  strain  condition,  we  obtain 

^l"l 

do^  = KqB  de^  + I G de^ 


(377) 


Equation  377  can  now  be  integrated  to  relate  vertical  stress  to 
vertical  strain  . Using  the  initial  conditions  that  = 0 when 
0^  = 0 , integration  of  Equation  377  results  in 


(378) 


Similarly,  we  can  obtain  an  expression  for  the  radial  stress 


(379) 


Combining  Equations  379 » 378,  and  37*+,  we  obtain  the  following  relation 
for  the  stress  path  associated  with  the  condition  of  uniaxial  strain 
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(380) 
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Equations  378  through  380  indicate  that  a constant-shear-modulus  model, 
with  a bulk  modulus  which  is  linearly  dependent  on  pressure  (Equa- 
tion 371 ),  can  predict  nonlinear  stress-strain  relationships  under  uni- 
axial strain  condition. 

103.  We  will  next  examine  the  behavior  of  Equation  372  under 
conditions  of  uniaxial  stress.  For  unieixial  state  of  stress,  =. 
and  Equation  372  results  in  the  following  expression  for  the  increment 
of  axial  stress  da^^ 


= (''o  * ^ r)  "'m.  * 2=  - I ■i^nn) 

The  increment  of  volumetric  strain  de  can  be  eliminated  from  Equa- 

nn 

tion  381  by  using  Equation  373  (this  can  be  done  because  volumetric 
strain  is  a function  of  pressure  alone);  thus. 


i 

i 


+ 2G  de^  - 


2Gda^ 

3(3Kq  + iL^o^) 


(382) 


Integrating  Equation  382  (with  initial  conditions  that 
= 0)  , we  obtain 


0^=0  when 


(383) 


Equation  383  relates  axial  strain  to  axial  stress  . Similarly, 
we  can  obtain  an  expression  for  radied  strain  Eg  = in  terms  of 
axial  stress 


(364) 
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Again  it  is  noted  that  highly  nonlinear  stress-strain  curves  can  result 
from  a constant-shear-modulus  constitutive  model. 

lOlt.  Since  the  shear  modulus  is  constant.  Equation  372  can  be 
integrated  to  yield  a total  stress-strain  relation  similar  to  that  for 
quasi-lineax  elastic  models.  In  terms  of  the  stress  and  strain  devia- 
tion tensors,  Equation  372  results  in  (assuming  zero  initia3.  conditions 
or  both  stresses  and  strains) 


(385) 


In  view  of  the  definitions  of  and  E^^  (Equations  122  and  l6l, 
respectively).  Equation  385  can  be  written  in  terms  of  the  stress  and 
strain  tensors 


(386) 


It  is  noted  that  Equation  386  is  a special  caise  of  quasi-linear  elastic 
material  (Equation  26U).  Substituting  Equation  37^**  for  J^/3  > into 
Equation  386,  we  obtain  the  following  quasi-linear  stress-strain  rela- 
tionship for  the  constant-shear-modulus  model 


It  is  noted  that  Equation  307  satisfies  Equation  31*6,  indicating  the 
existence  of  a strain  energy  function  for  the  constant-shear-modulus 
model.  Using  Equations  3l*3  and  387,  we  can  derive  the  following  expres- 
sion for  the  strain  energy  function 


It  can  easily  be  verified  that  if  we  substitute  Equation  388  into  Equa- 
tion 31*1  we  obtain  Equation  367. 
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Constant-Polsi8on*s-ratio  model 

105.  Another  type  of  variable  moduli  model  is  the  constant- 
Poisson* 8-ratio  model  where,  by  analogy  to  linear  elastic  materials,  it 
is  assumed  that  the  ratio  G/K  is  constant.  In  terms  of  elastic 
Poisson’s  ratio  v , this  ratio  is  given  by  (see  Equation  220) 

'539) 

From  Equation  389  it  is  obvious  that  K and  G will  have  similar  func- 
tional forms.  In  order  to  examine  the  consequence  of  this  condition, 
let  us  assume  that  the  functional  form  of  K is  given  by  Equation  371. 
As  was  pointed  out  previously.  Equation  371  indicates  that  volumetric 
strains  are  caused  only  by  changes  in  pressure.  In  view  of  Equa- 
tions 389  and  371  the  constitutive  equation  ttikes  the  following  form 

■‘-’IJ  “ ("0  ♦ *1  r)  "'nr.  hj 

It  can  also  be  shown  that  Equation  390  is  a special  case  of  incremental 
constitutive  equation  of  grade-one  (Equation  361). 

106.  The  behavior  of  Equation  390  under  hydrostatic  state  of 

stress  is  identical  with  that  of  Equation  372.  However,  the  behavior 
of  Equation  390  under  deviatoric  state  of  stress  is  quite  different 
than  that  of  Equation  372-  For  example,  consider  a series  of  constant- 
pressure  shear  tests  (see  Figure  6c)  conducted  at  P^  , Pg  , and  P^  , 
where  ^3  ^ ^ deviatoric  response  of  Equation  390  for  these 

constant-pressure  tests  is  depicted  in  Figure  lla.  It  is  observed  from 
Figure  lie  that  the  deviatoric  response  is  dependent  on  the  super- 
imposed hydrostatic  state  of  stress  (in  this  case  G veo'ies  linearly 
with  pressure).  In  the  case  of  a constant-shear-modulus  model  (Equa- 
tion 372),  on  the  other  hand,  the  deviatoric  response  does  not  depend 
on  the  superimposed  hydrostatic  stress.  Let  us  now  consider  two 
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Figure  11.  Behavior  of  constant-Poisson's-retio  variahle-ooduli 
model  during  constant-pressure  shear  test 

combined  ccnstant-pressure  shear  and  hydrostatic  tests  defined  by  stress 
paths  abed  and  abe'd'  shown  in  Figure  11b.  The  deviator ic  responses  of 
Equation  390  for  these  two  tests  are  also  depicted  in  Figure  Ub.  For 
the  stress  path  abed  it  is  observed  that  the  stress-strain  curve  during 
unloading  (line  cd)  is  above  the  loading  curve  (line  ab).  This  type  of 
behavior  results  in  an  energy^generating  loop  which  is  contrary  to  the 
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observed  behavior  of  real  material.  For  the  stress  path  abc'd*  the 
stress-strain  relation  during  unloading  is  along  path  c*d*,  which  re- 
sults in  permanent  deformation.  Furthermore,  if  the  stress  cycle 
abc'd' a is  repeated  several  times,  the  stress-strain  response  will 
restilt  in  an  unrealistically  excessive  amount  of  deformation  as  shown 
by  the  dashed  lines  in  Figure  11b.  In  order  to  avoid  these  -undesirable 
behaviors,  the  use  of  -variable-moduli  models  for  which  G is  a function 
of  hydrostatic  stress  shotild  be  restricted  to  stress  paths  where 
and  Jg  ranain  constejit  or  increase.  Variable-moduli  models  have  been 

used  to  simulate  the  hysteretic  behavior  of  earth  materials  during 

0 

cyclic  or  near-cyclic  conditions.  For  this  type  of  problem,  two  sets 
of  expressions  are  usually  specified  for  the  moduli  K and  G ; one 
set  for  loading  and  one  set  for  unloading.  A set  of  criteria  or  logics 
are  also  specified  to  determine  whether  the  material  under  considera- 
tion is  being  loaded  or  \inloaded  so  that  the  proper  set  of  mod-uli  can 
be  used.  Application  of  variable-moduli  models  for  treating  hysteretic 
effects  -will  be  discussed  later. 

107.  Let  us  next  consider  the  behavior  of  a constant-Poisson's- 
ratio  model  under  conditions  of  uniaxial  strain.  For  this  state  of 
deformation.  Equation  390  results  in 


. (1  . I B ) (ko  * 


(391) 


(392) 


Ve  can  eliminate  J^/3  from  Equations  391  and  392  by  using  Equa- 
tion 37^.  Substituting  Equation  yjU  in  Equations  391  and  392  and  noting 
that  ^ condition  of  uniaxial  strain,  ve  obtain 


(393) 


d02=(i-f  ilv 


dc. 


(39U) 


Integrating  Equations  393  and  39^  and  using  the  initial  conditions  that 
0^  = Og  = 0 when  = 0 results  in 


(395) 


(396) 


Equations  395 » 396,  and  37^  can  he  combined  to  obtain  an  expression  for 
the  stress  path  associated  with  the  condition  of  uniaxial  strain, 

i.e.. 


As  anticipated,  the  uniaxial  strain  stress  path  for  a constant-Poisson's- 
ratio  model  is  linear  (see  Equations  226  and  227)*  The  constant-shear- 
modulus  and  constant'-Poisson's-ratio  models  are  elementary  versions  of 
variable-moduli  models.  More  complicated  forms  of  variable-moduli 
models  are  discussed  below.  These  models  are  referred  to  as  nonlinear 
variable-mod\ili  models  in  that  both  the  shearing  and  the  volumetric 
responses  are  nonlinear  and  are  represented  by  different  functional 
forms. 

Nonlinear  variable-moduli  models 

108.  More  complicated,  and  perhaps  physically  more  realistic, 
forms  of  variable-moduli  models  can  be  formulated  by  expressing  K and 
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G as  separate  polynomial  functions  of  stress  invariants.  If  we  maiyA 
the  assumption  that  there  is  no  coupling  between  volumetric  strain 
deviatoric  stresses,  then  bulk  modxilus  K becomes  a function  of  J^/3 
alone.  For  a first-order  approximation  we  will  adopt  Equation  371  for 
the  bulk  modulus.  A different  functional  relation  csua  be  postulated  for 
the  shear  modulus  G . As  a first-order  approximation,  we  will  assume 
that  G is  linearly  related  to  the  first  and  second  invariants  of 
stress  tensor,  i.e.. 


G 


(398) 


where  Gq  (an  initial  shear  modulus),  G^  , and  Gg  are  material  con- 
stants. In  view  of  Equations  398  and  371  the  constitutive  relationship 
for  this  type  of  vsuriable-moduli  model  becomes 


de 


nn 


6 


ij 


+ 2 


+ G,  + G„J„ 
1 3 2 2 


(399) 


It  is  noted  that  Eqxiation  399  reduces  to  Equation  372  (constitutive 
equation  of  constant-shear-modulus  model)  if  G^  and  G^  are  set  to 
zero.  Also,  it  can  readily  be  shown  that  Equation  399  is  a special  case 
of  incremental  constitutive  equation  of  grade-two  (Equation  369). 

109.  The  behavior  of  Equation  399  under  hydrostatic  state  of 
stress  is  identical  with  the  behavior  of  constant-shear-modulus  and 
constant-Poisson's-ratio  models,  since  the  functional  form  of  bulk 
modulus  is  the  same  for  all  these  models.  For  deviatoric  states  of 
stress,  however,  the  behavior  of  Equation  399  differs  considerably  from 
that  of  constant-shear-oodulus  and  constant-Poisson's-ratio  models. 

For  example,  consider  a state  of  simple  shearing  deformation  defined  by 
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Equation  366.  For  this  state  of  deformation.  Equation  399  results  in 
(assuming  zero  initial  state  of  stress  and  deformation) 


do^^  = dOgg  = ^0^2  = 0 (liOOa) 

do^2  = 2(Gq  + GgO^gj  de^g  (UoOh) 


Integrating  Equation  UoOh  we  obtain 


(1*01) 


It  is  noted  from  Equations  372  and  390  that  for  this  state  of  deformation 
the  constant-shear-modulus  and  constant-Poisson’s-ratio  models  predict  a 
linear  relationship  between  the  shearing  strain  e^^g  and  the  sheeiring 
stress  o^g  . It  should  also  be  noted  that  within  the  framework  of 
variable-moduli  models  normal  stresses  are  not  required  in  order  to 
maintain  a state  of  simple  shearing  deformation.  As  was  shown  pre- 
viously, to  maintain  a simple  shearing  deformation  in  the  case  of  in- 
cremental constitutive  equations  of  grade-one  (or  higher  grades),  normal 
stresses  must  be  applied  to  the  boundaries  of  the  specimen  (see  Equa- 
tion 367). 

110.  A more  useful  description  for  shear  modulus,  especially  for 

g 

modeling  the  stress-strain  behavior  of  soil,  is  to  express  G in  terms 
of  Jj^/3  and^^  rather  than  Jg  . For  sxample,  for  a first-order 
approximation  we  can  express  G as 

0 •=  Gq  ♦ 0^  ^ ♦ Gg^J^  (U02) 

where  Gg  is  a material  constant.  The  sign  of  Gg  will  determine 
whether  the  material  softens  (Gg  < 0)  or  stiffens  (6g  > O)  during  shear. 
For  this  description  of  shear  modulus  the  constitutive  equation  becomes 
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)("^3  - k *'m.  ^j) 


(1^03) 


For  a state  of  simple  shearing  deformation  (Equation  366),  Equation  1*03 
results  in  (assuming  zero  initial  state  of  stress  and  deformation) 


do^^  = dOgg  = do^^  = 0 


do. 


-12  ' =(=0  * Vl2)^'l2 
Integrating  Equation  UoUb  we  obtain 


(U0l*a) 

(UOlih) 


_ 1 1 I °o  ■*  Via 

‘l2  “ 2 J *“ 

c 


(1*05) 


Eqxiation  U05  also  describes  a nonlinear  relation  between  the  shearing 
strain  e^g  ^1*^  shearing  stress  o^g  . 

111.  Let  us  next  consider  the  behavior  of  Equation  1*03  under 
cylindrical  state  of  strain  utilizing  the  stress  path  depicted  in 
Figure  6d.  For  this  stress  path  vn  = (J^/3  ~ ° * vhere  o^  = Og 

is  the  confining  stress  which  remains  constant  during  the  deformation 
process.  From  Equation  U03  the  increment  of  stress  difference  associ- 
ated with  cylindrical  state  of  strain  becomes 


dOj  - dOj  = 


(itoei 


Ve  can  eliminate  J^/3  fror.  luation  1*06  by  using  the  equation  for  the 
stress  path,  and  noting  that  VI  = (Oj^  - o^)//3  , we  obtain 
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dei  - deg  = 


da^  - dCg 


(1*07) 


For  a given  value  of  confining  siress  Og  » Equation  1*07  can  be  inte- 
grated to  yield  a relation  "between  strain*  difference  e^  - Sg  and 
stress  difference  <^3^  “ Og  * confining  stress  Og  by  0^ 

euid  carrying  on  the  integration,  we  obtain 


Cl  - Cg  = 


f¥ 
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°o-Vc-<°i-VVr-i 

=0  * Ve 


(li08) 


It  is  noted  from  Equation  U08  that  tue  shearing  response  of  the  material 
is  dependent  on  the  magnitude  of  the  confining  stress  0^  . If  is 
set  to  zero  in  Equation  1*08,  this  dependency  disappears  (i.e.,  shear 
modulus  bec<mies  independent  of  hydrostatic  stress,  see  Equation  1*02). 

112.  Next  we  will  consider  the  behavior  of  Equation  1*03  under 
conditions  of  uniaxial  strain.  For  this  state  of  deformation, 
dE^^  = dEj^  and  Equation  1*03  yields  the  following  expression  for  the 


increment  of  vertical  stress  do 


1 * 


*”1  “ ('=0  * ^ r ) 


(1*09) 


Ve  can  eliminate  J^/3  and  from  Equation  1*09  by  using  Equa- 
tion 37I*  and  noting  that  ” Jj^/3)/2  ; thus, 

3 Vl  N,  3%  Ah  ^h-  3 2)[ 
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It  is  noted  that  Equation  UlO  reduces  to  Equation  377  (the  corresponding 
expression  for  a constant-shear-modulus  model)  vhen  Og  and  are 
set  to  yero.  Equation  UlO  is  a first-order  differential  equation  and 
has  the  following  solution  i&ssming  zero  initial  state  of  stress  and 
deformation) 


0 


1 


Again,  it  is  noted  that  Equation  Ull  reduces  to  Equation  378  when  Gg 
and  Gj^  are  set  to  zero.  The  radial  stress  Og  required  to  pi'event 
radial  strain  can  he  determined  hy  direct  integration  of  Equation  U03 
(similar  to  the  procedure  followed  to  obtain  o^)  or  hy  using  Equa- 
tions 37^  and  Ull.  From  Equation  37^*  it  follows  that 


^^2 


(Ul2) 


In  order  to  relate  to  axial  strain  we  sitiply  eliminate 
from  Equation  Ul2  hy  using  Equation  Ull.  The  uniaxial  strain  stress 
path  for  this  model  can  he  determined  hy  using  Equations  Ull  and  37^. 
Equation  37**  can  he  inverted  to  relate  1^^  in  this  case)  to  , 

The  resulting  relation  can  then  he  substituted  into  Equation  **11  to 
express  in  terms  of  J^/3  which,  in  conjunction  with  the 
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expression^  Jg  = (o^  - J^/3)/2  , will  result  in  the  stress  path. 

113.  More  complicated  forms  of  nonlinear  variable-moduli  models 
can  be  developed  and  analyzed  as  for  the  preceding  models.  The  choice 
of  any  particular  model,  however,  must  be  based  on  the  experimental 
observation  of  the  stress-strain  behavion  of  the  material  of  interest. 
Treatment  of  hysteretic  behavior 

111*.  As  was  pointed  out  previously,  variable-moduli  models  have 
been  used  to  simulate  the  hysteretic  behavior  of  earth  materials  during 
cyclic  or  near-cyclic  loading  conditions.  To  show  the  procedure  by 
which  the  hysteretic  behavior  is  simulated,  we  express  the  basic  consti- 
tutive relation  of  the  variable-modxili  models  (Equation  370)  in  terms 
of  the  hydrostatic  and  deviatoric  components,  i.e.. 


dJ 

-T  = ^ 


(l*13a' 


dS^j  = 2G  dE^j  (J4l3b) 

where  dS^^  and  dE^^  are,  respectively,  the  deviatoric  stress  and 
strain  increment  tensors.  It  is  postulated  that  the  basic  form  of 
Equation  U13  is  valid  for  all  loading  conditions  (initial  loading,  un- 
loading, and  reloading).  However,  the  functional  forms  of  K and  G 
change  depending  on  whether  the  material  under  consideration  is  being 
loaded  or  unloaded.  As  depicted  graphically  in  Figure  12,  for  initial 
loadings  the  response  of  the  material  is  governed  by 


K,  dl, 


(itlUa) 


= “t  “ij 


Cl.ll4b) 


where  and  are,  respectively,  the  bulk  and  shear  moduli  etssoci- 
ated  with  iniMal  loading.  During  unloading  and  reloading  we  assume  that 
the  response  of  the  material  is  governed  by 
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dl^  (4l5a) 

dS.  = 2G  dE.,  (4l5b) 

u ij 

where  and  are,  respectively,  the  hulk  and  shear  moduli  as- 
sociated with  unloading  and  reloading. 

115.  To  complete  the  specification  of  the  model  we  need  a cri- 
terion to  define  the  loading  condition  during  a deformation  process. 

We  adopt  as  our  criterion  the  quantity 

aw  = cjj  dcy  (ki6) 


which  defines  the  rate  at  which  the  stresses  do  work  during  the  deforma- 
tion process.  According  to  this  criterion,  dW  > 0 defines  loading 
(initial  loading  or  reloading)  and  dw  < 0 defines  unloading.  The 
condition  dW  = 0 is  referred  to  as  neutral  loading.  The  neutral 
states  of  loading  associated  with  the  rate  of  work  criterion  impose 
c...’tain  restrictions  on  the  material  constants  in  the  constitutive  equa- 
tions for  loading  and  unloading  and  require  special  considerations. 

The  material  constants  must  he  chosen  sc  that  the  loading  and  unloading 
constitutive  equations  become  identiceuL  whenever  dW  = 0 , i.e.,  neutral 
loading.  This  requirement  must  he  met  in  order  to  obtain  a unique 
solution  for  a boundary-value  problem  involving  cyclic  loading  condi- 
tions. From  Equations  HlH  and  itl5  it  is  apparent  that  variable-moduli 
models,  in  general,  do  not  satisfy  this  requirement. 

116.  In  view  of  the  definition  of  deviatoric  stiain  increment 
tensor  we  can  express  Equation  U16  in  the  following  form 


or 


(HlTa) 


(l*17b) 


We  can  eliminate  emd  dl^  from  Equation  UlTb  by  using  Eqiaa- 

tion  ltl3;  thus. 


S dS  J,  dJ^ 
dw  = -22rr-22.  + ^ 


2G  9K 

Since  J*  = 1/2^S  S V it  foUovs  that 
2 \ mn  mn/ 


(1*18) 


dJl  = dS^„ 
2 mn  mn 


(1*19) 


where  dJ^  is  the  incrattent  of  the  second  invariant  of  stress  devia- 
tion tensor.  Equation  1*18  can»  therefore,  be  written  as 


dJi  J,  dJ 


(1*20) 


For  variable-mod’xLi  models  in  which  G = G(J^)  and  K = K(Jj^)  the  rate 
of  work  can  be  separated  into  hydrostatic  and  deviatoric  components. 
Equation  1*20  can  then  be  used  with  the  interpretation  that  dJ^/2G  is 
rate  of  work  due  to  deviatoric  stresses  and  dJ^/9K  is  rate  of  work 
due  to  hydrostatic  stress.  Accordingly,  two  criteria  for  defining 
V€o:ious  loading  conditions  can  be  prescribed.  For  the  deviatoric  part 
of  deformation,  loading  and  unloading  are  defined  according  to  whether 
dJg  is  positive  or  negative,  respectively.  For  the  hydrostatic  part 
of  deformation,  loading  and  unloading  are  defined  according  to  whether 
dJ^  is  positive  or  negative,  respectively.  In  this  manner,  it  is 
possible  for  the  material  to  xinload  in  shear  while  loading  in  pressure 
or  vice  versa.  It  should  be  pointed  out  that  these  criteria  also  do 
not  satisfy  the  requirement  of  neutral  loading.  Because  of  the  require- 
ment of  neutral  loading,  the  validity  of  variable-moduli  models  for 

9 

•creating  hysteretic  effects  has  been  questioned.^  Hysteretic  effects 
and  permanent  deformation  can  be  treated  within  the  framework  of  in- 
cremental theory  of  plasticity  without  violating  the  requirement  of 
neu'cral  loading. 
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PART  VI:  CONSTITUTIVE  EQUATIONS  OF 
SIMPLE  VISCOELASTIC  MATERIALS 


117*  For  viscoelastic  materials,  the  state  of  stress  is  a func- 
tion of  both  the  state  of  strain  and  time  rate  of  strain.  The  stress 
tensor  can,  therefore,  be  expressed  as 


a..  * F. .(e^^  , i) 

ij  ij  mn  ra 


(1^21) 


where  F.  = viscoelastic  response  function.  In  view  of  Equation  ll^i, 
the  response  function  F^j  takes  the  following  form 

■’iJ  = "o^lj  * 

* ^ * ‘tk=kj) 

* *^6  (^ik^kp^pj  * ^ik^kp^pj) 

* "TC^ikVrJ  * "ikVpj) 

* ''8('lkVptHj  * ^ik'kp'pt'tj) 

where  the  response  coefficients  Hq  , . • • , n0  are  polynomial  functions 

of  the  invariants  of  e and  e and  the  following  Joint  invariants 

mn  rs 


"l“ 

£ -uC-u 

ab  ba 

(423a) 

“2  = 

e , e,  e 
ab  be  ca 

(423b) 

"3  = 

e . e,  e 
ab  be  ca 

(423c) 

II 

^ab'ijc^cd^da 

(423d) 

It  is  noted  that  when  dependence  on  c disappears.  Equation  422  re- 
duces  to  the  constitutive  equation  of  Cauchy  elastic  material  (Equa- 
tion 197).  Various  classes  of  viscoelastic  materials  can  be  described 
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by  Equation  422  by  proper  selection  of  the  response  coefficients 

tIq  , . . . j n0  . We  will  be  dealing  with  simple  viscoelastic  materials 

in  this  report. 

Kelvin-Voigt  Material 

118.  The  constitutive  equation  of  Kelvin-Voigt  material,  which 
describes  the  simplest  type  of  viscoelastic  material,  can  be  obteiined 
from  Equation  422  by  allowing  the  coefficients  rig  , » Hg  » 

, and  rig  to  vaiish  and  the  remaining  coefficients  to  take  the  fol- 
lowing forms 


^0  = 

xi^  + X i, 

1 V 1 

(424a) 

2G 

(424b) 

"3' 

2G 

V 

(424c) 

Accordingly,  the  constitutive  equation  of  Kelvin-Voigt  materieil  becomes 


"ij 


V 1 ij 


(425) 


The  coefficients  X^  and  G^,  , analogous  to  the  elastic  constants  X 
and  G , are  the  dilatationed  €ind  shesir  viscosity  coefficients,  respec- 
tively. Again  it  is  noted  that'  if  X^  and  G^  are  set  to  zero  Equa- 
tion 425  reduces  to  the  constitutive  equation  of  linear  elastic  material. 

119.  Let  us  examine  the  behavior  of  Kelvin-Voigt  material  under 
uniaxial  state  of  strain  defined  by  Equation  221.  For  this  state  of 
deformation.  Equation  425  results  in 

= Xe^  ♦ X^Ej^  + 2Ge^  + SG^Cj^  (426a) 

°2  * ^®1  ^v^l  (426b) 

where  is  the  strain  in  the  direction  of  motion  in  a uniaxial  strain 
configuration  and  Og  is  the  lateral  stress  required  to  preven-  lateral 
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strain.  Equation  U26a  can  be  written  in  a more  compact  form  by  collect- 
ing texms 

o^  = (X  + 2G)e^  + (X^  + 2G^)e^  (U27) 

The  te^iB  X + 2G  is  recognized  as  the  constrained  modulus  M (see 
Equation  225).  Analogous  to  M , we  denote  the  term  X^  + 2G^  by  , 
which  is  the  viscosity  coefficient  associated  with  the  conditions  of 
uniaxial  strain.  Equation  1*27  may  now  be  written  as 

o,  = Me.  + M e,  (1*28) 

J.  X V J. 

Equation  1*28  is  the  counterpart  of  the  differential  equation  of  motion 
for  a linear  spring  and  a linear  dashpot  connected  in  paurallel  (Fig- 
ure 13). 

120.  V7e  will  now  exaaiine  the  behairior  of  Equation  1*28  for  an  ap- 
plied constant  stress  of  magnitude  Oq  . For  an  applied  constant  stress 
Oq  , Equation  1*28  becomes 


Cq  = Me^ 


Vl 


(1*29) 


Integrating  Equation  1*29  (with  initial  conditions  that  e^  = 0 at 
t = 0)  we  obtain 


e 


1 


-(M/M^)t 

e 


(1*30) 


Equation  1*30  describes  the  strain-time  response  in  uniaxial  strain  con- 
dition due  to  an  applied  constant  stress  of  magnitude  (Figure  13). 
In  view  of  Equations  1*30  and  l*26b,  the  lateral  stress-time  response  dur- 
ing application  of  becomes 


^2 


X 

M 


♦ 


M “ m)® 
k V / 


(1*31) 


121.  Let  MB  next  examine  the  behavior  of  Equation  1*28  due  to  ap- 
plication of  a time- varying  stress  condition  defined  by 


KELVIN- VOIGT  liATERIAL 


a.  APPUEO  STRESS 


b.  STRAIN-TiyE  RESPONSE 

Figure  13*  Behavior  of  Kelvin-Voigt  material  ia 
uniaxial  strain  configuration 


Oi  “ Oq(1  - e"®^) 


(1*32) 


^ere  c is  a constant.  Substituting  Equation  1*32  in  Equation  1*28  and 
integrating  the  resulting  expression  (with  initial  conditions  that 
* 0 at  t = 0)  we  obtain 


, ,°oL  -"v 

^1  M ^ (M  - cM^)  ® I 


K - cM 


(1*33) 


Equation  ^33  describes  the  strain-tine  response-  in  uniaxial  strain  condi- 
tion due  to  application  of  stress-time  history  given  by  Eqmtion  U32. 

We  can  eliminate  the  time  t between  Equations  ^32  and  U33  to  obtain 
the  following  strain-stress  relationship 


e 


1 


/ 

M(M  - cM.^) 


co-M 
0 V 


{h3h) 


Equation  U3it  indicates  that  the  stress-strain  response  is  not  unique  and 
is  dependent  on  the  constant  c . 

122.  Analogous  to  nonlinear  eleistic  material,  we  can  construct 
nonlinear  viscoelastic  constitutive  relationships  by  retaining  some  of 
the  second-order  terms  in  Equation  422.  For  example,  taking  = Hg 
= = ^8  ® 0 Equation  422,  a second-order  viscoelastic  constitutive 

relationship,  often  referred  to  as  the  nonlinear  Kelvin  solid,  esm  re- 
sult, i.e. , 


“ij  ' VlJ  * ViJ  * Vik\j  * "i^lj  * <'•35  > 

It  is  noted  that  in  Equation  435  there  is  no  tensorial  coupling  between 
the  strain  and  strain-rate  tensors.  Various  classes  of  nonlinear  Kelvin 
solid  can  be  developed,  similar  to  nonlinear  elastic  material,  by  expand- 
ing the  response  coefficients  , . . . , in  terms  of  the  invariants 
of  the  strain  and  strain-rate  tensors. 


Maxwell  Material 


123.  Rate-dependent  constitutive  relationships  can  also  be  ex- 
pressed as 


F,,(e 
ij  mn 


o ) 

rs 


(436) 


Equation  436  is  identical  with  Equation  349*  The  ax>st  general  form  of 
Equation  436  is  given  by  Equation  350  'md  contains  nine  response 
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coefficients.  If  ve  neglect  all  second-  and  higher-order  terms  in  Equa- 
tion 350  (i.e. , let  hg  = = hg  = = hg  = O),  and  assume  that 

the  remaining  coefficients  take  the  following  forms 


^0  = “Vl  * 


X I, 
m 1 


(U3Ta) 


\ = 2Gm 


^3  = -2®m 


(U37b) 

(li3Tc) 


Where  a , X , G , aud  0 are  material  constants,  we  obtain 
m m m m 


®ij  * ^-^m'^1  * Vl^^iJ  * 2G^^iJ  - 


ij 


(438) 


Equation  438  is  the  constitutive  equation  of  Maxwell  material.  It  con- 

tedns  four  material  constants.  It  is  noted  that  if  we  set  a and  0 

m m 

to  zero  and  integrate  Equation  438,  we  obtain  the  constitutive  equation 

of  linear  elastic  material  (replacing  G by  G and  X by  X). 

m m 

124.  Let  us  consider  the  behavior  of  Maxwell  material  under  uni- 
axial state  of  stress  defined  ty  Equation  215.  For  this  state  of  stress. 
Equation  438  results  in 


0 • (-“A  ♦ Vl>  * “m'2 


(li39.) 

(ii39b) 


where  is  axieil  stress  and  and  axial  and  lateral 

strain-rate  components,  respectively.  Since  for  uniaxial  state  of 
stress  and  , Equation  439  can  be  written  in  the 

following  form 


6.  ♦ (a  ♦ 20„)o.  = (X^  20^  )e.  ♦ 

1 m ml  m mi  me 


(440a) 


a 0, 
m 1 


Vl* 


(2X_  t 


m 


(UtOb) 


from  Equation  U5i0a  by  using  Equation  UUOb;  thus. 


We  can  elimii:ate  e. 


^1  = 


\ ^ G 
m m 

^3G  X + 2G 
‘mm  m 


30  X + a G + 20  G 
mm  mm  m m 

3G  X + 2G^ 
mm  m 


Denoting  (3G  X + 2G^)/ (X  + G ) by  E (elastic  Young's  modulus) 
\mm  m/^m  m m 

and  (3G  X + 2Gf  y(30„X  + a G + 20  G ) by  n , Equation  UlU  can  be 
\ m m m//  \ mm  mm  mm/  m 


written  as 


• -fix'll 

^1  “ E n 
m m 


(4U2) 


Equation  UU2  is  the  counterpart  of  the  differential  equation  of  motion 
for  a linear  spring  and  a linear  dashpot  connected  in  series  (Figure  lU). 
If  0^  is  suddenly  applied  and  then  held  constant  d^  = 0 , and  integra- 
tion of  Equation  i^U2  results  in  a steady  linear  increase  of  axial  strain 
with  time  (Figure  lU). 

125.  Analogous  to  nonlinear  Kelvin  solid,  we  may  develop  a non- 
linear Maxwell  solid  by  taking  = ng  = Hy  = bg  = 0 in  Equation  350. 
The  constitutive  equation  of  nonlinear  Maxwell  solid  then  becomes 


"o^ij  * "i^ij  * * ''s’ij  * Vik'n) 


(W3) 


Equation  Uk3  can  be  used  to  construct  various  classes  of  nonlinear 
Maxwell  solid  by  expanding  the  response  coefficients  tIq  , . • • , in 
terms  of  the  invariants  of  the  stress  and  strain-rate  tensors. 


Standard-Linear  Material 


126.  Rate-dependent  constitutive  relationships  can  also  be  ex- 


pressed as 


(UUU) 


Equation  UUU  reduces  to  Equation  U36  if  dependence  on  e disappears. 

i^q 

The  simplest  fora  of  Equation  Uhh,  which  is  referred  to  as  the 


MAXIELL  MATERIAL 


a.  APPLIED  STRESS 


ELASTIC  RESPONSE -o^/E, 


tL  STRAIN'TIIE  RESPONSE 


Figure  lU.  Behavior  of  Maxwell  material  In 
uniaxial  stress  configuration 

constitutive  equation  of  standard-linear  materials*  is  expressed  as 


d . . = /-a  o-  + X I,  ♦ X I.  Jd.,  + 2G  e, , ♦ 2G  e . , - 

ij  \sl  cl  sl/ij  sij  sij 


(Us) 


where  , X^  , X^  , » Gg  » and  Bg  are  material  constants.  As 

expected*  if  G and  X are  set  to  zero*  Equation  takes  the  form 
8 8 

of  the  constitutive  equation  of  Maxwell  material  (Equation  U36). 

127 . Let  us  consider  the  deviator ic  response  of  standard-linear 


material.  From  Equation  it  follows  that 


S.. 

J-j 


2G  E., 
s ij 


+ 2G  E., 
s ij 


Equation  1|1|6  is  the  counterpart  of  the  differential  equation  of  motion 
for  a linear  spring  and  a Kelvin-Voigt  element  connected  in  series  (Fig- 
ure 15 ) . It  can  be  shown  that  the  materiail  constants  G , G , and 

s s 

B and  the  parameters  of  the  corresponding  spring  and  dashpot  model  are 
s 

related  through  the  following  relationships 


2n  q 
s^s 

q + p 


(l*i^7a) 


‘Is^ 


G 

s 


2p  q 
q + p 


(J*U7c) 


If  the  deviatoric  stress  (say  S^g)  is  suddenly  applied  and  then  held 
constant  = 0 , and  integration  of  Equation  results  in  the  fol- 
lowing deviatoric  strain-time  response 


°s  \ 

26  0 j 

s s/ 


(Uh8) 


Equation  1*48  is  depicted  graphically  in  Figure  15.  As  indicated  in  Fig- 
ure 15,  the  material  exhibits  an  initieil  elastic  response  (similar  to 
Maxwell  material) 


s 


(449) 


and  an  asymptotic  elastic  behavior  (similar  to  Kelvin-Voigt  material) 


Vs  ^ Ps’ 


(450) 
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128.  Simileur  to  the  nonlinear  constitutive  equations  of  Kelvin 
and  Maxwell  materials,  we  may  develop  a constitutive  equation  for  a non- 
linear standard  solid  containing  second-order  terms  in  stress,  strain, 
and  strain-rate  tensors.  This  can  he  accomplished  by  combining  Equa- 
tions h35  (nonlinear  Kelvin  solid)  smd  Uh3  (nonlinear  Maxwell  solid), 
i.e. , 


= "o^lj  * VlJ  * ''2'il'kJ  * "3H3 

Various  classes  of  nonlinear  standard  solid  can  be  developed  by  expand- 
ing the  response  coefficients  , . . . , Hg  in  terms  of  the  invariants 
of  the  strain,  stredn-rate,  and  stress  tensors. 

129.  Constitutive  equations  cf  viscoelastic  materials  can  also  be 
formulated  in  integral  forms  or  series  forms  with  differential  oper- 
ators.^® Discussion  of  these  types  of  constitutive  equations  is  beyond 
the  scope  of  this  report. 
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PART  VII:  CONSTITUTIVE  EQUATIONS  OF  PLASTICITY 


S 


i' 

a 


I 


iiiWwiiiwiAi'i?6t»wa:?g! 


130.  Constitutive  equations  of  plasticity  are  designed  to  de- 
scribe the  stress-strain  behavior  of  hysteretic  materials.  The  basic 
assumption  employed  in  developing  these  equations  is  that  for  each  load- 
ing increment  the  corresponding  strain  increment  can  be  considered  as 
being  the  sum  of  the  plastic  (permanent)  and  elastic  (recoverable) 
strains.  Mathematically,  the  strain  increment  tensor  de. , is  ex- 


pressed as 


a 


Q ‘T' 

where  de^j  and  de^j  are,  respectively,  the  elastic  and  plastic 
strain  increment  tensors.  The  elastic  strain  increment  tensor  is  given 
in  tewas  of  incremeuteil  elastic  relation 


de:.=5i^ 


e 

"a 


2G 


dJi 

9r'« 


(1*53) 


The  elastic  moduli  G and  K can  be  assumed  to  be  constant  or  fionc- 
tions  of  str-^88  invarients,  as  dictated  by  test  data.  However,  to  be 
consistent  with  path  dependency  of  elastic  materials,  and  to  eliminate 
any  possibility  of  energy  generation  or  hysteretic  behavior  during 
elastic  deformation  (see  Figure  11 ),  the  forms  of  G and  K should  be 
restricted  to 


G = G(Jp 


(1*5W) 


K = K(J^)  (U5Ub) 

During  unloading,  or  during  loading  where  the  state  of  8ti*ess  is  below 
a specified  state  referred  to  as  yield  stress,  the  behavior  of  the  mate- 
rj.al  is  defined  coB^letely  by  Equation  U33.  At  the  onset  of  yielding 
(when  the  state  of  stress  is  such  that  the  yield  stress  is  reached) , and 
during  sidjsequent  loading,  the  material  will  experience  both  elastic  and 


plustic  deformation  and  Equation  U52  vill  govern  the  behavior  of  the 

material.  Therefore,  for  a complete  description  of  the  material  we 

need  to  specify  the  form  of  the  plastic  strain  increment  tensor  dc^j  . 

Guidelines  as  to  how  the  plastic  strain  increment  tensor  can  be  spec- 

11 

ified  were  established  by  Drucker  by  introducing  the  concept  of  mate- 
rial stability.  For  a stable  material,  the  work  done  by  a set  of  stress 
increments  wheii  applied  on  a specimen  of  the  material  is  positive. 
Furthermore,  if  the  stress  increments  are  I'emoved,  the  net  work  per- 
formed by  them  d ring  the  load-unload  cycle  is  zero  or  positive.  If  we 
denote  the  set  of  stress  increments  by  da^^  , and  denote  the  correspond- 
ing change  in  the  state  of  strain  due  to  application  of  dct^^  by  de^^  , 
then  the  first  condition  of  stability  caji  be  stated  as 

d=ij  dSy  > 0 (>.55) 

In  view  of  Equation  ^52,  the  second  condition  of  stability  can  be  ex- 
pressed as 


'^=ij  ' '^“ij  '^=13  - “ <‘56) 

Equations  1*55  and  1*56  provide  guidelines  for  determining  the  form  of 
the  plastic  strain  increment  tensor.  We  also  need  to  specify  the  form 
of  the  yield  function  ^ , which  defines  the  limit  of  elastic  behavior. 
Depending  on  the  specification  of  ^ , various  types  of  plastic  consti- 
tutive equations  can  be  established.  In  general,  we  will  be  concerned 
with  ideal  and  work-hardening  plastic  materials. 

Ideal  Plastic  Material 


131.  For  an  idesfl  plastic  material  the  yield  f)inction  (J  (or 
yield  surface)  is  fixed  in  the  principal  stress  space,  i.e.,  it  does  not 
move  cr  expand  during  plastic  deformation.  The  yield  sui'face  is  only  a 
function  of  stress  tensor,  or  function  of  invariants  of  stress  tensor 
for  an  isotropic  material.  Unlimited  plastic  flow  takes  place  when 


1U8 


= k , where  k is  a material  constant  defining  the  onset  of 
yielding.  During  plastic  deformation  then 

The  stability  condition  for  an  ideal  plastic  material  is  specified  'ey 


(U58) 


For  uniaxial  stress  configuration,  for  example.  Equation  lt58  indicates 

that  during  plastic  deformation  the  stress  remains  constant  while  the 

litrain  increases.  Since  all  admissible  stress  increments  do. . satisfy- 

* J 

ing  Equation  1*57  must  also  satisfy  the  stability  postulate  given  in 
Equation  1*58,  it  follows  that 


s A -jjL- 


(1*59) 


where  A is  a positive  scalar  factor  of  proportionality  and  is  depen- 
dent on  the  particular  form  of  the  yield  function  ij  . Equation  1*59  is 
often  referred  to  as  the  plastic  flow  rule.  Inherent  in  Equation  1*59 
is  the  normality  condition  which  indicates  that  the  plastic  strain  in- 
crement (viewed  as  vector)  is  normal  to  the  yield  surface  i . Another 
consequence  of  the  stability  postulate  is  the  convexity  condition,  which 
requires  that  the  yield  surface  ^ must  bo  convex  in  the  principal 
stress  space.  In  view  of  Equations  1*53  and  1*59,  the  complete  expression 
for  the  strain  increment  tensor  becomes 


20  9K 


'ij 


3o 


ij 


(1»60/ 


Equation  1*60  prevails  in  the  plastic  range  (dj  = 0).  In  the  elastic 
range,  and  during  unloading  from  a point  on  the  yield  surface  (d^  < O), 
the  behavior  of  the  material  is  governed  by  Equation  1*53. 

132.  In  order  to  use  Equation  i*60  we  must  determine  the  form  of 
the  proportionality  factor  A . This  con  be  accouqtlished  by  combining 


ll*9 


Equations  457  an<i  460.  From  Equation  460  we  can  determine  the  stress 
increment  tensor  da. 


"“ij ' “ 35^  - i) 


(461) 


Substituting  Equation  46l  in  Equation  457  we  obtsdn 


We  can  eliminate  dJ^  from  Equation  462  by  using  Equation  460..  J^um 
Equation  460  it  follows  that 


dJ,  = 3K  (dr  - A ^ 

1 \ 1 30^^^  ij^ 


(463) 


In  view  of  Equations  463  and  462,  the  proportionality  factor  A takes 
the  following  form 


H . . 3K  - 20  ,,  H . 

, 3sr  ‘‘'ij  * -ns—  ^ *ij 

A = ^ 


(464) 


^‘^ij  ®°ij  ■ 


It  is  iiJioted  that  all  indices  in  Equation  464  are  dummy  indices,  indicat*- 
ing  the  scalar  character  of  A . Equations  464  and  46o  cem  now  be  ccmt- 
bined  to  give  an  egression  for  the  strain  increment  tensor 


dSj,  dJ, 

‘■"tj  ° lo  * # *1J  * 
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(465) 


It  follows  from  Equations  46l,  463,  and  464  that  the  stress  increment 
tensor  takes  the  following  representation 


150 


iiiiiii 


T 


'iO 


? :<  1 
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ir~  ^ -Tr~  3^  V 

do  = 2G  dE.,  + K dl  6 - 
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In  order  to  use  Equations  465  and  466  ve  only  need  to  define  the  form 
of  the  yield  function  ^ for  a particular  material  of  interest.  For 
a number  of  engineering  materials,  particularly  soils,  the  yield  func- 


(467) 


tion  is  generally  expressed  in  terms  of  and  VJ^  i-e.. 


l5(aij)  = (((ji  , V^)  = k 


For  the  above  specification  of  ^ it  follows  that 

^ ^ . -ilL.  I!(S  . IL  5 5L-S 
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Application  of  Equation  468  in  Equation?  465  and  466  results  in 
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(469) 


K dl^  6^^ 


3K  dl 


-^=  S dE 

^ Vj^  sV^  ““ 


9K 


(^)*‘(^) 


(h70) 


Equations  469  and  470  are,  therefore,  special  cases  of  Equations  465 
and  466,  respectively,  where  the  yield  function  <J  is  restricted  by 
Equation  46?.  In  the  next  section  we  will  discuss  the  procedure  by 
which  these  equations  can  be  utilized  for  specific  yield  functions. 
Prandtl-Reuss  material 

133.  Prandtl-Reuss  material  is  the  most  widely  used,  and  perhaps 
the  simplest,  ideal  elastic-plastic  material.  The  yield  condition  asso- 
ciated with  the  Prandtl-Reuss  material  is  the  well-known  Von  Mises  cri- 
terion given  by 


(un) 


Equation  471  describes  a right-circular  cylinder  in  the  principal  stress 
space  with  its  central  axis  the  line  of  hydrostatic  stress  as  shown  in 
Figure  l6.  When  the  state  of  stress  is  such  that  Equation  471  is  satis- 
fied, the  material  would  flow  plastically,  undergoing  plastic  as  well  as 
elastic  strains.  When  the  stresses  are  less  than  those  satisfying  Equa- 
tion 471*  the  materieil  will  undergo  elastic  strains  only. 

134.  In  order  to  obtain  the  constitutive  equation  of  Pra  dtl- 
Reuss  materieil  we  simply  substitute  Equation  471,  for  the  yield  function 
^ , into  Equation  469.  Completing  the  substitution,  and  considering  the 
fact  that  during  plastic  deformation  Vj^  ” ^ • we  obtain 


dS, , dJ,  S dE 

de  » — — i.  « ♦ -H£ ML  s 

‘^^ij  2G  ^ 9K  ®ij  * ^2  ^ij 


(472) 


Figiire  l6.  Vcn  Mises  yield  surface  in  principal  stress  space 

Similarly,  substituting  Equation  U71  in  Equation  U70  we  obtain  the  fol- 
lowing e35)ression  for  stress  increment  tensor 

GS  dE 

"“ij  = “ =lj  "•73) 

The  quantity  dE^  in  Equations  U72  and  1*73  is  recognized  as  the 
rate  of  work  due  to  distortion.  Expanding  this  quantity  with  respect  to 
the  plastic  and  elastic  components  we  obtain 

(dE®  + (1*71*) 

mn  mn  mn  \ mn  mn/ 

Since  dE®  = dS  /2G  (see  Equation  Ul3b),  Equation  U7I*  becomes 

mri  J5JJ  • 

“to  = <-  <‘‘75> 


The  quantity  S dS  is  the  increment  of  the  second  invariant  of 

• ^ mn  mn 


stress  deviation  tensor  (see  Equation  4l9)  and  is  zero  for  the  Von  Mises 
yield  criterion.  Equation  reduces  to 


S dE  = S dE^ 
rnri  nui  nin  nin 


(ki6) 


indicating  that  in  the  plastic  range  the  rate  of  distortional  work  is 
only  due  to  plastic  deformation.  Also,  from  Equations  1^72  and  '+73  it 
follows  that 


dJ, 


, _ 1 - e 

^"ii  = W = '^^ii 


(^77) 


In  view  of  Equation  1+52,  Equation  1+77  indicates  that 


dcPj  = 0 


(1*78) 


That  is,  no  plastic  volume  change  can  occur  in  the  plastic  range  for 
Praudtl-Reuss  material. 

135*  We  can  now  summarize  the  Prandtl-Reuss  equation  in  the  fol- 
lowing manner.  During  elastic  loading  < emd  during  unloading 

( (3j(/3o.  Jdo. , < o),  the  elastic  constitutive  equation  (Equation  1+53) 
prevails.  In  the  plastic  range  = k and  (3((/30j^j  )dOj^j  = 0^,  Equa- 

tion 1+72  (or  Equation  1+73)  governs  the  behavior  of  the  material.  The 
Prandtl-Reuss  constitutive  equation  can  then  be  expressed  as 


* K dl^  dy 


when 


V3T  < k or 


. da. , < 0 

SOij  ij 


(l+79a) 
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when  VjI”  = h and  da . . = 0 

2 ij 


(1+791.) 


136.  In  order  to  demonstrate  the  application  of  Equation  1+79,  we 
will  examine  the  behavior  of  Prandtl-Reuss  material  under  conditions  of 


151+ 


uniaxial  strain.  For  uniaxial  strain  conditions  the  strain  increment 
and  strain  deviation  increment  tensors  are  given  as 


de^  0 0 


0 0 0 
0 0 0 


(ItSOa) 
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- 3*^^! 


(480b) 


In  the  elastic  range  the  behavior  of  the  materiel  is  governed  by  Equa- 
tion J*79a. 

dOi  = + I oj  de^  = M de^  » — ) dJ^^  (i*8la) 


2G 

dOj^  - dOg  = do^  - do^  = 2G  de^^  = dJ^ 


(48lb) 


For  virgin  loading  in  the  elastic  reuage.  Equation  l*8l  governs  the  be- 
havior of  the  material.  It  should  be  noted  that  if  the  initial  state  of 
stress  and  surain  is  zero,  for  virgin  loading,  Equation  U8l  can  be  used 
in  terms  of  total  rather  than  incremental  quantities.  In  the  uniaxial 
strain  test 

VF=^(a  - a)  (li82) 

V3 


Thxid  the  materieLL  will  yield  when 
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In  viev  of  Equations  and  U83,  the  value  of  vertical  stress  at 
yield  becomes 

Thus,  when  o^  reaches  the  value  given  by  Equation  U81;,  the  material 
yields  and  continued  application  of  vertical  stress  causes  the  material 
to  move  along  the  yield  surface,  undergoing  plastic  as  well  as  elastic 
strains.  In  the  plastic  range  Equation  li79a  no  longer  applies  and  re- 
course to  Equation  U79b  is  necessary.  According  to  Equation  U79b,  the 
deviator  stress  increment  dS.  in  the  plastic  range  is  given  by 

GS  dE 

dS^  = 2G  dE^ (U85) 

& 

The  rate  of  work  S dE  for  conditions  of  uniaxial  strain  reduces  to 

Tnyt  Tim 

= S,  dE,  + 2S-  dE-  (486) 

mn  mn  11  2 2 


Utilizing  the  fact  that  = dE^^  = 0 Equation  486  reduces  to 


S dE  =^S,dE- 
mn  mn  2 1 1 


(487) 


In  view  of  Equations  480b  and  487,  Equation  485  becomes 


(488) 


2 — 3 2 

In  the  plastic  range  k “ ^2  ~ ¥ ®1  Equation  488  reduces  to 


dSj^  = 0 


1489) 


Since  dSj^  ° 0 * Equation  489  indicates  that  dS.,  = 0 also,  and 


do^  ^ dS^  * dJ^/3  = dJ^/3  (490a1 

= dSg  ♦ dJ^/3  = dJj^/3  (490b) 
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Equation  U90  indicates  that  the  material  behaves  as  though  it  vere  a 
fluid  once  it  has  reached  its  limiting  shear  resistance.  From  Equa- 
tions U79b  and  U80a,  it  follows  that 


dJ^/3  = K de^  (U91) 

Substituting  Equation  U91  into  Equation  U90a,  the  vertical  stress-stredn 
Increment  relation  in  the  plastic  range  becomes 

do^  = K de^  (1*92) 

Thus,  the  loading  slope  of  the  versus  cxirve  breaks,  or  softens, 
when  yielding  occurs  and  becomes  equal  to  the  bulk  modulxis.  Accordixigly, 
the  loading  slopes  of  the  principal  stress  difference-pressure  curve  and 
the  principal  stress  difference-strain  difference  cvirve  become  zero. 

Since  de^  = 0 , the  slope  of  the  pressure-volumetric  strain  c’orve 
remains  constant.  Once  the  material  unloads,  it  behaves  as  a linear 
elastic  solid  again,  satisfying  Equation  U81.  If  unloading  is  continued 
until  the  lower  yield  surface  corresponding  to 


(o  - o ) = -t.  (I493) 


is  reached,  the  material  flows  plastically  again  and  Equation  U79b 
governs  the  behavior  of  the  material.  The  foregoing  analyses  are 
depicted  schematically  in  Figure  17.  From  Figure  17  it  can  readily  be 
seen  that  for  a Prandtl-Reuss  material,  the  vertical  stress-stredn 
curve  associated  with  uniaxial  strain  configuration  would  break  or 
soften  when  yielding  occurs  and  would  remain  concave  to  the  strain  axis 
with  continued  application  of  vertical  stress. 

137.  bet  us  next  examine  the  behavior  of  Prandtl-Reuss  material 
under  a plane  stress  condition  defined  by 


0 0 
0 0 0 

0 0 o 


(U9U) 
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Figure  17.  Behavior  of  Prandtl-Reuss  material  under 

conditions  of  uniaxia: 

L strain 

For  this  state  of  stress 


1/2^2  \ 

= 3 V^l  ^ <^3  - 


(h95) 


Thus,  the  material  will  yield  when 


2 2 2 
<^1  ^ ^^3  - ‘^1^3  “ 


(496) 


Equation  496  describes  an  ellipse  in  the  coordinate  system 

(Figure  l8).  We  will  consider  a stress  path  where  is  held  constant 
at  k while  cf^  is  increased.  At  the  start  of  the  test  (assume  a 
compression  test),  point  A in  Figure  l8,  = 0 According  to  Equa- 

tion 496  the  material  yields  when  = 2k  , point  B in  Figure  l3. 

Prior  to  yield  the  behavior  of  the  material  is  governed  by  Eouation  479a 


AB_AN0  AC  ARE  STRESS  PATHS  VJk  I 

de'*  - PLASTIC  STRAIN  INCREMENT  VECTOR  H 

Figure  18.  Von  Mises  yield  curve  for  special  plane  stress  condition 
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de^  = E de^ 
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(it9Ta) 

(h9rb) 

(h9Tc) 


At  point  B the  material  yields  and  it  follows  from  Equation  1+59  that 

dfiP  = 0 (U98a) 

de|  = -deP  (l+98h) 


Unlimited  plastic  deformation  takes  place  at  yield.  It  is  noted  from 
Equation  ^+98  that,  as  expected,  = 0 . If  we  now  repeat  the  same 

test  and  change  the  direction  of  (i.e.,  a tension  test),  we  find 

that  the  material  yields  when  = -k  , point  C in  Figtire  l8.  At 
point  C the  material  yields  in  tension  and  from  Equation  U59  it  follows 
that 

de|  = 0 (li99a) 

dej  = -de^  (l*99b) 


The  concept  of  normality  can  he  demonstrated  from  this  simple  example  by 
superimposing  the  plastic  strain  coordinates  on  the  stress  coordinates 
; n Figure  I8.  As  shown  in  Figure  I8,  in  the  case  of  the  compression 
test  de^  *=  0 and  the  plastic  strain  increment  vector  de^  is  per- 
pendicular to  the  yield  surface  at  point  B.  In  the  case  of  the  tension 
test,  on  the  other  hand,  de^  = -de^  indicating  that  the  plastic  strain 
incrsnent  vector  is  perpendicular  to  the  yield  surface  at  point  C. 
Drucker-Prager  material 

138.  The  Von  Mises  yield  condition  was  modified  by  Drucker  and 

12 

Prager  to  include  the  effects  of  the  hydrostatic  stress  on  the  shear- 
ing resistance  of  the  material . The  yield  function  ^ was  assumed  to 
tedte  the  following  form 


i = k 


(500) 
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where  , a positive  material  constant,  represents  the  frictional 
strength  of  the  material.  Equation  500  describes  a right  circular  cone 
in  the  principal  stress  space  (Figure  19).  Substituting  Equation  500 


Figure  19.  Drucker-Prager  yield  surface  in 
principal  stress  space 


into  Equation  U69  we  obtain  the  following  stress-strain  relationship 
associated  with  the  Drucker-Prager  yield  function 


— M 


' HM 
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■m 


* 1 
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S dE  - 3Ka . dl^ 
mn  isn  f 1 


9KaJ;  + G 


(502) 


indicating  that  for  Drucker-Pre^er  material,  as  a consequence  of  depen- 
dency of  yield  function  on  hydrostatic  stress,  plastic  deformation  is 
accompanied  hy  volime  expansion  (it  is  noted  from  Equation  502  that  if 
= 0 the  plastic  volumetric  strain  is  zero  also).  The  increment  of 
total  volumetric  strain  dl^  can  he  determined  from  Equations  501  and 
500.  From  Equation  501  we  have 


d£  - J.  dlw'3)  - 3Ka  dl 
—f  \ mn  mn  11/  f 1 


“i  = 3r 


VjT 


(503) 


9Ko^  + G 


Solving  for  dl^  and  cons! lering  the  fact  that  during  plastic  deforma- 
tion vn  - a.J,  = k (Equation  500),  we  obtain 


VjJdJl/  2 ^ 30 

dl,  = — zy.-’  (gKtt^  + G)  - a de 

•1  3KGk  \ f / k mn  mn 


(50lt) 


The  increment  of  plastic  volumetric  strain  de^  then  becomes 


Vdp  / p \ 


0 de  - 
mn  mn  3K 


- (505) 


It  should  be  pointed  out  that  the  volume  change  is  due  to  scalar  non- 
linearity and  represents  uniform  dilatation.  For  example,  consider  a 
simple  shearing  stress  defined  by  the  following  stress  increment  tensor 


(507) 


d"l  = deg  = <ie3 


That  is,  there  are  no  normal  deviatoric  strains  associated  vith  the 
volume  expansion. 

139.  Substituting  Equation  500  in  Equation  U70  we  obtain  the 
following  relationship  for  the  stress  increment  tensor  for  the  Drucker- 
Pi’ager  material 


do.  = 2G  dE.,  + K dl,  6.,  - 
i.3  ij  1 ij 


dE„^  - 3Ka . dl, 
mn  mn  f 1 


9KaJ  + G 


S.  . - 3Ka-6., 
ij  f iJ 


(508) 


Equation  508  (or  Equation  501)  governs  the  behavior  of  Drucker-Prager 
material.  The  effect  of  the  dependency  of  the  yield  function  on  hydro- 
static stress  can  be  further  demonstrated  by  examining  the  behavior  of 
Drucker-Prager  material  under  unia:tial  state  of  strain  (Eqxxaticn  U80). 
The  elastic  behavior  of  the  material  is  given  by  Equation  48l.  The  mate- 
rial yields  when 


^ (0^  - Oj)  - a^(o^  . 20j)  = k 


(509) 


In  view  of  Equations  U8I  and  509,  the  value  of  vertical  stress  0^  at 
yield  becomes 

iilajLisi  k = — (510) 

^ 6g  - 9'^  Ka^  2G  - 3/3  Ka^ 


It  is  noted  that  if  is  set  to  zero.  Equation  510  reduces  to  Equa- 
tion U8U,  the  corresponding  expression  for  Prandtl-Reuss  material.  The 
effect  of  in  this  case  is  to  increase  the  value  of  the  vertical 
stress  0^  at  yield.  Wlien  reaches  the  value  given  by  Equation  510, 
the  material  yields.  Continued  application  of  vertical  stress  causes 
the  material  to  move  along  the  yield  surface,  undergoing  both  elastic 
and  plastic  deformation.  From  Equation  508  the  incremental  relation  be- 
' 'oen  vertical  stress  and  vertical  strain  becomes 


I 

i 


(511) 


Again  it  is  noted  that  when  is  set  to  zero  Equation  511  reduces  to 
Equation  H92,  the  corresponding  expression  for  Prandtl-Reuss  material. 

As  was  pointed  out  previously,  for  Drucker-Prager  material  plastic  de- 
formation is  accompanied  hy  volume  expansion  (see  Equation  502).  Accord- 
ingly, using  Equations  50^  and  509,  we  obtain  the  following  incremental 
relation  for  volumetric  strain  in  the  case  of  uniaxial  strain  test 


9Ko, 


dJ3_  = 


° ) 


9Ka^  + G 


dl^  + 3K  dl^ 


(512) 


When  is  set  to  zero.  Equation  512  reduces  to  the  corresponding 
expression  for  elastic  material.  The  increment  of  plastic  volumetric 
strain  then  hecomes 


^ a„(9Kct.  - 2/3  G) 

deP  1 dJ.  (513) 

^ 3KG(1  + 2v^  a^) 


In  order  for  the  uniaxial  strain-stress  path  to  reach  the  yield  surface, 
the  following  condition  should  hold 

> 3«,  (Slk) 

Therefore,  as  expected,  the  increment  of  plastic  vol’imetric  strain  is 
negative  (expansion). 

Work-Hardening  Plastic  Material 


1^*0.  In  the  case  of  work-heu*dening  plastic  material,  the  yield 
surface  ^ is  not  fixed  hut  expands,  or  translates,  as  plastic  defor- 
mation takes  place.  The  material  can  then  sustain  stresses  beyond  those 
required  to  reach  the  initial  yield  condition.  Therefore,  ve  can  use  a 
loading  concept  in  the  case  of  work-hardening  plastic  material  according 
to  the  direction  of  the  stress  increment  tensor  do^j  (viewed  as  vector). 
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During  loading  from  a point  on  a given  yield  surface  the  stress  vector 
is  pointing  outvard  and  thus  (3jJ/9a^^)dcf^j  > 0 • During  unloading  the 
stress  vector  is  pointing  inward  and  thus  • Accord- 
ingly for  work-hardening  plastic  material  we  define 


_3jL 

3o. 

> 0 

loading 

(5158.) 

-iiL 

'"“ij 

< 0 

'mloadiag 

(515b) 

The  condition 

= 0 

(i.e. , 

when  da . , 
ij 

is  tangent  to  yield 

surface)  is  known  as  neutral  loading  and  produces  no  plastic  deformation 
in  the  case  of  work-hardening  material.  The  stability  condition  for 
work-hardening  plastic  material  is  given  as 


da . . de . . > 0 
ij  ij 

(516a) 

da . , de? . > 0 
ij  *- 

(516b) 

where,  xmlike  ideal  plastic  material,  the  equality  sign  in  Equa- 


tion 516b  holds  only  when  de?.  = 0 . For  work-hardening  plastic  mate- 
11 

rial  Drucker  has  shown  that  the  e:cpression  for  plastic  strain  incre- 


ment tensor  is  simileu:'  to  Equation  1*59  where  the  proportionality  factor 
A depends  on  stress,  plastic  deformation,  and  history  of  plastic  de- 
formation. We  can,  therefore,  use  Equation  i*60,  in  conjunction  with 
the  loading  conditions  given  by  Equation  515,  for  calculating  the  strain 
increment  tensor.  During  loading  from  a point  on  the  yield  surface 
(0^/3o  > 0),  Equation  k60  governs  the  behavior  of  the  material 

In  the  elastic  range,  and  during  unloading  from  a point  on  the  yield 
surface  ((3|5/3o^y  dc^  ^ < o),  the  behavior  of  the  material  is  governed  by 
Equation  U53.  wqien  ( ® (neutral  loading),  deV^  = 0 and 
Equations  U60  and  1*53  become  identical  (thus  establishing  continuity  at 
a load-unload  interface), 

lUl.  We  now  adopt  a yield  condition  of  tne  following  type 


(517) 
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for  strain-hardening  material.  Equation  517  indicates  that  the  yield 
surface  is  not  fixed  in  the  principal  space  and  that  it  changes  as 
plastic  deformation  takes  place.  We  further  assume  that  k is  a con- 
stant. Following  the  same  procedure  as  was  used  to  derive  an  expression 
for  A in  the  case  of  ideal  plastic  material  we  obtain 


A = 


3.1,  6G  3.^, 


_ii.  _y_  + 3K-JG 
3o,,  30.,  6G 
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Equation  518  is  the  espression  for  the  proportionality  factor  A asso- 
ciated with  the  strain-hardening  yield  condition  given  by  Equation  517. 

It  is  ‘noted  that  Equation  518  reduces  to  Equation  when  the  dependency 
of  the  yield  function  on  the  plastic  strain  disappears  (i.e.,  ideal 
plastic  material).  In  view  of  Equations  518  and  WO,  the  strain  incre- 
ment tensor  associated  with  the  yield  condition  of  Equation  517  becomes 


(519) 


We  can  also  dei'ive  an  expression  for  the  stress  increment  tensor 
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For  a mm'ber  of  engineering  materials,  soils  in  particular,  the  yield 


function  ^ is  expressed  in  terms  of  , , anc 


and  , i.e.. 


(521) 


For  the  above  specification  of  , Equations  519  and  520  become 
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It  is  noted  that  Equations  522  and  523  reduce  to  Equations  U69  and  UTO, 
respectively,  when  the  dependency  of  the  yield  function  on 
disappears. 

llt2.  In  order  to  demonstrate  the  application  of  Equation  522  (or 
Equation  523)  let  us  consider  an  elliptic  yield  function  defined  by  the 
following  equation  (Figure  20) 


<(•'1  • = 


(52U) 


Equation  52^*  has  been  used  successfully  for  modeling  the  stress-strain 
behavior  of  earth  materials. For  a first-order  approximation,  the 
variable  Y , which  controls  the  expansion  of  the  yield  surface,  is 
assumed  to  take  the  form 


Y = AeL 


(525) 
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lk3.  Let  ns  now  examine  the  behavior  of  Equation  527  under  l^dro- 
static  state  of  stress  (Figure  6a).  For  hydrostatic  state  of  stress 
= Y and  Equation  527  becomes 


(3K  * A) 


(528) 


It  should  be  noted  that  for  this  state  of  stress  the  same  results  could 
have  been  obtained  directly  from  Equations  1+52,  li53»  and  525  without 
recourse  to  Equation  527.  For  virgin  loading.  Equation  528  can  be  in- 
tegrated to  yield  (assuming  zero  initial  pressure  and  volumetric  streiin) 


J,  = K I, 
1 el 


During  purely  elastic  deformation  (Equation  1*53) 


(529) 


dJ^  = 3K  dl^  (530) 

Since  Kg  < 3K  (Equati'^n  528),  it  follows  that  plastic  compaction  pro- 
duces an  apparent  softening  of  the  effective  bulk  modulus.  Figure  21 
depicts  the  behavior  of  the  material  under  hydrostatic  state  of  stress. 
The  behavior  of  the  material  from  iJoint  1 to  point  2 is  governed  by 
Equation  529  (the  material  undergoes  plastic  as  well  as  elastic  deforma- 
tion). If  the  material  is  unloaded  froui  point  2 to  point  3,  and  then 
reloaded  from  point  3 to  point  2,  the  behavior  is  elastic  and  the  re- 
sponse of  the  material  is  governed  by  Equation  530. 

lUU.  Let  us  next  examine  the  behavior  of  Equation  527  under  a 
constant-pressure  shear  test  (Figure  6c).  The  queditative  behavior  of 
the  model  is  depicted  in  Figure  22.  The  material  is  first  hydro- 
statically loaded  from  point  1 to  point  2.  The  response  of  the  material 
from  point  1 to  point  2 is  governed  by  Equation  529  and  is  identical  to 
that  shown  in  Figure  21  (the  material  undergoes  both  piaBt--*<'  and  elastic 
deformation).  The  material  is  then  sheared  from  point  2 to  point  3 by 
increasing  while  is  kept  constant.  Since  is  kept  con- 

stant, all  volume  changes  from  point  2 to  point  3 are  plastic.  From 
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Figure  22  Behavior  of  work-hardening  elastic-plai’+ic 
material  during  constant-pressure  shear  test 

Equation  U59  it  follows  that  the  increment  of  plastic  volumetric  strain 
is  given  as 

de^  ~ 3A  (531) 

In  view  of  Equation  52U,  Equation  531  hecomes 


= 3A(2J^  - Y) 


(532) 


Since  A is  iKjsitive,  Equation  532  indicates  that  the  plastic  volumet- 
ric strain  during  the  shearing  process  is  positive  (compaction).  At 
point  3 ® (normality  condition)  and  the  yield  surface  ceases  tc 


expand.  The  shearing  response  of  material « expressed  in  terns  of 
versus  Vf  t then  reaches  its  icaxlmun  value  (for  the  particular  value 


of  at  point  2)  asymptotically  at  point  3*  As  shown  hy  the  dashed 
lines  in  Figure  22,  if  the  material  were  to  linload.  from  any  point  dizring 
the  shearing  process  it  will  behave  as  a linear  elastic  material.  This 
single  example  points  out  the  basic  difference  between  ideal  and  work- 
hardening  plastic  materials.  That  is,  for  ideal  plsistic  materials  the 
yield  surface  is  fixed  and  does  not  e:q)and  during  plastic  deformation. 
Unlimited  plastic  flow  takes  place  at  the  onset  of  yielding.  In  the 
csise  of  work-hardening  material,  on  the  other  heind,  the  yield  surface 
moves,  or  expands,  causing  the  material  to  harden  as  plastic  deformation 
takes  place. 
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